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This thesis describes a new finite element 
structural analysis computer program "ASTRA" developed 
by the author for use on the APPLE II microcomputer. 

The program is written in Applesoft Basic and is based 
on the displacement method of matrix structural 
analysis. It has a capacity of twenty five nodes and 
elements, and two types of structural elements: a space 
rod and a space beam. Along with a review of the basic 
applicable structural theorems needed, a thorough 
discussion of the program is presented. Description of 
special purpose subroutines, methods of data handling/ 
storage, and suggestions on possible improvement are 
also included. It is hoped that this program can serve 
as a base from which more advanced programs can be 
developed for use with the APPLE II microcomputer. 
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Chapter 1 


INTRODUCTION 

The rapid advances of science in the field of . 
data processing, particularly the introduction of 
microcomputers, has made available low cost computer 
equipment to engineers and scientists. These computers 
can greatly enhance and speed the solution of complex 
engineering problems. 

In the field of structural analysis, the idea of 
utilizing computers to obtain the solution of complex 
structural problems is not new, but the majority of 
modern finite elements structural analysis programs are 
available only on large and expensive computer systems. 
It is felt that modern microcomputers are powerful 
enough to warrant the creation of software for the 
analysis of structures using the finite element method. 
This study is intended to validate this concept by 
developing a finite element structural analysis program 
for use on the APPLE II Plus microcomputer. The program 
named ASTRA, an abbreviation for "Analyzing Structures 
with APPLE," is not intended to compete with the more 
advanced structural programs available on main-frame 
computers. Rather, it has been designed to provide an 
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alternative program for use by small companies which do 
not have the capital to employ one of the larger 
programs, by field engineers to use in remote locations 
where access to main-frame computers is limited, or by 
educational institutions teaching the finite element 
method. 

ASTRA is made up the pre-processor IASTRA (Input 
to ASTRA) and the main program ASTRA. IASTRA was 
written to obtain all the information required for the 
analysis of structures and to generate the input data 
for use in ASTRA. In this manner, the amount of time 
needed to train a future user of the program should be 
minimized and errors reduced. Extra care was also 
devoted to make the program ASTRA user oriented, i.e., 
the input of information necessary for the running of 
the program was made as easy as possible in order to 
make the program convenient to use with little formal 
training. 



Chapter 2 


BASIC STRUCTURAL THEOREMS 

Introduction 

The methods of analysis discussed here are 
limited in application to elastic structures. The 
following assumptions are made for each element: 

1. The body is elastic 

2. The material properties of the body are 
homogeneous and isotropic 

3. The deformations are assumed to be 
infinitesimal 

The Theory of Elasticity forms the basis for 
performing analysis of structures. In order to define 
the various terms involved in the displacement and force 
methods of analysis, the fundamentals of the Theory of 
Elasticity and basic structural theorems will first be 
quickly reviewed. 

Three conditions must always be met in any 
analysis of the internal forces and deformation in a 
structure, these are: 
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1. Equilibrium of forces must exist 

2. Laws of material behavior must be obeyed, and 

3. The displacement results must be compatible 

The first condition requires that the internal 
forces balance the applied external loads. This 
condition alone is sufficient to enable the solution of 
statically determinate proolems, however, it yields 
insufficient information to enable the analysis of 
redundant structures. Under those circumstances the 
laws of material behavior which in problems of linear 
elasticity, neglecting temperature changes, reduce to 
Hooke's Law, and the conditions of compatioility must be 
invoked in order to complete the analysis. 

EQUATIONS OF EQUILIBRIUM 

If we consider the cube of Figure 1 and express 
the equilibrium of forces acting on it, in the x 
direction one has: 

Xdxdydz — a x dydz + -f dx j dydz — z, x dzdx 

4- (r yx + ^dy^jdxdz - t lx dydx + + ^dz^jdydx = 0 


( 2 . 1 ) 


5 

where X is the force in the x direction per unit 
volume. Cancelling dx ay az we obtain: 


- 1 - 4-^=0 

ox ay oz 


( 2 . 2 ) 


Similar equations can be obtained for the equilibrium of 
forces in the y and z directions. 


'&+&+&+ 5 '- 0 


(2.3) 
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+ dy 


^i + Z = 0 

az 


(2.4) 


Equations (2.2) through (2.4) are Navier's equations of 
equilibrium for an elastic solid. 
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Stress-Strain Relationships 

The relationship between stress and strain, or 
the law of material behavior for a linear elastic body 
below its yield point are given as: 

^fox “t" ^i)! 

~ v (tf, + OJ (2.5) 

f, = -f <7,)] 


and 



( 2 . 6 ) 


I 
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Equations of Compatibility 

The components of strain in terms of the 
components of elastic displacement can be presented as 
follows: 




du 

~Tx 




(2.7) 


and 


_ du . dv 
- Ty + Jx 

v _dv , 

7 " “ Tz + 

_ dw du 

7 ” ~ Tx + Tz 


( 2 . 8 ) 
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The components of strain are not independent but 
are related by relationships called compatibility 
equations. From Equations (2.7) and (2.8) it follows 
that 


, d*e, = d 2 y x , 

dy 2 dx 2 6xdy 


5z 2 ^ 5v 2 


dydz 


d'e, 

IF 


= 5 i y„ 

dz 1 dzdx 


(2.9) 


and 


5 

5J 


d y*> #y„ dy» 
FT Fx ~ Fv 


= 2 


d 2 e x 

575r 


d 

Ty 

d 

Tz 


d y» dy„ , dy„ 

JT - Ff + ~FI 

dy„ _ dy„ 5>’„1 

' FF] 


FT~FT 


- 2 dlf r 

~"57c£ 

5x57' 


( 2 . 10 ) 
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Substituting Equation (2.5) into Equations (2.8) and 
(2.9) the Beltrami-Michell compatibility equations are 
obtained as follow. Let © = a x + a, + a, and A 
represent the Laplacian operator, then 


Ac x -+- 

1 

_ <? J © _ 

v I 

<dX , 

dY 

dZ\ 

-2 d * 

1 + 

v dx 2 

1 - v' 

\&x 

Ty 

Ur/ 

l Tx 

Act, + 

1 

d 2 Q _ 


f dX , 

dY 

dZ\ 

-2 d -X 

1 + 

ri\ 

% 

1 

1 - v' 

\tx 

Jy ‘ T 

o: ) 

U y 

Ac, 4- 

1 

d 2 Q _ 

. v __ | 

<dX . 

dY 

dZ\ 

-2 d J 

1 + 

v d: 2 

1 - v' 


Ty 

Uz/ 

Ur 


( 2 , 


For a more extensive derivation see Volterra and 
Gaines (1971). 
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Chapter 3 


FINITE ELEMENT METHOD OF STRUCTURAL ANALYSIS 

Introduction 

The finite element method of structural analysis 
extends the matrix methods of structural analysis for 
application on high speed digital computers. It is 
based on the concept of replacing a continuous structure 
by an equivalent mathematical model made up of 
structural elements having known material and geometric 
properties expressed in matrix form. These individual 
matrices are then assembled following a set of rules 
derived from the Theory of Elasticity, to provide the 
geometric and material properties and static and dynamic 
characteristics of the structure. 

The basic aims of structural analysis is to 

ascertain the stresses and deformations of a structure 

as it is subject to varied conditions of loading, 

temperature changes, applied deformations, etc.. As the 

structures become more complex, the finite element 

method is the only alternative to performing an accurate 

structural analysis. The most common matrix methods 

utilized to obtain the approximate solution to linearly 

elastic structures are the unit load method, the unit 
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displacement method, and the basic variational form of 
the displacement method derived from minimum potential 
energy considerations. The concise and systematic 
notations of matrix algebra utilized in these methods 
are ideally suited for programming digital computers to 
solve the large number of equations generated in the 
analysis of highly redundant structures. 

Matrix methods of structural analysis are 
classified as matrix force or matrix displacement 
methods depending on whether forces or displacements are 
the problem coordinates. These two methods will be 
presented in detail in the following sections. 

One dimensional elements (rods) have been chosen 
here to illustrate the two methods because exact 
solutions exist for the one dimensional problem. In 
addition, they are particularly useful because they 
allow the various steps of the solution process to be 
clearly examined. 

Structure Displacement 

In the finite element method the structure is 
subdivided into separate substructures known as elements 
which are joined and loaded at a number of points or 
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nodes. Each of these nodes can have up to six 
components of joint displacements and rotations 
(u, v, w, y, fi, a) called degrees of freedom. This 
idealized structure is known as a mesh or mathematical 
model. 

The notations used for the force and 
displacement methods used in the following sections are 
now introduced. [Q] will be the column vector that 
describes the internal forces while [S] will be a column 
vector describing the corresponding internal 
displacements. A subscript will be used to show the 
member and direction [Q] and [S] belong to. For 
example, [Q n ] indicates the internal forces for the Nth 
degree of freedom which in turn denote a particular 
member and direction of the internal force. If we let 
[Q] and [S] denote respectively the nodal forces and 
nodal displacements, then the external work done and the 
internal strain energy may be expressed as 

[We] = [Q] [S] (3.1) 


and 


[Wi] = [Q] [S] 


(3.2) 
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Equating [We] and [Wi] we obtain 

IQ] [S'] = [Q] [S] (3.3) 

For a statically determinate structure, the 
member forces and deflections may be expressed in terms 
of the external nodal loads 

IQ] = [RJ IQ] (3.4) 

[S] = [R] [S] (3.5) 

where [R] is a matrix that transforms the external 
forces to the elements forces. For rod and beam 
elements, it has the physical meaning of a rotation 
matrix that maps the global or structure system into the 
local or element coordinate systems. 




fll 

ml 

nl 

[R] 

= 

: 12 

m2 

n2 



[13 

m3 

J 


(3.6) 
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where the origins of X Y Z (global) and X' Y' Z' (local) 
systems are the same and 11, ml, nl, 12, m2, n2, 13, m3, 
n3 are the direction cosines of the X' Y' Z' axes 
relative to the X Y Z axes, respectively. 

The Matrix Force Method 


In the matrix force method of structural 
analysis the internal loads and external reactions are 
the unknowns. The correct system of loads is that which 
satisfies the minimum energy condition. 

The flexibility coefficient constitutes a 
relationship between the elements displacements and 
forces on a structure and are calculated from the 
geometric and material properties of the elements. This 
coefficient is an expression of the reduced structure 
and each element is defined as follow: 


/ 


k 1 
i j 


is the displacement necessary at node i acting 
in the k direction required to sustain a unit 
force at node j in the 1 direction 


The element flexibility coefficients can be 
determined from deformation geometry, virtual work, or 
strain energy methods. The element flexibility 
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matrix [F] is assembled from the corresponding element 
flexibility coefficients. For example, for a plane rod 
element with two degrees of freedom per node, i.e., two 
translations in the X and Y directions at each end of 
the element, the flexibility matrix is: 


XX 

f 11 

XV 

*ii 

XX 
f 12 

f Xy 
f 12 

f n 

f 12 

f 13 

f 14 

f yx 
*11 

f y y 

r n 

f Yx 

M2 

1 12 

j 

; f 21 

f 22 

f 2 3 

f 24 

1 

rXX 
f 21 

f 21 

rXX 

t 22 

f xy 

t 22 

i 

| f 3 1 

| 

f 32 

f 33 

1 

! 

f 34 

_y x 

i f o i 

U " 1 

f yY 

f 21 

yx 
f o o 

yy 

f 2 2 

| 

Lf 41 

f 42 

f 43 

f 4 4_j 


The matrix [F] in Equation (3.7) is shown using 
two different notations to identify the individual 
matrix components. A little experimentation will show 
that the use of two subscripts/superscripts to identify 
the element flexibility coefficients is impractical when 
dealing with multi-element structures or writing 
computer codes. For this reason, degrees of freedom 
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numbers are used as subscripts here in this thesis to 
explain the theory of the force and displacement methods 
of analyses and in the eventual programming of the ASTRA 
computer code. 

The structure flexibility matrix [^] is 
assembled by the addition of all applicable individual 
element flexibility matrices [F] that are applicable. 


[F s ] = ± [R-j^ [F] (l) [R] (l) (3.8) 

i = 1 

Here the total number of elements in the structure is n, 
and [R] is the force transformation matrix. 

The deformation of a member having n degrees of 
freedom (DOF) may be expressed using the flexibility 
coefficients in terms of the separate influences of the 
set of member forces [Q]. 


T1 


S 1 


*11 

f 12 

f 1 3 * 

• 

• 

f ln ! 

j 

1 Q 

i ■ 

s 2 


i f 2 1 

[ 

f 2 2 

f 2 3 * 

• 

• 

f 2n 

Q; 

S 3 


|*31 

1 

f 3 2 

f 3 3 * 

• 

• 

f 3n | 

Q; 

• 

— 

1 

! • 

• 

• • 

• 

• 

• 

• 


^nl ^n2 ^n3 


nn 


► | 


\ i * 

1 


Qn 


(3.9) 
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Equation (3.9) can be more concisely expressed in the 
local or member coordinate system as 

[S] = [F] [Q] (3.10) 

and in the structure or global coordinate systems 

[R] [S] = [F] [R] [Q] (3.11) 

The relationsnip between the nodal displacements 
[S] and the external forces [Q - ] may be derived by 
expanding Equation (3.11) to include more than one 
element as follows: 

[S] = [R] T [F] [R] [Q] (3.12) 

[S] = [Fs] [Q] (3.13) 

where [Fs] is the structure flexibility matrix and can 
be found using Equation (3.8). 

For statically determinate structures, Equation 
(3.13) gives a direct solution for all the nodal 
displacements in terms of the external forces acting at 


the nodes. 
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If the structure is statistically indeterminate, 
then the work performed by the redundant forces must be 
included and the compatibility conditions invoked. For 
these cases, Equation (3.13) becomes 


[S] = [F] [Q] - [S r ] 


(3.14) 


where [Sr] are the displacements of the primary 
structure at the points of redundancy. The conditions 
of compatibility state the the displacements at all the 
releases of redundant points caused by the applied loads 
and redundant forces must be made to vanish. 

Equation (3.14) may be expressed in matrix form 
as 

(3.15) 


s s l 


i F 




ss 


r s 


sr 


r r 



Q s~ 


5r 


where [Q £ ], [S„] stand for the structure and [Q^], 

[S3 for the redundant forces and deflections. The 
compatibility condition is expressed as 

l F ts ] [Q S J + [F rr ] [Q r ] = [S r ] 

A 
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however, if the redundant coordinates are fixed in 
space, [S r ] = 0 and 

[F rs 3 [Q £ ] + [F rr ] [Q r ] = 0 (3.17) 

from which we obtain 

IQ r l =-lF rr [ 1 [F rs ] IQ S ] (3.18) 

which describes the solution for the redundant loads. 

The solution procedure for redundant structures 

then is: 

1. Define the external loads [Q] and specify 
the redundancies. 

2. Calculate the force transformation matrix 
[Rl and the element transformation matrix [F] for all 
the elements. 

3. Assemble the structure flexibility matrix 

4. Solve for the redundant forces [Q r ]. 

5. Calculate the nodal displacements [S ]. 

6. Solve for the member forces [Q] . 

For a more extensive derivation of the force 
method see Yan-Yu Hsieh (1970) . 
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The Displacement Method 

In the displacement method of structural 
analysis, the structures displacements are the basic 
unknowns. A set of equations of equilibrium equal to 
the degree of kinematic indeterminacy (number of unknown 
diplacements) has to be solved in order to determine 
these displacements. 

In this method, the compatibility conditions are 
first satisfied by correlating the external nodal 
displacements to the displacements of its members. Then 
a force-displacements relationship is established 
between the member's end forces and deformations. 
Finally, the equilibrium equations are used to calculate 
the nodal displacements, member forces, and reactions of 
the structure. 

The compatibility conditions used in the 
displacement method express the condition that the 
member deformations [S] must be consistently related to 
the nodal displacements [S]. Let represent the 

value of the member deformations caused by a unit 
nodal displacement Sj. The value of the total member 
deformations [S] caused by all the nodal displacements 
[S] may be written as 


21 
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R 11 

R 12 

R 13 
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In 
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R 21 
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R 23 
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R 31 

| 

R 32 

R 33 

* * ‘ R 3n ! 
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! 
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t 

| • 
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R nl R n2 R n3 


R. 


nn 


r*x] 

! s„ 1 


i ! 


n 


(3.19) 


or 


[S] = [R] [S] 


(3.5) 


and similarly we have 


[Q] = [R] [Q] 


(3.4) 


where [Q] is the vector of element forces and [Q] the 
vector of nodal forces. Tne matrix [R] is the 
displacement transformation matrix and relates the 
internal member deformations to the external nodal 
displacements. For beam and rod elements it can be 
physically represented as a geometric transformation of 
coordinates from the structure coordinate system'to the 
element coordinate system as mentioned previously. 
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For all structural elements a direct 
relationship exists between the forces applied at the 
boundaries and the boundary displacements. This 
relationship can be expressed in matrix notation as 

[Q] = [K] [S] (3.20) 

where [Q] is the vector of element loads, [S] is a 
vector of element displacements and [Kj is the element 
stiffness and is defined as follows: 


is the force at node i acting in the m direction 

required to sustain a unit displacement at node 

i J 

j in the 1 direction 


and, for example, in the case of a two dimensional rod 


[K] 


. xx , xy xx . xY 
k ll k ll k 12 k 12 

, yx . yy . yx yy 

k l 1 k l1 k 12 k 1 2 


1. 


! k 11 k 1 2 k 1 3 k 14 
k 2 1 k 22 k 2 3 k 24 


1 


k XX k X ^ k XX k X ^ k V k k (3.21) 

! k 21 k 21 k 22 K 22! K 31 K 32 K 33 K 34 I 


I ..y* k yy w v x k yy; 

i*^21 k 21 ^22 K 22i 


^41 k 42 k 43 k 44J 


Tne element stiffness matrix [K] and the structure 
stiffness matrix [Ks] are assembled by the addition of 
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all the individual stiffness coefficients and element's 
stiffness matrices respectively in a manner similar to 
that employed to asssemble the element and structure 
flexibility matrices. 

For a structure composed of one element having n 
degrees of freedom (DOF) we have 



(3.22) 


which expressed in abbreviated form becomes 


IQ] = U] [S] 


(3.20) 
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For an idealized structure composed of multiple 
members, the nodal forces and nodal displacements can be 
expressed in terms of the element forces and 
displacements by substitution of Equation (3.4) and 
(3.5) into Equation (3.20) to obtain 

[R] IQ] = [K] [R] [S] (3.23) 

or 

[Q] = [R] T [K] [R] [S] 

which reduce to 

[Q] = [Ks] [S] (3.24) 

where [Ks] is the structure stiffness and is assembled 
from the sum of all the individual element stiffness 
matrices expressed in the global or structure coordinate 
system. 

[Ks] T(l) [K] (1) [R] (1) (3.25) 

i = i 

The solution procedure of the displacement 
method is straight forward and can be summarized as 


follows: 
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1. Compile the basic data of the structure 
(material properties, geometry, load conditions, etc.) 

2. Construct the idealized structure; define 
the location of nodes, element orientation (topology), 
element material and geometric characteristics. 

3. Calculate the element stiffness [K]. 

4. Calculate the element transformation 
matrices [R]. 

5. Calculate the structure stiffness 

[KS] = £ [R] T(l ' [K] (l) [R] (l) (3.2 

i = l 

6. Eliminate degrees of freedom (DOF) not used 
to establish the reduced structures stiffness [Kr], 
loads [Q ], and displacement [S r ] matrices. Loads 

and displacement matrices for the redundancies are 
related by 


[Q r ] = [K r ] [S r ] 


(3.26) 
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7. Partition all matrices. Partitioning of 
matrices is used to lessen the work in matrix 
multiplication and inversion. 



(3.27) 



(3.28) 



(3.29) 


Substituting Equation (3.27) through (3.29) into 
Equation (3.26) 



K a/J 



(3.30) 
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Note that 


[K a * ] = [K^ a ] T 

and [S a ] is the vector representing the degrees of 
freedom which are free to move or displace. [S^ ] is 

the vector of degrees of freedom that are fixed in 
space. If the loads are expressed as applied forces and 
not as given displacement [S^ ] =0. Similarly, [Q a ] 
and [Q^ ] denote the nodal loads and reactions, 
respectively. 

8. Calculate the unknown displacements 

[K a a i 1 [Q 0 ] = [S a ) (3.31) 

9. Calculate the reactions 

IQ ] = IK^ ] [S a ] (3.32) 

10. Calculate the element forces and stresses 

[Q] (l) = [K] (i) [R] (l) [S 0 ] (3.33) 

Examples of structural problems using the 
displacement method are shown in Appendix E. 
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Comparison of The Force and Displacement Methods 

When comparing the force and the displacement 
methods of analysis, it must be emphasized that the same 
input information and element properties can be used 
regardless of the method to be employed in the 
analysis. Both methods lead to the same theoretical 
result, one being the inverse of the other. The 
computational path, however, leading to the calculation 
of displacements is different in each method. Because 
of possible ill-conditioning of the equations, different 
rounding-off errors, and the actual numerical analysis, 
the numerical results may differ for the two methods. 

If the matrix operations for both methods are 
compared, it is found that in the force method, through 
the use of the Jordanian elimination technique 
(Przemieniecki, 1968), the sequence of matrix operations 
is considerably more complicated than that for the 
displacement method. When writing computer programs for 
the solution of structural problems, except for 
structures which involve many joint displacements and 
few redundants, displacement methods or a combination of 
displacement and force methods are usually preferred. 

The displacement method uses the same procedures for 
analyzing statically determinate and statically 
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indeterminate structures, operates directly on the 
complete structure, and does not require a reduced 
structure as does the force method. In addition, in the 
displacement method each unit load is applied to each 
degree of freedom while all others are held zero. The 
equations also tend to be better conditioned, i.e., the 
largest terms lie on the main diagonal. In the force 
method a well conditioned flexibility matrix depends 
upon a good choice of redundants. Finally, for 
vibration problems the stiffness matrix [Ks] of the 
displacement method is determined directly and is better 
conditioned and easier to work with than the flexibility 
matrix [Fs] of the force method. 


Chapter 4 


THE COMPUTER PROGRAM "ASTRA" 

Introduction 

The computer program "Analyzing Structures with 
the APPLE" or "ASTRA" for short, is a finite element 
program based on the displacement method using the 
APPLE II or a similar microcomputer. The intented range 
of applications of the program extends to the solution 
of static structural problems in a plane, and in space, 
using a combination of beam and truss members with 
concentrated point loads and moments at the joints. 

The program has been written in such a way so as 
to allow conversion to microcomputer systems other than 
the APPLE II or to more advanced programming languages. 
Eventual expansion to include more advanced element 
types, and increases in the number of elements and nodes 
that it can handle has also been facilitated. This has 
been accomplished by the use of subroutines and comment 
statement. Although the use of comments involves a 
substantial use of available RAM, their inclusion in the 
program makes it much easier to understand and to follow 
the course of the matrix operations and use of the 
subroutines. 


30 



31 


When working with microcomputers, the programmer 
must make a decision as to how to best employ the amount 
of RAM and disk storage available. In tailoring ASTRA 
to the APPLE computer system, it was found that although 
program running time is critical to the efficient use of 
the computer, the limiting factor was the amount of 
memory available. Thus, the heart of of the ASTRA 
program is the efficient use of available disk storage, 
and the use of special purpose subroutines to speed up 
matrix operations. This is accomplished primarily in 
the calculation, storage, and retrieval of the structure 
stiffness matrix [Ks] and the partitioned stiffness 
matrices [K aa ] ana [K a ^ ]. Only the lower half of the 
bandwidth below the matrix diagonal is stored. Although 
this feature increases the running time and complexity 
of the program, it allows ASTRA to be run on a 
microcomputer by using the disk drive as the main 
storage device. 

While developing ASTRA, maximum thought was 
given to making the computer program "user oriented", 
i.e., a program that would be easy to use by anyone 
having some background in structural analysis with a 
minimum amount of training. This was accomplished by 
the use of the pre-processor "IASTRA" or "Input to 
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ASTRA." All that the user has to do to prepare or 
change input information required to run ASTRA is to 
answer the questions formulated by the computer. It is 
felt that this approach eliminates many of the problems 
associated with running a computer program and thus 
frees the user to concentrate on the analysis of the 
structure. 

The first version of the computer program ASTRA 
is the result of two years of research and experience 
with other finite element computer programs commercially 
available on large computers. It is hoped that it will 
become a flexible and efficient tool for the analysis of 
structures. The program presently contains the 
following types of elements: 

a) Three-dimensional truss element (rod) 

b) Three-dimensional bar element (beam) 

The structure to be statically analyzed may be 
composed of a combination of a number of these elements. 

The capacity of the program depends on the type 
and number of storage devices used, computer language, 
and time available for the analysis. Today, it is felt 
that a comfortable number of elements capable of being 
handled by the program is twenty-five. This is limited 
by the use of an APPLE II equipped with 48K of RAM and 




33 


116K mini-floppy disk. However, should a hard disk 
drive (10 megabytes) and an APPLE III with Pascal and 
124K RAM become available, the capacity of the program 
could be tremendously improved and the running time 
reduced by using matrix bandwidth storage techniques. 

The purpose of this section is to present the 
concept and the theory of the program ASTRA, and show 
the formulation of the element stiffness and 
transformation matrices and the solution procedures. 

Method of Analysis 

ASTRA is a finite element program for the 
solution of static problems based on the displacement 
method of structural analysis. The finite element 
method and general equations which govern the 
equilibrium of the system are given in Chapter 3. 
However, for completeness, the equations and solutions 
procedures are reviewed in the following discussion. 

The relationship between the forces applied at 
the boundary and the displacements in a structural 
element can be expressed as follows: 


[Q] = [K] [S] 


(3.20) 
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where [Q] is the vector of elements loads, [K] is the 
element stiffness and [S] is the vector of element 
displacements. The member deformations [S] and forces 
[Q] can be related to the nodal displacements [S] and 
nodal forces ["q] . This can be expressed in matrix 
notation as 


[Q] = [R] [Q] 


(3.4) 


[SJ = [R] [S] 


(3.5) 


The same relation can be defined globally for 
the complete structure. 

[Q] = [Ks] [SJ (3.24) 


where [Ks] is the structure,stiffness matrix as defined 
by Equation (3.25). 


= £ [H]™ [K] (i) [R] (i) 

i = i 


[KS] 


(3.25) 
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After the model is constructed, the element 
transformation and stiffness matrices are calculated and 
the structure stiffness matrix is assembled. Then, as 
expressed in Equation (3.26), the redundant 
displacements are eliminated to establish the reduced 
structure stiffness, displacements and loads 

fQr] = [Kr] [Sr] (3.26) 

and consequently partitioned per Equation (3.27) through 
(3.29). 


M = 



(3.27) 



K aa 

Ho Hfi 


(3.28) 




S 0 

[ s 3 

ii 






36 

Substituting Equation (3.27) through (3.29) into 
Equation (3.26) yields 


V 


K a a 

K«/ 






(3.30) 


If there are no induced displacement in the 
structure [S ^ ] = 0 and the unknown displacements, 
reactions and element forces are calculated 

(3.31) 

(3.32) 

(3.33) 

For a more complete derivation of the displacement 
methods see Hsieh (1970). 

Structural Elements 


[K aa ] _1 [Q a ] = [S a ] 

IQ 0 ) = [K/r a ] ISa ! 
[Q] (i) = [K] (i) [R] W [£„] 


The present version of the computer program 
ASTRA contains two structural elements. These are 
three-dimensional rods and beam elements. These two 
element types were selected because their stiffness and 
displacement transformation matrices are easy to 
calculate and exact solutions exist to check the final 





37 


results. In addition, their use enables the various 
steps of the solution procedure to be examined and it 
simplifies the programmer's task of debugging the 
computer program. 

The element description, characteristics, and 
calculation of the element stiffness and displacement 
transformation matrices are presented in the following 
sections. For a more detailed discussion see Hsieh 
(1970), and Przemieniecki (1968). 

The rod element shown in Figure 2 is a simple 
tension-compression bar element pin jointed at the ends 
and with uniform cros.s-sectional area. The beam element 
illustrated in Figure 3 is assumed to be a straight bar 
of uniform cross-section capable of resisting axial 
forces, bending moments about the two principal axes in 
the plane of its cross-section, and torsion about its 
centroidal axes. 

The stiffness matrices for these elements are of 
order 12 x 12. This matrix size was chosen as each node 
has six degrees of freedom (DOF) and there are two nodes 
per element. Also, it is much easier to eliminate 
redundant degrees of freedom in the structure stiffness 
matrix [Ks], as shown in Equation (3.26), than to keep 
track of which node has how many degrees of freedom. 
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If the local axes are chosen to coincide with 
the principal axes of the cross-section, it is possible 
to construct the 12 x 12 stiffness matrix from 2x2 and 
4x4 submatrices. The force displacement relationships 
for the uniform beam element will now be derived 
directly from the differential equations for beam 
displacements. The stiffness coefficients thus obtained 
will be exact within the limits of the assumptions in 
the general theory of beams subjected to loads. 

The differential equations for the axial 
displacements of the beam shown in Figure 4 have the 
boundary conditions that at X = 0 the displacement is SI 
while at X = L, S ? = 0. 



which can be integrated directly. For the above 
boundary conditions 


Q x = (EA/L)S 1 (4.2) 

from the equations of equilibrium in the X direction it 
follows that 


(4.3) 



39 


then 


k l, l = k 7, 7 = EA/L (4.4) 

and employing equilibrium conditions 

k l, 7 = k 7, 1 =_EA / L (4.5) 

The differential equation for the twist on 
the beam shown in Figure 5 is 

Q =-GJ 4^- (4.6) 

where GJ is the torsional stiffness of the beam cross 
section. By integrating Equation (4.6), Equation (4.7) 
is obtained. 


Q 4 x=-GJ£ + C i (4.7) 

With the boundary conditions Q =0 at X = L, C i is 
calculated 


C 1 ■ °4 L 


(4.8) 
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and since £ =S 4 from Equations (4.7) and (4.8) we 
obtain 

Q 4 = (GJ/L)S 4 (4 

For the twisting moments we have 

Q 10 ~0 4 (4 


Hence 


k 4, 4 " k 10, 10 " GJ/L 


and 


k 4, 10 * k 10, 4 =- GJ/L < 4 

The lateral deflection S on the beam shown in 
Figure 6, subjected to shearing forces and associated 
moments, is given by 


9) 


. 10 ) 


. 11 ) 


12 ) 


S = Sb + Ss 


(4.13) 
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where Sb is the lateral deflection due to bending 
strains and Ss is the additional deflections due to 
shearing strains. These deflections can be expressed 
respectively as: 


EI z (dSb/dx) = Q 2 *-Q 6 (4.14) 


and 


(4.15) 


Integration of Equations (4.14) and (4.15) gives 


3 2 

= Q2 X _ Qa X _ Q 2 EI z X 


+ CjX+C 2 (4.16) 


where As represents the effective beam shear area. The 
constants of integration and C 2 can be found from 
the boundary conditions shown in Figure 6. 

At X = 0 and X = L 


d S 
d x 


d Ss 
d x 


(4.17) 
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and at X = L 


S = 0 


Equation (4.16) then becomes 




(4.18) 
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where 


Q 


6 = 


Q^L 

2 


(4.19) 


and 


* = 


i2. E , I .z, 
G A s l3 


(4.20) 


The boundary conditions for the built-in end are taken 
as & - 0, i.e., the slope due to bending deformation is 
equal to zero. The remaining forces acting on the beam 
can be determined from equilibrium. 

Q 8 -Q 2 (4.21) 


and 


°12 - Q 6 + ° 2 L 


(4.22) 


At X = 0 and S = S 7 , we obtain from Equation (4.18) 


t Q?L J 


S n - (1 + <p )- 

2 V 12£I Z 


(4.23) 
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By using Equations (4.19) and (4.21) through (4.23) and 
by the use of the conditions of symmetry, or the 
differential equations for the beam deflections, it can 
be shown that 

k 2' 2 =K 8 ' 8 = ~ k 8 f 2 =(12EI)/(1+ * )]j3 

(4.24) 

and 


= (6EI) / (1+ 


) L 


(4.25) 


A beam subjected to bending moments and 
associated shears is shown in Figure 7. In order to 
calculate the deflections Equation (4.16) is evaluated 
with a new set of boundary conditions. Let S = 0 at X = 
0 and X = L and at X = L 


d S _ d S s 

d x d x 


(4.26) 


then Equation (4.16) becomes 


ei z s =^- 2 ( x3 - l2x ) + 


(4.27) 
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and 


^2 = 



(4.28) 


The remaining forces acting on the beam can be 
determined from Equations (4.21) and (4.22). At x = 0 

dS b _ dS dS s _ 

17" " d7“ " < 4 - 29 ) 


and Sg becomes 


= Q6( 1 + 0 )L 
EI Z ( 4 + $ ) 


(4.30) 


Thus, the stiffness coefficients are obtained from 
Equations (4.21), (4.22), (4.28) and (4.30). 


(4 + 0)EI z 


k 6,6 = k 12 , 12 


( 1+0)L 


(4.31) 
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-6EI z 

k8 ' 6 - (1 +*)l2 <4.32) 


k 1 2,6 


(2- 0)EI 2 
( 1+ <P)L 


(4.33) 


The stiffness coefficients associated with the 
deflections and can be derived following a 
similar approach. 

The element stiffness matrix [K] is assembled 
from the individual stiffness coefficients k and for a 
beam element are shown in Table 1. If the shear 
parameter 0 representing shear deformation is taken as 
zero in Table 1, Table 2 is produced. The member 
stiffness matrix for a rod element is similarly 
constructed by assembling the individual stiffness 
coefficients. Equations (4.1) through (4.5), as shown in 


Table 3 
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TRANSFORMATION MATRICES 

The member stiffness matrix [K] must be converted 
from the member to the structure coordinate system in 
order to add the contribution of each member stiffness 
to the structure stiffness matrix [Ksj as shown in 
Equation (3.25), 

[Ks] = £ tRjT(l) [K](l) [R]W (3.25) 

i = l 

Consider the body shown in Figure 8. Two sets 
of orthogonal axes with origins at 0 are shown. The Xm, 
Ym, and Zm axes will be taken as a set of member 
oriented axes and referred to as the member or local 
coordinate system. The Xs, Ys, and Zs axes are assumed 
to be parallel to a set of structure reference axes 
identified as the global or structure coordinate 
system. Let ^ij represent the direction cosines of the 
Xm, Ym, and Zm axes with respect to the Xs, Ys, and Zs 
axes, i.e., the cosines of the angles between the member 
and the structure axes. Let [Am] be the set of 
coordinates, written in matrix format, for a point 
expressed in the member coordinate system and [As] the 
set of coordinates for the same point expressed in the 
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structure coordinate system. The relationship between 
[Am] and [As] can be expressed as follows: 


and 


[Am] 


Xm 

Ym 

Zm 


[As] 



(4.34) 


(4.35) 


[Am] = [R] [As] 


(4.36) 


where [R] is a transformation matrix of the form 


* 11 

*12 

*13 

j*21 

*22 

*23 

| 

| *31 

*32 

*33 


(4.37) 
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and *ij represents the direction cosines with the first 
subscript referring to the Xm, Ym, and Zm axes ana the 
second subscript referring to the Xs, Ys, and Zs axes. 

The rotation matrix [R] for a space rod is of 
indefinite form because the orientation of its principal 
axes is unspecified. Let the Xm axis of the member 
coordinates be chosen to coincide with the longitudinal 
axis of the member and the member Zm axis be taken along 
the global Zs axis in order to define its position in 
all cases, including that of a vertical member. The 
rotation matrix [R] for the space rod then takes the 
general form shown in Equation (4.38), where Cx, Cy, and 
Cz are the direction cosines of the member axis with 
respect to the global Xs axis. 





c x 

Cy 

Cz 



[R] 

= 

*21 

*22 

*2 3 

(4.38 




L * 31 

*3 2 

*3 3 



The 

direction 

cosines 

for the last two rows of 

[R] in 

Equation (4.38) 

can be 

found 

directly from 

Figure 

9 by 

geometric 

considerations 

. Alternatively, 
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the transformation may be considered to take place in 
two steps; A rotation about the Ys axis through an 
angle /S and a rotation through an angle y about the Zs 
axis. The matrix [R ^ ] is composed of the direction 
cosines of the intermediate axes with respect to the 
structure axes. 


[R 


0 


C x 

v/cf+clF 


0 



0 1 0 
-c z C x 

\ 0 v/cI+CJ 


(4.39) 


and the matrix [R y ] consists of the direction cosines 
of the member axes with respect to the ft axes. 
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[R 


Y 


'ylc\ + Cl Cr o" 

-Cy VcT+c| 0 

0 0 1_ 


(4.40) 


The rotation matrix [RJ can now be obtained as the 
product of the matrices [R y ] and [R p ] as follows: 

[R] = [R y ] [R p ] (4.41) 

and 


[R] * 


" C x 

— CxCy 
\C.T + Cl 
—Cz 

-y'ci + ci 


Cy 

VcTFci 

o 


Cz 

— CyCz 
\'C? Y + Cl 
Cx 

Vc\ +~ci. 


(4.42) 


with Cx f Cy, Cz being the direction cosines of the 
member axis. The transformation matrix [R] shown in 
Equation (4.42) is valid for all positions of the member 
except when it is vertical. When the member is vertical 
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the direction cosines of the member axes with respect to 
the structure can be determined by inspection and [R] 
becomes 

(4.43) 


[R] = 


0 Cy 0 
-Cy 0 0 

0 0 1 


In this matrix Cy = 1 when the member points "up" and 
equals to -1 when the member points "down". 

The transformation matrix for a space beam may, 
under special circumstances, be the same as that for a 
space rod. However, it generally has its principal axes 
Ym and Zm in skew directions and a more complicated 
transformation matrix need be generated. A space beam 
with the rotations y , p , and a is shown in 
Figure 10. The rotation through the angle a is further 
illustrated in Figure 11 and requires the introduction 
of the rotation matrix [R a ] in which the elements are 
the directions cosines between the member axes Xm, Ym, 
and Zm and the axes: 


[R ct 


10 0 
0 costt sino 
0 -sino coso 


(4.44) 
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The required rotation matrix [R] for the space beam 
element can be obtained as the product of [R y ], 
[Rp ]r and [R a ]: 


[R] = [R y ] [R p ] [R a ] (4.45) 


and 



C x C Y C z 


— CxCr cos a — C z sin a 
\£ 'v C'z 

CxCy sin a —Cz cos a 

\CVTC| 


\C 3 y + cos a 


— CyCz cos a-f-r'.v sin o 


— % C 2 X C 2 Z sin a 


\c- A . + c\ 

CyCz sin a-j-Cx cos o 

\cf+c| 


(4.46) 


where Cx, Cy, Cz are the direction cosines of the member 
axes with respect to the global Xs axis and a is the 
angle between the Y y and Ym axis. 
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For the special case of a vertical member there 
is a rotation about the Zm axis through the angle y 
which may be either 90 or 270. The second rotation is 
about the Xm axis and through an angle a . The 
rotation matrix for the vertical beam member can thus be 
obtained by inspection and consists of the direction 
cosines of the member axes Xm, Ym and Zm with respect to 
the structure axes Xs, Ys, and Zs. 


[RJ * 


0 

•Cy cos a 
Cy sino- 


Cy 0 
0 sino 
0 cos® 


(4.47) 


Again, in tnis case, CY = +1 when the member points 
"up", and CY = -1 when the member point "down". 

Equation (4.47) is valid for all possible vertical 
orientations of the member. 

Tne orientation of a particular member may be 
such that the angle a may not be readily available and 
thus a tecnnique for calculating the angle of rotation 
must be developed. A suitable method for specifying 
is to use the coordinates of a point that lies in one of 
the principal planes (Ym-Xm), but not on the axis of the 
member itself. This point and the Xm axis will define 
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without ambiguity a plane in space. All that is 
necessary now is to obtain an expression for the angle 
of rotation a , shown in Figures 10 and 11, in terms of 
the coordinates of the given points and the coordinates 
of the ends of the members. Let Xps, Yps, and Zps be 
the coordinates of point P with respect to the structure 
axis. Xp, Yp, and Zp are the point coordinates of P, as 
shown in Figure 11, and Xj, Yj , and Zj represent tne 
coordinates at the end of the member (node j), then 

Xps = Xp - Xj 

Yps = Yp - Yj (4.48) 

Xps = Zp - Zj 


The coordinates of P with respect to the axes Xs, Ys, Zs 
(principal axes in the structure or global system) can 
be obtained by 


Xpy 



YPY 

I ZpY 

= [R y 



Xps | 

Yps 

Zps 


(4.49) 
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Then, from the geometry of the beam shown in Figure 11, 
the following expressions for the sine and cosine of the 
rotation angle a in terms of the coordinates of the 
given point and those of the ends of the member are: 

Z PY 

Sina = - (4.50) 

\/ r P Y +z PY 

and 

y PY 

cos a = - (4.51) 

7 y p/ z py 

The preceding discussion dealt with non-vertical 
members. In the case of the vertical member, it is 
possible to calculate cos a and sino directly from 
the coordinates of the point P 


Zp s 

v/xfl+z’ps 


s in a 


(4.52) 
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and 


-A ps 

coso = ' , ~ Cy (4.53) 

>«P8 +Z P S 

Equations(4.47) , (4.52), and (4.53) are general formulas 
valid regardless of whether the member is pointing "up 
or "down". 

It can be summarized that the rotation [R a ] 
for a space beam member can be generated using either 
the rotation angle a directly, or from the coordinates 
of a point P that helps to identify a principal plane 
(in ASTRA the Xm-Ym plane). If the angle a is zero, 

[R] reduces to one of the expressions generated for a 
space rod member. Equations (4.42) and (4.43). If a is 
not zero then [R] takes the form of Equation (4.46) 
or (4.47). 



Chapter 5 


A PROGRAMMERS GUIDE TO "ASTRA" 

With the exception of very elementary problems, 
the application of the finite element method necessarily 
implies a considerable amount of computation, thus, the 
use of a digital computer is implicit. In this chapter 
the computer program ASTRA will be described. Its 
sequence of computer operations, handling and storage of 
matrices, and the special purpose subroutines are 
listed. A general program flowchart can be found in 
Appendix D. 


The Program 

In order to circumvent the memory size 
limitation of the APPLE II microcomputer (48k bytes), 
the program ASTRA was divided into a number of 
subprograms and subroutines. The difference between 
subprograms and subroutines as defined is that a 
subroutine is a program called for from a subprogram, 
while a subprogram is a complete computer program that 
is a part of the main program. ASTRA has been divided 
into 4 separate subprograms called ASTRA-1 through -4. 

Each of these subprograms must be executed in sequence. 

58 
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In the computer system only one floppy disk was 
available. This forces the user to separately load and 
execute each of the four ASTRA subprograms. Should 
another floppy disk drive or hard disk drive become 
available to store both the programs and the data, the 
programs could all be set to load and run automatically 
by the use of a text editor file. This file is the IBM 
equivalent of a CLIST type file and allows for 
instructions to the computer to be stored and executed 
as though it were a regular computer program. Should 
Pascal, a higher level computer language, be used 
instead of the Applesoft Basic interpreted language used 
on the APPLE II, considerable amount of core storage 
could be saved and program execution would be speeded 
considerably. 

The program ASTRA assigns six degrees of 
freedoms (DOF) to each nodal point. These degrees of 
freedom, 1 through 6, are equivalent to three 
translations and three rotations in and about the global 
X, Y, and Z axis respectively. In turn, the assigned 
degrees of freedom are not activated until the 
individual nodes are called for by an individual 
element. This feature allows for faster program 
execution by reducing the equilibrium equations that 
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must be solved for. In addition, it allows the use of 
dummy nodes, i.e., nodes which have been given 
coordinates in the program but are actually not used by 
any of the elements. This provides for the addition and 
deletion of elements and nodes without the necessity for 
renumbering all the nodes in the structure. 

For example, consider the structure shown in 
Figure 13 and its corresponding model. The individual 
numbers shown are the node numbers and those enclosed in 
parenthesis are the element numbers. Let elements 1 
and 2 be rod elements, and element 3 be a beam element. 
ASTRA would assign 18 (6 DOF x 3 nodes) degrees of 
freedom to the structure. Of these, only 15 degrees of 
freedom would be activated by the elements. Six DOF 
each at nodes 2 and 3 by the beam element and three DOF 
at node 1 by the rod elements would be activated. 

In order to check whether a DOF is active or 
not, the user is referred to the ASTRA print out when 
DOFIU (degrees of freedom in use) in the computer 
program is displayed. 

In ASTRA the high speed available memory of the 
APPLE II (RAM) is used primarily to store the arrays 
currently in use, and to provide a buffer for the 
solution of equations. The low speed memory provided by 
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the floppy disk drive is used for the storage of element 
properties, structure stiffness and transformation 
matrices, input data, and any large array. The total 
capacity of the program is thus limited by the amount of 
RAM and disk storage available while its speed of 
execution is governed by the total number of DOF in the 
structure. Because the speed of execution of ASTRA is 
proportional to the number of equilibrium equations to 
be solved, i.e., the number of active DOF in the 
structure, to increase the speed of execution it is 
recommended that all DOF not compatible with the 
elements connected to the nodal points be suppressed. 
This will reduce the number of degrees of freedom at 
each node and thus speed the execution of the program by 
decreasing the order of the structure stiffness matrix 
[Ks], and allowing for the analysis of larger structures 
in less amount of time. For example, in a plane 
structure in the X-Y plane, suppress all translations in 
the Z direction and rotations about the X axis. This 
will reduce the number of degrees of freedom at each 
node by a factor of two and the order of the matrix [Ks] 
by a factor of 2n, where n is the number of nodes in the 


structure. 
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Often interesting problems involving large 
matrices cannot be solved because they are too large and 
are either impossible or very expensive to invert with 
available computer storage. If the matrices involved 
are sparse matrices, there are ways of optimizing their 
storage so as to allow for the solution of larger 
problems. Sparse matrices are matrices having only a 
small percentage of nonzero elements. If these elements 
are grouped in a narrow band on both sides of the matrix 
diagonal, then they can be processed in packed form with 
only the bandwidth, or 1/2 the bandwidth stored with the 
necessary indexing information. Packed storage confers 
the following advantages: 

1) Larger matrices can be stored and handled in 
internal storage. 

2) Access time is faster with RAM; therefore, 
the packed form is preferred since it relieves the need 
for external storage. 

3) A substantial amount of time is saved if the 
operation involving elements outside the bandwidth (zero 
value elements) are not performed. 
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4) The inverse of a given matrix expressed as a 
product of elementary matrices stored in packed form 
usually needs less storage than the explicit inverse of 
a matrix. 

The ideal storage would be one that minimizes 
both the total storage used and the total computation 
time. In general, these two requirements are 
incompatible and a trade-off must be made. 

In ASTRA the technique of packed storage was 
used to store the structure stiffness matrix [Ks] at the 
expense of increased running time. Instead of storing 
all the terms in the lower triangular matrix, only those 
terms within the lower half bandwidth are stored. This 
was necessitated by the limited amount of storage space 
(RAM) within the computer and the single floppy disk 
drive available. It was also done in order to reduce as 
much as possible the amount of disk storage. According 
to the Applesoft programming manual [6], it takes 10 
bytes of computer memory to store a real number. With 
the floppy disk drive, a minimum of 16 bytes of storage 
were required to store a real number while preventing 
records from overwrite. In addition, to save space, 
some matrices made up of submatrices were selected for 
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storage while others were kept in RAM. An example of 
this is the displacement transformation matrix [R] 
described by Equation (3.4). 


[ r J - i 

N N T N 
[_N N N T_ 

In this matrix, [T] is a transformation matrix 
and [N] a null matrix. Only the non-zero terms of the 
primary submatrix [TJ are stored on disk storage. When 
the matrix [R] is required by the program, [T] is 
recalled from storage and [R] is generated. 

The features outlined above greatly extend the 
time required for running ASTRA, however, they allow the 
implementation of the program on a microcomputer. 

Unless the computer programs are written in machine 
language and optimized for the system and/or a second 
floppy disk drive, a hard disk drive or bubble memory 
becomes available as a source of storage, this is the 
only way in which even simple problems can be solved on 
the APPLE II without extensive software development and 



addition of RAM. 
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Two different types of files can be created on 
the APPLE II disk system; a program file and a text 
file. The program file is used to store the computer 
programs (example: ASTRA-1). The text file is more 
flexible and can be used to store data in sequential or 
random modes, and computer programs to be run from the 
disk. 

The random storage data file are used 
extensively in ASTRA to handle the large matrices 
generated. A special feature of some of the storage 
subroutines in ASTRA is the storage of only non-zero 
real numbers. Through experimentation, it was found 
that it requires a minimum of 16 bytes on disk to store 
a real number using random access files without the risk 
of records overwriting one another. This applies 
whether the number is a large quantity or a zero. Since 
the floppy disk presently available can handle only 116k 
bytes of storage, it was imperative that the amount of 
data to be stored be reduced to the absolute minimum. 

A feature of most matrices handled by ASTRA 
during the execution of the program is that they are 
large, sparse, symmetric matrices. By taking advantage 
of these characteristics, it is possible to write 
computer subroutines that store and retrieve the 
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non-zero real numbers from the lower triangular matrix 
only. In order to achieve this feature, advantage is 
taken of two features of the APPLE disk system, namely, 
data is stored in a random storage data file and the "on 
error goto" command. 

The number access files on the APPLE disk are 
like a honeycomb, i.e., a collection of equal size 
cells. Each time that a record is stored in these 
cells, the system energizes a 250 byte block and writes 
the required information in it until that particular 
block is filled. This enables it to write information 
in random files, address them, and retrieve them. If in 
a five by five matrix the information from row three, 
column one is to be stored, the corresponding disk 
address number will be two (rows) times five 
(records/row) plus one (column) or disk address number 
eleven. For retrieval the disk will search directly the 
given address. When the computer tries to retrieve a 
record from an address that has not been used, a fatal 
error will be generated. When this occurs, the use of 
the "on error goto" command, similar to the fortran EXIT 
statement, will enable the program to bypass the error 
and continue with the program execution. The "on error 
goto" is used to avoid having an error message printed 
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and program execution halted. The command sets a flag 
in the computer which causes an unconditional jump to a 
program line indicated by a pre-established line 
number. Although the "on error goto" command is very 
helpful, extreme care must be taken when it is being 
used as it has side effects on the working of the 
computer. The use of the command causes the do loops 
and subroutine pointers to be disturbed. It can also 
cause problems with the buffers which handle the 
internal mathematical calculations and affect the 
available memory and working of the program. For these 
reasons it is recommended that it is used only when 
absolutely necessary. 

The element stiffness and displacement 
transformation matrices are used in two different 
sequences of operation in ASTRA. They are first used in 
the calculation and assemblance of the structure of a 
stiffness matrix 

[Ks] = £ [R] T(l) [K] (l) [R] (l) (3.25) 

1=] 

and in the calculation of the element forces 
[Q] W = [K] (l) [R] (i) [S J 


(3.33) 
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In order to speed the execution of the program, 

(i) (i) 

it is necessary to store both matrices [K] and [R] 
for each element. However, the storage of these 
matrices would take up a considerable amount of the 
available disk storage space. Advantage was taken of 
the special characteristics of each matrix for more 
efficient operation. For example, the element 
transformation matrix is made up of a series of real and 
null matrices 


and 


M 


1 T N N N j 
: N N N N j 

In n t n 

: N N N N i 


(rod elements) 


m 


T N N N 

N T N N' (beam elements) 
IN N T N ! 

|_N N N Tj 


The submatrix [T] is a three by three real matrix while 
[N] is a null matrix. In order to save on storage 


space, only the non-zero elements of the transformation 
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matrix [T] are stored on the disk. Separate subroutines 
are then used to generate the displacement transformation 
matrix R according to the type of element being 
used. 

A similar technique is used in the storage of 
the element stiffness matrix [K], For a rod element, 
see Table III, only a few non-zero terms need to be 
stored. The complete matrix can be generated from a 
single term AE/L which is stored on the disk memory. 

For beam elements, the stiffness matrix is of a 
more complex format andmust be stored differently. A 
subroutine was written that stores the non-zero terms of 
the lower triangular matrix. This procedure takes a 
somewhat longer retrieval time than that used for 
storage of the full matrix as the full matrix must be 
generated from the values stored. However, this 
procedure has the advantage that a substantial amount of 
disk memory is saved when compared to storing the 
complete twelve by twelve matrix. 


MULTIPLICATION SUBROUTINES 

Multiplication of two matrices [C] = [A] [B], 

where [A] is an "1 x m" matrix and [B] is an "m x n" 
matrix, is defined by 
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m 

c ij ■ Z A ik B kj 

k = l 

The number of mathematical operations required to 
complete the multiplication is of the order of 
"2 x 1 x m x n". Since many of the matrices with which 
ASTRA deals with are made up of a combination of real 
and null submatrices it is convenient to write specific 
matrix multiplication subroutines for them. For 
example, in the calculation of the structure stiffness 
and load matrices for the different elements, 

[K g ] = 2 [R]T(i) [K](1) tR](l) (3.2! 


[Q] (l) = [K] (i) [R] (l) [S a 3 (3.33) 

special subroutines were written instead of using a 
general routine. The added software required is more 
than offset by substantial time savings in the 
multiplication operations when compared to a full 
multiplication using all terms. 
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Solutions of Simultaneous Equations 

The method used in ASTRA for the solution of 
simultaneous equations is the well known Cholesky's 
method of triangular decomposition. This method was 
chosen because row or column interchanges are not 
required, keeping round-off errors small. In addition, 
this method considerably reduces the number of 
mathematical operations required to accomplish the 
matrix inversion. For most applications of finite 
element theory, the stiffness matrix [K] has the 
important characteristics of being positive definite and 
symmetric. With these properties, Cholesky’s method is 
optimal as it requires only r2 + n 

operations,where n is the order of the matrix, which is 
approximately one fourth the number required by 
Gaussian Elimination. 

Node Point Numbering to Exploit Sparsity 

For a given structure all numbering schemes lead 
to the same size of stiffness matrix [K], however, 
different numbering schemes lead to different 
arrangements of non-zero terms. In order to obtain the 
smallest possible bandwidth it is recommended that the 
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structure be numbered across the "short axis," i.e., the 
side with the least amount of nodes, as this will 
produce the narrowest banded diagonal possible. For an 
example, see Figure 12. In ASTRA the members are also 
automatically redefined by the CPU with node I set equal 
to the number of the lower node and node J set equal to 
the higher node value. This is done in order to 
facilitate the calculation and to avoid adding a 
stiffness element above the diagonal. The node number 
is not changed, only the order of input of the nodes. 
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RECOMMENDED IMPROVEMENTS TO ASTRA 

Recommended improvements to this first version 
of the computer program ASTRA can be categorized into 
four general areas: 

1. Equipment 

2. Programming language 

3. Computer program improvements 

4. Pre and post processor subroutines 

Equipment 

The minimum equipment that is required to run 
ASTRA in any other than a laboratory environment is a 
48K APPLE II or similar microcomputer, two floppy disk 
drives, a printer, and a TV monitor. Two floppy disk 
drives are a minimum because one should be used to store 
the programs and input information, and the other should 
be used to store the data generated. 

At the time of this writing, the APPLE Computer 
Corporation has announced the introduction of the 
APPLE III series of microcomputers. This would be an 
ideal computer around which to build a new system as the 

amount of memory would be increased from 48K bytes in 
73 
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the APPLE II to 120K bytes. The APPLE III has an 
integral disk operating system (DOS 3.3) with a single 
floppy disk that has a 20 percent (16 versus 13 sectors/ 
track) increased capacity over the old disk drive 
(DOS 3.2). In addition, it has built-in series and 
parrallel interface boards and allows printing of 80 
characters on the monitor versus 40 on the APPLE II. 

Upgrading the data storage/retrieval system is 
also an area of consideration for improvements to 
ASTRA. The addition of more low speed memory in the 
form of disk drives would greatly increase the capacity 
of the program to handle larger structures and speed the 
program execution by allowing the storage of more 
complete matrices and other generated data. This can be 
achieved in any of three ways: addition of more floppy 
disk drives (up to a maximum of 8), replacement of the 
floppy disk drives by a hard disk, or a bubble memory 
unit. The addition of more floppy disk drives is not 
recommended as it would cost a substantial amount of 
money, increase the complexity of the system, require 
software changes, and would not result in a worthwhile 
improvement in the size of memory as could be obtained 
by other means. Instead, the preferred alternative is 
the replacement of the floppy disk drive by a hard disk 
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or bubble memory. Although the initial cost may be 
greater, it is believed that in the long run it would be 
more economical and trouble free. Memory storage would 
then be increased to 10M bytes allowing the storage of 
the complete element stiffness and transformation 
matrices and other generated data. This would increase 
the speed of execution of the program and improve the 
quality of the program printout by allowing more 
information regarding the analysis completed to be shown. 

The printer should be a commercially available 
matrix dot printer capable of both printing and drawing 
graphics so as to increase its versatility and drive 
down the costs and complexity associated with a separate 
printer and plotter unit. Graphics-capable printers are 
not as accurate as table plotters but are quite adequate 
for most applications. Table plotters are more accurate 
and neat and have the added option that graphics can be 
plotted in color. However, they are expensive to 
purchase and maintain and usually require the 
development of specialized software for their 
operation. A thermal type printer is not recommended 
unless it has graphics capability and is used as a 
plotter only since the quality of print is not durable 
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and the cost of the paper makes it prohibitive to print 
large quantities of computer output. 

The TV monitor should be a color monitor as it 
can be used to advantage when implementing computer 
graphics. Also, a worthwhile addition to any computer 
system would be a digitizer. The digitizer would allow 
for geometric data to be entered directly from the board 
by means of an analog pencil and would eliminate the 
need for typing in the data. 

Programming Language 

The change in programming language from 
"Applesoft" Basic to UCSD Pascal is highly desirable. 
"Applesoft" Basic is a good lower level programming 
language. However, programs written in "Applesoft" are, 
when compared to other higher level programming 
languages such as Pascal or FORTRAN, bulky to store and 
highly inefficient in computing, output handling, and 
editing operations. 

The substitution of "Applesoft" with Pascal 
would accomplish the following changes: 

1. Increase amount of available RAM by 16K 
bytes through the use of a language card. 
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2. Provide for better utilization of RAM. 
Programs written in Pascal are said to use 50 percent 
less RAM to store the program than those written in 
"Applesoft." 

3. Reduction in running time. Pascal is 
reputed to run up to twenty times as fast as "Applesoft." 

4. Increased program editing and output format 
capacity. 

5. Compatibility of the program with other 
systems having Pascal as a user's language. 

Program Improvements 

This first version of ASTRA is simple in scope 
and needs improvements in order to enhance its ability 
to handle more complex problems faster. In order to 
improve the program, it is first suggested that the 
computer system used be upgraded to include a minimum of 
two floppy disk drives or a hard disk drive, and UCSD 
Pascal programming language. If these changes are not 
implemented, it is believed that the improvements to be 
accomplished will not be worth the time spent in 
carrying them out. 
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These are five areas in which ASTRA can be 
improved: 

1. The number and types of structural elements 
available 

2. The capability to conduct thermal and 
dynamic analysis 

3. The capacity of the program to handle more 
degrees of freedom 

4. The program execution speed 

5. Pre and post-processing 

The addition of more types of structural 
elements to include membrane, plate, and three 
dimensional isoparametric elements would greatly 
increase the versatility of the program. In this 
present version, ASTRA can only handle rod and beam 
elements. Many structures must thus be idealized to a 
large extent in order to allow their analysis using 
ASTRA. The addition of the above mentioned elements 
would allow for the analysis of more varied types of 
structures while allowing for more accurate analysis to 
be accomplished. 

A worthwhile research topic would be to explore 
the additions of options to perform thermal, dynamic and 
non-linear analysis. Adding the capability to ASTRA of 
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performing thermostructural analysis should not present 
much trouble as the equations are simple in nature. The 
addition of dynamic and non-linear analysis would be a 
more difficult project as it is unknown at this time 
whether the accuracy of the APPLE II would be sufficient 
to solve the necessary equations and extract the 
appropriate eigenvalues in other than very elementary 
structures. 

In order to increase the capacity of the program 
to handle more degrees of freedom, the first steps 
mandated would be to convert the USCD Pascal and employ 
a hard disk. It is not known at this time what the 
absolute capacity of the APPLE II would be in terms of 
degrees of freedom when using Pascal. How much 
remaining RAM would be available after the program is 
compiled for the storage of matrices is the key 
question. When using a hard disk, the amount of time 
necessary to solve the problem rather than the amount of 
available memory may well become the deciding factor in 
limiting the program capacity. 

Data retrieval and disk space requirements in 
data storage/retrieval operations can also be reduced by 
storing the arrays using a combination of sequential and 
random access files. Data arrays which do not require 
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access at a particular location can be stored in disk 
using sequential file storage. This eliminates the 
problem of records overwriting reduces the number of 
bytes required for storage by not having to use a 
specified record length and simplifies the software by 
eliminating the need for calculating the addresses of 
records stored in the disk. 

ASTRA is sadly deficient in the availability of 
pre and postprocessor subroutines. These are 
subroutines that would help the user prepare and check 
the input data and then interpret the program results. 
These tasks are usually accomplished in the larger 
finite element programs by element and nodal point 
coordinate generating subroutines, and plotting programs 
that draw the structure and its deformed shape, and 
temperature and stress profiles. When programming 
plotting subroutines, the user is warned that due to the 
allocation of memory by the APPLE II to graphics, 
overwrites between the graphics data and information 
stored in arrays may occur, and some of the information 
stored in the arrays may be lost. Thus, care must be 
exercised when using the APPLE II high-resolution 
graphics capabilities. 
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USER’S GUIDE TO ASTRA 

The Subroutine IASTRA 

The subroutine IASTRA (Input to ASTRA) was 
written for the specific purpose of helping the user 
prepare the input data to ASTRA. This section of the 
thesis is intended to serve as a user's guide to both 
the preprocessor IASTRA and the program ASTRA. As these 
are general computer codes intended for use with a 
variety of computer systems, no attempt will be made to 
tailor the user's instructions to a particular system. 
Instead, a comprehensive review of the preprocessor, 
IASTRA presented together with a general outline for 
running ASTRA will be described. 

In order to prepare the input data to ASTRA, all 
that is required of the user is that he answers the 
questions posed by the preprocessor IASTRA. The first 
step in running ASTRA then is to load and run IASTRA 
(enter the input data). IASTRA will then require the 
user to select an option by entering the needed option 
number (and then hitting the return key), there are five 
user's options in the preprocessor: 
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1. Create a new data file 

2. List input data 

3. Modify an existing data file 

4. Copy an existing data file 

5. Exit 

After the option is completed, the program will 
return to wait for the next option number. Option five, 
"EXIT," releases the program. At that time ASTRAl, 
ASTRA2, ASTRA3, and ASTRA4 should be loaded and run. 

Option No. 1: Create a New Data File 

As the name implies, this option is used to form 
a new set of input data. The following is a list of 
terms asked by the program when using the option. 

1. DWG ID CODE? 

This is a code number or name under which 
the input data will be stored on the disk. The code can 
consist of a combination of alphanumeric characters less 
than 10 in number. 

2. DATE? 

Self explanatory. Date < 10 characters 

3. YOUR NAME? 

Self explanatory. Name < 10 characters 
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4. NUMBER OF ELEMENT CARDS? 

Total number of element generating cards 

5. DATA FOR ELEMENT CARDS? 

The element generator cards are used to 
enter information about the structure elements. One 
card can be used to generate more than one element. The 
information on the element cards includes eight fields 
of data as shown below andin Table (IV). 

Card #ABCDEFG 


where 


A = element type (1 = rod, 2 = beam) 

B = group number 

C = number of elements to be generated by 
the card 
D = node i 

E = increment for node i to be used in the 
generation of other elements 
F = node j 

G = increment for node j to be used in the 
generation of other elements 
Example: Generate the element cards for 

the structure shown in Figure (16b). 
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Let the material and geometry of elements 2 
through 5 be the same as element 1, and those of 
elements 7 through 13 be the same as element 6 
respectively. This structure can be divided into two 
material groups. Group one with elements 1 through 5, 
and group 2 with elements 6 through 13 respectively. 
The element cards will look as shown in Table 4 in the 
Appendix. 

6. ENTER MAXIMUM NUMBER OF NODES? 

Maximum number of nodes in the mesh 

7. NODAL POINTS COORDINATES 

This card consists of five fields as shown 

below 

Card # A B C D 

where 


A = nodal point number 
B = X coordinate 
C = Y coordinate 
D = Z coordinate 

All coordinates are in the global or 
structure coordinate system. 
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8. NUMBER OF SUPPRESSION CARDS? 

Each node number that requires to have some 
or all of its degrees of freedom restricted must have 
one suppression card. The number required at this step 
is the total number of suppression cards. 

9. INPUT SUPPRESSIONS? 

The computer will ask for the node number 
and what degrees of freedom are to be suppressed. Enter 
the node number and answer "Y" or "N" to the questions. 
Default is "N" (free to deflect). 

10. INPUT MATERIALS AND SECTION PROPERTIES? 

This part requires that the user enters the 
total number of material groups and then the specific 
material properties required for each group. In each 
structure there can be a number of structural elements 
possessing the same material and geometric 
characteristics. In order to avoid entering one 
material and geometric data card for each element, those 
having the same properties are given the same group 
number and corresponding mechanical and geometric 
properties entered for the group. 
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11. TOTAL LOAD CARDS? 

For every node that has one or more forces 
or moments acting on it, a load card must be prepared. 
Loads are references on the global or structure 
coordinate system. A right hand system of coordinates 
is used with the applied load being positive when acting 
in the positive direction of the corresponding 
coordinate system. 1 

At this point the xnput option is complete and 
the program saves the information in the proper format 
in the disk. For a description of the procedure used to 
correct errors see option NO. 3, "MODIFY OPTION." 

Option No. 2: List Option 

This option will print the input information of 
any given file on the system printer. After the option 
has been completed, the printer must be reset manually. 
In order to run this option, the printer is first 
energized and then the option is executed. 

Option No. 3; Modify-Edit Option 

As the name implies, this option is used to 
change an existing data file. Its uses include . 
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correction of errors, and/or additions and deletions of 
nodes, element material properties, nodal suppressions, 
and nodal loads. This option is easy to use and all 
that is usually required is reentering the complete 
data card to be changed or added to the data file. The 
procedure for deleting data, for example an element 
card, is somewhat tricky and will be discussed below. 

In order to delete a data card, the program will discard 
the last data card for that particular section. For 
instance, let us imagine a set of element cards that 
contain five cards. If the delete option is exercised, 
the program will discard element No. 5. However, if the 
card that we want to discard is card No. 2, the 
procedure is somewhat different and must be executed in 
two steps as follow: 

1. Change element card No. 2 by entering the 
last card in the section as card No. 2. 

2. Exercise the delete option to discard the 
last card in the section. 

The procedure described was necessitated by the 
limited amount of time available to develop the 
appropriate software. It is recommended that this 
procedure be simplified in the future in order to avoid 


confusion. 
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Option No. 4; Data File Duplication 

This option is used to duplicate an existing 
data file onto another disk. The file will be saved 
with the existing file name. The computer will ask the 
user to place in the disk drive the disk containing 
the file to be duplicated, and the file name. Once this 
is done it will ask that a new disk be inserted in the 
disk drive and, upon hitting return, it will duplicate 
the file in the new floppy disk that was just inserted. 
Once the duplicating procedure is completed the program 
can be restarted by hitting return. It is recommended 
that all data files in use be duplicated in order to 
avoid loss of data if the disk containing the original 
file is lost or its information destroyed. 

Option No. 5: Exit Option 

On the completion of any of the previously 
discussed options, the program executes a loop and 
restarts itself. The EXIT option causes the program to 
exit the loop and thus terminates the execution of 


IASTRA. 
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Stiffness Matrix K for a Beam Element 
with Shear Deformations 
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without Shear Deformations 
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Table 3 

Stiffness Matrix for a Space Rod Element 


10 0-100 


0 0 0 0 0 0 



0 0 0 

-1 0 0 


0 0 
1 0 


0 

0 


0 0 0 0 0 0 


0 0 0 0 0 0 


Table 4 

Sample Element Generating Cards 


CARD 

ELEMENT 

GROUP 

# OF 

NODE 

INCREMENT 

NODE 

INCREMENT 

# 

TYPE 

# 

ELMTS 

i 

IN i 

J 

IN J 

1 

2 

1 

2 

4 

1 

5 

1 

2 

2 

1 

2 

7 

1 

8 

1 

3 

2 

1 

1 

10 

0 

11 

0 

4 

2 

2 

3 

1 

3 

4 

3 

5 

2 

2 

3 

2 

3 

5 

3 

6 

2 

2 

3 

3 

3 

6 

3 


COMPANY MODEL RAM ROM TERMINAL STORAGE SOFTWARE COMMENTS 
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RADIO TRS-80 32K 12K VIDEO CASSETTE MONITOR Z80 BOARDr 
SHACK 64X1A CHAR 5 INCH BASIC DOCUMENTATION AVAILABLE 

KEYBOARD FLOPPY FORTRAN WIDE RANGE OF SOFTWARE 
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3LIST 

10 REM 6 DEC 1980 
20 REM PROGRAM IASTRA 
30 HOME t CLEAR 

40 GOSUB 3001 REM ALLOCATE ARRAYS 
50 HTAB (10)*. PRINT 'PROGRAM IASTRA" 

60 PRINT 

70 HTAB (5): PRINT "ANALISING STRUCTURES WITH APPLE" 

80 PRINT 

90 VTAB (5)! PRINT ‘THIS PROGRAM IS USED TO PREPARE THE INPUT DATA TO THE 
FINITE ELEMENTS ANALYSIS COMPUTER PROGRAM ASTRA" 

100 VTA£„<9>: PRINT "CHOOSE ONE OF THE AVAILABLE OPTIONS:" 

110 VTAB (11): PRINT TAB- 5>*"1. CREATE A NEW INPUT DATA FILE.' 

120 PRINT TAB< 5 if"2. LIST INPUT DATA. 

130 PRINT TAB( 5)» r 3. MODIFY AN EXISTING DATA FILE." 

HO PRINT TAB( 5>»'4. COPY OPTION." 

160 PRINT TAB( SWS. EXIT FROM THE PROGRAM.' 

170 VTAB (19): INPUT "ENTER OPTION NUMBER ?"PAS 

180 IF (A* < 1) OF: (AX > 6) THEN HOME I PRINT 'WRONG OPTION NUMBER* TRY 
AGAIN": GOTO 100 

190 HOME J IF AX = 1 THEN GOSUB 380: REM INPUT OPTION 

200 IF AX = 2 THEN GOSUB 3370*. REM LIST OPTION 

210 IF AX = 3 THEN GOSUB 5030: REM MODIFY OPTION 

220 IF AX = 4 THEN GOSUB 4920: REM COPY OPTION 

240 IF AX = 5 GOTO 2701 REM EXIT 

250 HOME : CLEAR 

260 GOTO 10 

270 HOME PRINT "END OF JOB' 

280 END 
290 REM 

* * * i** »* * «* **«* »***« 

»»**»*******»***»» 

30C REM ALLOCATION OF ARRAYS 
310 DIM AAID$( 5 >f ABEGEh’XC lQOf 9 ) 

320 DIM ACNPCO<200*3) 

330 DIM AGNDOF/.( 50*6 )i REM NODAL SUPPRESSIONS 

340 DIM AJPNODE( 200*6)! REM LOAD AT NODES 

350 DIM AKMSPROP( 10*9)! REM MAT t> SECTION PROPERTIES. 

360 RETURN 

370 rem mmmmimmmmmmm************** 

380 REM SUBROUTINE INPUT 
390 HOME : VTAB (3) 

400 GOSUB 3000: REM INPUT HEADING 

420 HOME : PRINT "ENTER THE FOLLOWING INFORMATION FOP. EACH ELEMENT GENERA 
TOR CARD IN FREE FORMAT." 

430 VTAB (7): PRINT "NOTE: ANY ERRORS NOT CORRECTED BEFORE THE LINES AR 

E ENTERED CAN BE TAKEN CARE OF AFTER ALL CARDS HAVE BEEN ENTERED.': PRINT 

440 PRINT "HIT RETURN TO CONTINUE" 

445 INPUT A« 

460 HOME 

470 GOSUB 3380: REM ELEMENT CARDS 

480 GOSUB 3680*. REM NPCOS 

490 GOSUB 560*. REM INPUT SUPPRESSIONS 

500 GOSUB 1010: REM INPUT MAT. i SECT. PROP. 

510 GOSUB 1390: REM INPUT LOADS 
520 HOME 

530 PRINT "INPUT SUBROUTINE COMPLETED" 

540 GOSUB 1750*. REM SAVE DATA 






550 RETURN 

560 rem mmmmmmmmmmmm****** 

570 REM INPUT SUPPRESSIONS 
580 HOME 

590 INPUT "NUMBER OF SUPPRESSIONS CARDS ? r ?AH% 

600 IF AH% > 100 THEN PRINT "TOO MANY NODES SUPPRESSED <tK*100)f TRY AGA 
IN": PRINT : GOTO 590 
610 HOME 

620 VTAB (5)J PRINT "ENTER OPTION NUMBER" 

630 VTAB (7): PRINT TAB( 3)? r 'l* ONLY ROD ELEMENTS ARE USED♦ 1 

640 PRINT TAB< 3)?'2. A MIXTURE OF BEAM AND ROD ELEMENTS ARE USED/ 

650 INPUT "ENTER OPTION NUMBER ? ,: >I1% 

660 IF Oil < 1) OR. (11% > 2) THEN HOME l PRINT "WRONG OPTION NUMBER? TR 
Y AGAIN ■: VTAB (3): GOTO 620 
670 121 = 1 
680 GOSUB 790 
690 HOME 

700 INPUT "DO YOU WANT TO ADD ANY MORE SUPPRESSIONS ?( Y OR N ) ’ >A$ 

710 IF ( A$ = "N r ) GOTO 760 

720 IF (A* < > "Y r ) THEN PRINT “URONG INPUT "\ GOTO 700 

730 INPUT "HOi* MANY EXTRA NODES?' ?I3% 

740 in = mi + i:ahx = mi + 135; 

750 GOTO 680 
760 RETURN 
770 REM 

780 REM INPUT SUPPRESSIONS-! 

790 FOR I = 111 TO AHZ 

800 HOME : PRINT "ENTER NODE NUMBER AND Y GR N TO THE QUESTIONS/ : VTAB ( 
3) 

810 PRINT 

820 PRINT "NODAl SUPPRESSIONS CARD * *♦. ‘ »I 
830 VTAB (6) 

840 INPUT "NODE NUMBER ? "?AG%(I?G) 

850 INPUT ’SUPPRESS X ? r »I$ 

860 IF (II = "Y r ) THEN AG%(I*1> = 1 
S7C INPUT "SUPPRESS Y ? 1 fIS 
880 IF <1$ * ”Y r ) THEN AG%(I»2) = 1 
890 INPUT "SUPPRESS Z ? B »I$ 

90C IF C IS = r Y 1 ’ ) THEN AG%(I?3) = 1 

910 IF 11% * 1 GOTO 980 

920 INPUT ’SUPPRESS THETA X ? r ?I* 

93C IF (1$ = r Y ! ) THEN AG%(I?4> = 1 
940 INPUT "SUPPRESS THETA Y ? "?I$ 

950 IF (1$ = "Y”) THEN AGZ(I»5) = 1 
960 INPUT "SUPPRESS THETA Z ? r ?I$ 

97C IF (1$ = "Y’ ) THEN AGttlfo) = 1 
980 NEXT 
990 RETURN 
100C REM 

1010 HOME : PRINT "INPUT MAT * i SECT* PROPERTIES" 

1020 VTAB (3 ) 

1030 INPUT "NUMBER OF GROUFS ? 1 imi 
1040 HOME 

1050 VTAB (3)J PRINT "ENTER OPTION NUMBER’ 

1060 VTAB ( 7)1 PRINT TAB( 3>? 2 1* ONLY ROD ELEMENTS ARE USED* 1 

1070 PRINT TAB( 3)J P 2* A MIXTURE OF RODS AND BEAM ELEMENTS ARE USED* r 

1080 VTAB (11) 

1090 INPUT "ENTER OPTION NUMBER ? r ?Il% 

HOC IF (11% < 1) OR (11% > 2) THEN HOME : PRINT "WRONG OPTION NUMBER? T 
RY AGAIN r : VTAB (3); GOTO 1050 
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1110 121 = 1 

1120 FOR I = 121 TO AF1 
1130 HOME 

1140 PRINT "NAT. 6 SECT. PROPERTIES FOR CROUP ♦ ' *1 
1150 IF Ill = 1 THEN AK(I»0) = 15 GOTO 1170 
1160 INPUT "ELEMENT TYPE ? SAK(I.O) 

1170 VTAB (3)1 INPUT 'CROSS-SECTIONAL AREA ? 'fAK(Ifl) 

1160 IF Ill = 1 GOTO 1220 

1190 PRINT "HOHENTS OF INERTIA ABOUT LOCAL AXIS"! PRINT 
1200 INPUT "I <Y'-Y' ) ? '!AK(Ir2X 
1210 INPUT *1 (l'-T ) ? ,: IAN(If3) 

1220 PRINT ! INPUT "MODULUS GF ELASTICITY E T "!AK(I»4> 

1230 IF Ill = 1 GOTO 1290 

1240 INPUT "POISSON S RATIO V ? '!AK(I,5> 

1250 INPUT "POLAR MOMENT OF INERTIA J r '»AK(I,6> 

1260 INPUT "X' ? "!AKU»7> 

1270 INPUT * Y' ? SAKUtS) 

1280 INPUT "Z' I "JAN! 1*9} 

1290 NEXT 
1300 HOME 

1310 INPUT "DC. YOU WANT TO ABB ANY MORE GROUPS ?( Y OR N) ’ !At 
1320 Ir (At = "N" ) GOTO 1370 

1330 IF (At < > "Y") THEN PRINT "WRONG INPUT': GCTC 1310 

1340 INPUT "HOW MANY EXTRA GROUPS ? "5131 

1350 121 = AF1 r 11AF1 = AF1 t 131 

1360 GOTO 1120 

1370 RETURN 

1380 REM 

1390 HOME : PRINT 'LOADS OPTION' 

1400 VTAB (311 PRINT "THE FOLLOWING INPUT PERTAINS TO THE lCAB GARBS.' 

1410 VTAB (6): INPUT "ENTER NUMBER OF LOAB GARBS r ’{All 

1420 HOME 

1430 VTAB (3): PRINT ’ LOAD OPTIONS. " 

1440 VTAB i7>: PRINT TAB' 3)!'i. ONLY RGB ELEMENTS ARE USED." 

1450 PRINT TABC 3>!"2. A MIXTURE OF ROBS ANB BEAM ELEMENTSS ARE USEB. 

1460 VTAB ill): INPUT "ENTER LOAB OPTION ? "PUS 

1470 IF (Ill v 1) OR (Ill > 2) THEN HOME ; PRINT “WRONG OPTION NUMBER! T 
RY AGAIN GOTO 1430 
1430 121 = 1 
1490 GO'SUB 1590 
1500 HOME 

1510 INPUT "BC YOU WANT TO ABB ANY MORE LCAB CARDS Y OS N) r !A$ 

1520 IF (At * -N’ > GOTO 1570 

1530 IF (At < > ' Y* ) THEN PRINT "WRONG INPUT 1 : GOTO I5i0 

1540 INPUT 'HOW MANY EXTRA CARDS?'!131 

1550 121 = All + 1IAI1 = All 4 131 

1560 GOTO 1490 

1570 RETURN 

1580 REM tt*r * m *»m m ?’»>>*»*** **>»****»*« **4** ** 

1590 REM INPUT LOAD-1 
1600 FOR I = 121 TO All 
1610 HOME 

1620 PRINT "LOAD FOR LOAB CARD * il 
1630 VTAB (6) 

1640 INPUT "ENTER NODE NUMBER ? "lAJil.O 
1650 INPUT "ENTER P-X ? ”!AJ(I,1) 

1660 INPUT "ENTER P-Y ? "!AJ(I,2> 

1670 INPUT "ENTER P-Z ? "!AJ(I»3) 

1680 IF (III = 1) THEN HOME I GOTO 1720 
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1690 INPUT 'ENTER N-X ? "»AJ(If4) 

1700 INPUT 'ENTER M-Y ? r rAJ(Ir5> 

1710 INPUT 'ENTER M-Z ? "iAJ(If6> 

1720 NEXT 
1730 RETURN 
1740 REM 

1750 REN SUBROUTINE SAVE DATA 
1760 It* = CHRt <4)1 REN CTRL B 
1770 PRINT IH}" OPEN n iAAt(C) 

1780 PRINT Dti "DELETE" iAAti 0 ) 

1790 PRINT Dti"GPEN iAAt< C S' , L12 ,: 

1800 PRINT Dti’URITE 'SAAt(0-j*S RSI 
1810 PRINT ADI 

1820 PRINT Dtf "WRITE " iAAtC 0 >5" ,R" »2 
1630 PRINT AEZ 

1840 PRINT Dtr’ WRITE "iAAti 0 )i r »R r !3 
1850 PRINT AFZ 

I860 PRINT Dti’WRITE ' 5AAt. C >5' , R r 54 
1870 PRINT AHL 

1880 PRINT Dti' WRITE ‘ iAAt'. 0 )i" ,RS5 

1890 PRINT All 

1900 REN SAVE ELEMENT CARDS 

1910 FCR 1=1 TO AE1 

1920 J * < <I - 1 . » 9 i t 1C 

1930 FOR K = 1 TO 9 

1940 PRINT Dti’WRITE’iAAtiCSSR iTJ + K ) 

1950 PRINT ABZ< I.K) 

I960 NEXT K 

1970 NEXT I 

1980 REN SAVE NFCC 

1990 A = *AEi » 5; + 10 

2000 FOR I = 1 TO ADI 

2010 J=((I-1)*3)+A 

2020 FCR K = 1 TO 3 

2030 PRINT Dti' WRITE r >AAtv 0 if " f R‘ »; J t K) 
2040 PRINT ACi IfN ) 

2050 NEXT K 

2060 NEXT I 

2070 REM SAVE SUPPRESSIONS 

2080 A = (AEZ * 9 j + (AD/. » 3) t 10 

2090 FOR I = 1 TO AHX 

2100 J = Ar«I-l)*7i 
2110 FOR K = 0 TO 6 

2120 PRINT Dti"WRITE’ i AAti 0 )i" fR' i.J + 1 + K> 
213C PRINT AGX( I «K > 

2140 NEXT K 
2150 NEXT I 

2160 REM SAVE MATERIAL i. SECTION PROPERTIES. 
2170 A = A + (AHX * 7) 

2180 FOR 1=1 TO AFZ 
2190 J = A + ((I - 1) * 10) 

2200 FOR K = 0 TO 9 

2210 PRINT Dti "WRITE" iAA*( C SSR'K J il Hi 
2220 PRINT AK(IfK) 

2230 NEXT K 
2240 NEXT I 
2250 REM SAVE LOADS 
2260 A = A + (AFZ * 10) 

2270 FOR I = 1 TO All 
2280 J=A+((I-1>»7) 
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2290 FOR K = C TC 6 

2300 PRINT D*5 ,: WRITE' »AA*< 0 IS* »R C S( J r 1 + K ) 

2310 PRINT AJ(IrK) 

2320 NEXT K 

2330 NEXT I 

2340 PRINT DS T CLOSE'S AAV 0) 

2350 PRINT 'STORAGE OF INPUT DATA COMPLETE" 

2360 RETURN 

2370 rem mmmmmm*mmmmm*»**mumt*******» 

2380 REM SUBROUTINE READ DATA 

2390 D* = CHR$ (4)1 REM CTRL D 

2400 GOSUB 3000 

2410 HOME. 

2420 PRINT D$f"OPEN ’ »AA*C>r r L12 ,: 

2430 PRINT D*? 11 READ " SAAV 0 :5’ »R r i 1 
2440 INPUT ADZ 

2450 PRINT D*S’ READ 1 SAAV 0 H‘ fk' .2 
2460 INPUT AEI. 

2470 PRINT Ii*? - READ ‘ i AAV 0 51” *R* i 3 
2480 INPUT AFZ 

2490 PRINT D* S’READ 1 SAA*< C )S 1 »R’S4 
2500 INPUT AHZ 

2510 PRINT D*r READ * SAAV C tt* *R' 15 

2520 INPUT AIZ 

2530 REM READ ELEMENT CARDS 

2540 FOR I = 1 TC AEZ 

2550 J = (<I - 1) * ?. t 1C 

2560 FGR K = 1 TC 9 

2570 PRINT D$?’ READ " SAAV G )f ’ >R r SiJ + K1 
2580 INPUT ABX(IiK) 

2590 NEXT K 

2600 NEXT I 

2ol0 REM READ NPCC 

2620 A = (AEZ * 9i + 1C 

2630 FOR: I = 1 TC ADZ 

2640 J * < (I - 1> * 3) t A 

2650 FOR K = 1 TO 3 

2660 PRINT D*S’ READ ’’.AAV 0 ii' rR* K J r U 5 
2670 INPUT AD I.K) 

2680 NEXT t; 

2690 NEXT I 

2700 REM READ SUPPRESSIONS 
2710 A = (AEZ * 9) + (ADZ * 3) + 1C 
2720 FOR I * 1 TO AHZ 
2730 J = A + UI - 1) * 7) 

2740 FOR K = C TC 6 

2750 PRINT D$ S ’’ READ “ SAAV CIS' »R’ Si J + I + K, 

2760 INPUT AGZi 1,10 
2770 NEXT K 
2780 NEXT I 

2790 REM READ MAT. t SECT. PROP. 

2800 A = A + (AHZ * 7) 

2810 FOR I » 1 TC AFZ 
2820 J = A + ((I - 1) * 10) 

2830 FOR K = 0 TC 9 

2840 PRINT D$S"READ " SAAV C )>' »R' Si J r 1 + K ) 

2850 INPUT AW 1,0 
2860 NEXT K 
2870 NEXT I 
2880 REM READ LOADS 







»;’ ,R r »t j + i + k/ 


2890 A = A t (AF% * 1C) 

290C FOR I = 1 TO AI% 

2910 J = A + <(I - I > * 7) 

292C FOR K - G TO 6 
2930 PRINT Stf"READ ,: ;AAt(C) 

2940 INPOT AC(I*K) 

2950 NEXT K 

2960 NEXT I 

2970 PRINT UtrCLOSE' JAA$( 0) 

298C RETURN 
2990 REM 

3C0C REN INPUT HEADING 
3010 HONE I VTAB (3, 

3C20 PRI»Tt "SUBROUTINE HEADING' 

3C3C VTAB (10) 

3040 INPUT "DUG ID CODE '?' }»Ai(Ci 
3050 INPUT "DATE V iAAIDtd ) 

3060 INPUT "TOUR NAHE ?"iAAIDt(2. 

3100 RETURN 

3370 REN #»#»>*»**» ************ »»**#* 

3380 REN INPUT ELEMENT CARDS 
3390 INPUT "NUMBER OF ELEMENT GENERATING CARD 
3400 12% = 1 

3410 FOR II = 12% TC AE% 

3420 HONE 

3430 PRINT "ENTER DATA FOR E~M T CARD # Ml 
3440 PRINT 

3450 VTAB < 5 It INPUT "1. ELEMENT TYPE » GET Cut 
3460 INPUT "2. GROUP NUMBER."?ABi.11»2 > 

34 7 0 INPUT "3. NUMBER OF ELEMENTS Il, 

3480 INPUT "4. NODE #1.">AB%(Ilf4 ) 

3490 IF AB%( 11*3) = 1 GOTO 3510 

3500 INPUT "5. INCREMENT#! V #A3%( IfcrS: 

3510 INPUT "6. NODE #2. r"tA8%(11*6) 

3520 IF AB%(I1»3) = 1 GCTC 3540 

3530 INPUT "7. INCREMENT #2 V JAB%CI1.7; 

3540 IF AB%( 11; 1) < 3 GOTO 3580; REM *1 

3550 INPUT "S. NODE *3 V ;AB%( 11*6> 

3560 IF AB%(I1»3) = 1 GOTO 3530 

3570 INPUT "9. INCREMENT *3 ?"iAB%(II t9 > 

3580 NEXT 
3590 HOME 

3600 INPUT "DC YOU WANT TC ADD MORE ELEMENT 0 
3610 IF (At = "N" ) GCTC 3670 
3620 IF (At < | "Y"> THEN PRINT "URCNG INFU 
3630 HOME ; INPUT "riCL MANY ? "M3% 

3640 12% * AE% + 1 
3650 AE% = AE% r 13% 

3660 GOTO 3410 
3670 RETURN 

3680 ren 

3690 HOME I REN INPUT OF NPCC 

3700 INPUT "ENTER MAX. NUMBER CF NODES (N<=5C 

3710 12% = 1 

3720 FOR II = 12% TC AD% 

3730 HOME : PRINT "ENTER THE NPCC FCR NODE * 
3740 VTAB (5)! INPUT "X COORDINATE r' JACdlrl 
3750 INPUT "Y COORDINATE ?"iAC(Il,2> 

3760 INPUT "2 COORDINATE ?'}AD. II .3 ) 

3770 NEXT 
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3780 HOME 

379C INPUT "DC YOU WANT TO ADD ANY MORE NODES <Y DR N > ? r 5A* 

380C IF (A* = -N-) GOTO 385C 

3810 IF (At < > "Y r ) THEN PRINT "WRONG INPUT 1 ; GCTC 3750 
3820 HOME I INPUT ''HOW MANY ? "513% 

383C 12% = ADJ. + 1JAD2 = ADI + 13;. 

3840 GOTO 3720 
3850 RETURN 

3860 reh mmmmmmmmmmmm*m*mmmmmm 
3870 HOME t PRINT "LIST OPTION" 

3880 VTAB <5K PRINT "TURN PRINTER ON' 

3890 GOSUB 23801 REM READ DATA 
3900 REM _ PRINTER CONTROL 
3910 P5* = ’"I REM LINE LENGTH 
3920 P6t = REM NORMAL MODE 
3930 P9t = " “ t KEH 10 C;HAk/t. INF 
3940 PA* = ""1 REM 12 CHAR/LINE 
3950 PB* * REM SELECT PRINTER 
3960 PR# 1 
3970 PRINT P5* 

398C PRINT P6t 
3990 PRINT PA* 

400C PRINT PB* 

4010 PRINT "HEADING INFORMATION" 

4020 PRINT t PRINT 

4030 PRINT "DWG ID CODE IS .... "lAAt(O) 

4050 PRINT "DATE IS ."»AA*( 1 ) 

4060 PRINT "PREPARED BY .">AA*(2; 

408C REM PRINT ELEMENTS 

4090 PRINT J PRINT 

4100 PRINT "ELEMENT CARDS" I PRINT 

4110 PRINT " LINE"} 

4120 HTAB (8): PRINT “ELMT GROUP * OF" 5 

4130 HTAB (27): PRINT "NODE INC NODE INC NODE INC 1 

4140 PRINT " 

4150 PRINT TAB< 4)4" TYPE"» TAB: 4);"*"; TAB: 6)1 "ELMT 1 5 TAB: 5)5 
4160 PRINT ”1*'5 TAB( 5) 5 * 1 c 5 TAB( 6)5"2 ; i TAB( 6>5’2 : ; 

4170 PRINT TABi 4)5 r 3"5 TAB( 4)5"3 : : PRINT 
4180 FOR I * 1 TC AE2 
4190 HTAB (4): PRINT I. 

4200 HTAB (9)1 PRINT AB%(I,i>5 

4210 HTAB (16): PRINT AB%( 1,2)5 

4220 HTAB (22): PRINT ABC: 1,3)5 

4230 HTAB (28): PRINT AB%< 1,4)5 

4240 HTAB (33): PRINT AB;.( 1,5)5 

4250 HTAB (39): PRINT AB%( 1,6)5 

4260 HTAB (1J: PRINT TAB: 6»»AB1CI,7*I 

4270 HTAB (1): PRINT TAB: 6)5 ABU 1,8)5 

4280 HTAB (1).‘ PRINT TAB: 6>5ABUI,9) 

4290 NEXT 

4300 REM PRINT NPCO 
4310 PRINT J PRINT 

4320 PRINT "NODE #* 5 TAB) 3 )5 ,: X-C00RD k 5 TAB: 4 )5" Y-COORD' 5 TAB. 4 >5’2-COO 

RD" 

4330 PRINT 

4340 FOR I = 1 TO ADD 
4345 PRINT IJ 

435C HTAB (UK PRINT AC( 1,1)5 
4360 HTAB (21): PRINT AC( 1,2)5 
4370 HTAB (31): PRINT AC: 1,3) 
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NEXT 

REM PRINT NODAL SUPPRESSIONS 
PRINT : PRINT 

PRINT “LINE 4 NODE *“ 5 TAB; 3 >5 
PRINT “X Y Z 
PRINT "RX RY RZ' 

PRINT 

FOR I = 1 TO AH2 
HTAB (2)! PRINT 15 
HTAB (8): PRINT AG%( 1.0)5 
HTAB (17)J PRINT AGZ; I»1)» 

HTAB (21)*. PRINT AG%< 1,2)5 
HTAB (26): PRINT AGZ( 1.3)5 
HTAlT (30): PRINT AG%(I,4>5 
HTAB (35)’. PRINT AGX(If5» 

PRINT TAB; 5>5AG3i; I.e) 

NEXT 

PRINT : PRINT 

REM PRINT MATERIAL PROPERTIES 
PRINT “MATERIAL PROPERTIES' 

PRINT 

FOR I = 1 TC AF1 STEP 3 

11 = I + 1 

12 = I + 2 

PRINT “LINE HTAB (8)1 PRINT I.II.12 
PRINT “GROUP"*: HTAB (8 PRINT 1*11*12 
PRINT “cLKT 1 i I HTAB ,3;; PRINT AK( I*C )*AK( 11*0 )*AKi I2*C ) 

PRINT • AREA r * I HTAB (8)1 PRINT AK< I *1 )*AK( 11*1 )*AK( 12*1) 

PRINT ’ I-YY' * I HTAB (3)1 PRINT AK( 1,2 )*AK( 11,2 *>AK( 12*2 ; 

PRINT "I-ZZ r ;: HTAB (8K PRINT AK( I*3)*AK(I1*3,*AK; 12*3) 

PRINT “E'5 5 HTAB (8K PRINT AK( 1*4 )»AK; 11*4 )*AK( 12*4 ) 

PRINT “V“*: HTAB (8): PRINT AK; I *5 )*AK( 11,5 >*AK; 12.5 ; 

PRINT “J“*: HTAB (8*: PRINT AK( 1,6 ;.AK( 11*6 ;.AK( 12,6 ) 

PRINT "X' “ *: HTAB (8)1 PRINT AK( I *7 ,»AM 11*7 )*Af.-. 12,7 > 

PRINT ■Y'"*: HTAB (8)J PRINT AK; 1 *8 )*AKi 11*3 .*,AK; 12,3 

PRINT "Z'- *: HTAB (8K PRINT AK( 1*9 ;*AK( 11*9 ;»AK( 12,9 ) 

PRINT : PRINT 
NEXT 

REM PRINT “LOADS" 

PRINT : PRINT “LOADS" 

PRINT “LINE * NODE ♦“ * TABt 3)5 
PRINT “PX PY PZ “5 

PRINT "MX MY MZ“ 

FOR I * 1 TO AI2 
HTAB (2): PRINT I* 

HTAB (Si! PRINT AJ< 1*0), 

HTAB (15): PRINT AJ(I,1>5 
HTAB (21): PRINT AJ( 1*2)5 
HTAB (28): PRINT AJ( 1,3)5 
HTAB (35)1 PRINT AJ< 1*4)5 
PRINT TAB( 12 )5AJ< 1*5)5 
PRINT TAB( 12 >5AJ( 1*6 ) 

NEXT 
PR* 1 
RETURN 

rem mmmmmmmmm**»*»»m***m***»*» 

REM CCFY CPTICN 

HOME : PRINT "INSERT FLOPPY CONTAINING DATA TO BE COPIED' 
PRINT I FRINT "HIT RETURN WHEN READY TO PROCEED' 

INPUT A* 
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4955 GOSUB 2380t REM REAL BATA 

4960 HOME I PRINT "INSERT FLOPPY TC BE COPIED INTO 1 

4970 PRINT ! PRINT "HIT RETURN WHEN READY TC PROCEED' 

4980 INPUT A* 

50CC GOSUB 1750t RErt SAVE DATA 
5C10 RETURN 

5020 ren mm**mt**m**m*m*m*u**m ************ 

5030 GOSUB 2390: REN READ BATA 

5050 HONE I INVERSE : PRINT "CHANGE OPTION": PRINT I NORN.AL 
5055 PRINT : PRINT "THIS ARE THE OPTIONS AVAILABLE!" 

5060 PRINT ! PRINT TAB: 5>i"l. ELEMENT GENERATOR CARDS." 

5070 PRINT TAB< 5)!"2. NPCDS." 

5080 PRINT TAB' 5){'3. SUPPRESSIONS." 

5090 PRINT TAB; 5)5"4. MATERIAL AND SECTION PROP.' 

5100 PRINT TAB( 5)5“5. LOADS." 

5110 PRINT TAB: 5)5"6. EXIT." 

5120 PRINT ! INPUT "ENTER OPTION NUMBER ?'5A2 

5130 IF (AZ < 1) OR (A* > 6) THEN HOME J PRINT ' WRONG OPTION, TRY AGAIN. 

"! GOTO 300 

5140 IF AS = 1 THEN GOSUB 5230 

5150 IF AZ = 2 THEN GCSUB 555C 

5160 IF AZ = 3 THEN GCSUB 5S0C 

5170 IF AZ = 4 THEN GOSUB 6150 

5180 IF AZ = 5 THEN GOSUB 6570 

5182 IF AS. < > 6 THEN GOTO 5050 

5185 GOSUB 1750! REM SAVE BATA 
5190 RETURN 

5210 rem mmmm»*m»*m*»**»***»**»»***»»*»**»*»******** 

5220 HOME 

5230 HOME I INVERSE ! PRINT "CHANGE ELEMENT GENERATOR CARDS' : NORr.AL 
5240 PRINT ! PRINT "THIS ARE THE OPTION AVAILABLE:' 

5250 PRINT I PRINT TAB: 5)5"1. CHANGE A LINE." 

5260 PRINT TAB( 5)5"2. ABB OR DELETE A LINE.' 

5270 PRINT TAB: 5>i"3. EXIT." 

5280 PRINT I INPUT "ENTER OPTION NUMBER ? " i AS. 

5290 IF (AX < 1 ) OR (AZ > 3) THEN HOME : PRINT 'WRONG OPTION, TRY AGAIN. 

°: GOTO 5240 

53OC IF AZ « 1 THEN INPUT "ENTER NUMBER OF LINE TC BE CHANGED ? ‘ill; GOTO 

5360 

5310 IF AZ = 3 THEN GOTO 553C 

5320 INPUT "DC YOU WANT TO 'ADD' OF. 'DELETE' A LINE ? '{AS 
5330 IF At * "ABB" THEN AES. = AEZ + 1111 = AEZ: GCTG 53s: 

5340 IF At = "DELETE- THEN HOME !AEZ = AES. - 1! PRINT “LAST LINE HAS BEE 
N DELETED.": GOTC 5230 

5350 HOME ! PRINT "WRONG KEYWORD HAS BEEN ENTERED, TRY AGAIN. : GCTG 5230 
5360 HOME 

5370 PRINT -ENTER BATA FOR ELMT CARD * 'ill 
5380 PRINT 

5390 VTAB <5>! INPUT "1. ELEMENT TYPE (OPTION *) ?" iAfcJti 11,1; 

5400 INPUT "2. GROUP NUMBER."iABZ<11,2 i 
5410 INPUT ’3. NUMBER OF ELEMENTS ?"iABZUl,3S 
5420 INPUT -4. NODE ♦1."}ABZ(11,4 ) 

54 30 IF ABZ( 11,3) = 1 GOTO 5450 

5440 INPUT "5. INCREMENT *1 ?"iABZUl,5) 

5450 INPUT '6. NODE *2. ?'iABZ;ilr6) 

5460 IF ABZ;il,3) = 1 GOTO 5480 

5470 INPUT "7. INCREMENT *2 ?"fABZ(Il,7) 

5480 IF AB%( 11,1 ) < > 3 GOTO 5530: REM 41 
54 90 INPUT -8. NODE *3 ?"JARZ(I1,8) 
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5500 IF ABZ( 11*3) = 1 GOTO 5530 

5510 INPUT "9. INCREMENT #3 ?'iABX(Il,9) 

5530 RETURN 

5540 ren m**mm***m*m*M 

5550 HOME 

5560 HTAB (8)1 INVERSE I PRINT ’CHANGE NPCCS"NORMAL 
5570 PRINT PRINT "THIS ARE THE OPTION AVAILABLE'.’ 

5580 PRINT : PRINT TAB( 5)5" 1. CHANGE A LINE.’ 

5590 PRINT TAB' 5)?”2. ADD OR DElETE A NODE.’ 

5600 PRINT TAB( 5>5"3. EXIT." 

5610 PRINT I INPUT "ENTER OPTION NUMBER ? ’iAX 

5620 IF (AX < 1 ) OR (AX > 3) THEN HOME I PRINT ’ WRCNG OPTION, TRY AGAIN. 

*: GftTO 5570 

5630 IF AX * 1 THEN HOME ! INPUT ’’CHANGE COORDINATES OF NODE * ? “Ill: GOTO 

5730 

5640 IF A% = 3 THEN GOTO 5780 

5650 INPUT "DC YOU WANT TO 'ADD' OR 'DELETE' A NODE ? ’ »A* 

5660 IF A$ = "ADD’ THEN ADI = ALL + II11 = ALXI GCTG 5730 

5670 IF (A* < > "DELETE”) THEN HOME : PRINT "WRONG KEYWORD HAS BEEi, ENT 

ERED, TRY AGAIN.": GOTO 5560 

5680 HOME : INPUT “ENTER NUMBER CF NODE TO BE DELETED ? ’»Ii 
5690 AC( 11*1) = 0 
5700 ACi 11*2) = 0 
5710 AC< 11*3> = 0 

5720 HOME 5 PRINT "NODE HAS’BEEN DELETED.’ : PRINT I GCTC 5560 

5730 HOME : PRINT "ENTER THE NPCO FOR NODE ♦ "fll 

5740 VTAB <5i! INPUT ‘X COORDINATE ?’5AC. 11,1) 

5750 INPUT "Y COORDINATE 7“ *AC< 11*2) 

5760 INPUT "Z COORDINATE 7“ *ACL II *3 > 

5770 HOME : GOTO 5550 
5780 RETURN 

5790 rem *mm***mmm***mm*mm**»*»*m»m**u****w* 

5800 rem mm ****** 

5810 HOME 

5820 INVERSE : PRINT "CHANGE NODAL SUPPRESSIONS’ : NORMAL 
5830 PRINT ! PRINT "THIS ARE THE OPTION AVAILAB-E: ’ 

5840 PRINT ! PRINT TAB( 5>»" 1. CHANGE A LINE.’ 

5850 PRINT TAB( 5)!"2. ADD OR DELETE A LINE.' 

5860 PRINT TAB' 5>?"3. EXIT.” 

5870 PRINT I INPUT "ENTER OPTICN NUMBER ? * »AX. 

5680 IF (AX < 1) OR (AX > 3) THEN HOME : PRINT 'WRONG OPTICN, TRY AGAIN. 

": GOTO 5830 

5890 IF AX = 1 THEN INPUT "ENTER NUMBER OF LINE TO BE CHANGED 7 *1: GCTG 
5950 

5900 IF AX = 3 THEN GOTO 6130 

5910 INPUT "DC YOU WANT TO ADD' OR 'DELETE' A LINE ? ’ iA* 

5920 IF A* = "ADD" THEN AHX = AHX + KI = AHi: GOTO 5950 

5930 IF A* = "DELETE" THEN HOME I AHX. = AHX - 1*. PRINT ’ LAST LINE HAS BEE 
N DELETED.": GOTO 5820 

5940 HOME : PRINT “WRONG KEYWORD HAS BEEN ENTERED, TRY AGAIN.’ : GCTC 5820 

5950 HOME I PRINT "ENTER NODE NUMBER AND Y OR N TO THE QUESTIONS.’: VTAB 

(3) 

5960 PRINT 

5970 PRINT "NCDAL SUPPRESSIONS CARD # ... ")I 
5980 VTAB (6) 

5990 INPUT "NODE NUMBER ? "iAGX(I,0> 

5992 FOR J * 1 TO 6 
5994 AGX(I.J) = 0 
5996 NEXT J 
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60OC INPUT “SUPPRESS X ? "fit 
6010 IF (1$ = "Y“ ) THEN AG2(I,1> = 1 
4020 INPUT "SUPPRESS Y ? “lit 
6030 IF (1$ = "Y* ) THEN AGX!I,2) = 1 
4040 INPUT “SUPPRESS Z ? "lit 
6050 IF (It = “Y*) THEN AEi< I»3> = 1 

6060 IF III = 1 THEN AGZ(I,4) = 11AGX(I*5> = 11AG2(I,6> = 11 GOTO 6130 

6070 INPUT “SUPPRESS THETA X ? “lit 

6080 IF (It = "Y* ) THEN AGX(I,4) = 1 

6090 INPUT “SUPPRESS THETA Y 7 'lit 

6100 IF (It = “Y“ ) THEN AG£(i,5) = 1 

6110 INPUT “SUPPRESS THETA I ? 'lit 

6120 IF (J .t = "Y“ ) THEN AG2(I,6) = 1 

6125 GOTO 5810 

6130 RETURN 

6140 REN 

6150 REN *t***6*******»*}6;»}6}.*M*}4 4}*6*66****6***4**66»64»6*4.}. 

4160 HOME 

617C INVERSE 1 PRINT “CHANGE MATERIAL t, SECTION PROP.' ; NORMAL 
6180 PRINT I PRINT “THIS ARE THE OPTION AVAILABLE 1' 

6190 PRINT I PRINT TAB; 5>1 1. CHANGE A LINE.' 

6200 PRINT TAB; 5H"2. ADI OR DELETE A .INEZ 

4210 PRINT TAB; 5)1 "3. EXIT.' 1 

6220 PRINT 1 INPUT "ENTER OPTION NUMBER ? 1 »A2 

6230 IF (A2 < 1) OF; ; fii > 3) THEN HOME i PRINT ' WRONG OPTION» TRY AGAIN. 
"1 GOTO 6180 

6240 IF A/, r 1 THEN INPUT “ENTER NUMBER Or LINE TO BE CHANGED ? ‘ lit GOTO 
6310 

625C IF A2 = 3 THEN GOTO 6550 

6260 INPUT “DC YOU WANT TO ADI' OR DELETE A LINE ? ' 1A* 

6270 IF At = “ADD' THEN AFi = AF2 + ill = AF21 GOTO 6310 

6280 IF At = "DELETE" THEN HOME 1AF2 r AFL - 11 PRINT 'LAST lINE HAS BEE 
N DELETED."1 GCTC 6170 

6290 HOME 1 PRINT “WRONG KEYWORD HAS BEEN ENTERED, TRY AGAIN. 1 GCTC oi/v 

6300 HOME 
6310 HOME 

6320 VTAB (3)1 PRINT “ENTER OPTION NUMBER' 

6330 VTAB (7)1 PRINT TAB; 3)S"1. ONLY ROD ELEMENTS ARE USED.' 

6340 PRINT TAB; 3)1*2. A MIXTURE OF RODS AND BEAM ELEMENTS ARE USED. 

6350 VTAB (11 ) 

6360 INPUT "ENTER OPTICN NUMBER I'lIlS 

6370 IF (112 < 1) OR (112 > 2) THEN HOME 1 PRINT “WRONG OPTION NUMBER 1 7 
RY AGAIN “1 VTAB (3)1 GOTO 6310 
6330 HOME 

6390 PRINT “MAT. 6 SECT. PROPERTIES FOR GROUP * II 
64 00 IF 112 = 1 THEN AM; 1,0) = 11 GOTO 6420 
6410 INPUT “ELEMENT TYPE ? “1AK; 1,0) 

64 20 VTAB (3)1 INPUT “CROSS-SECTIONAL AREA ? r 1 AN; 1,1) 

6430 IF 112 = 1 GOTO 6470 

6440 PRINT “MOMENTS OF INERTIA ABOUT LOCAL AXIS' 1 PRINT 
64 50 INPUT “I <Y'-Y' i ? “!AK( 1,2) 

6460 INPUT “I il'-Z- ) ? “!AK( 1,3) 

44 70 PRINT 1 INPUT “MODULUS OF ELASTICITY E ? “5AK(I,4> 

6430 IF 112 = 1 GOTO 6160 

64 90 INPUT "POISSON'S RATIO V ? “!AK(I,5) 

6500 INPUT “POLAR MOMENT OF INERTIA .J ? “1AK(I,6) 

6510 INPUT “X' ? “ !AKi 1,7) 

6520 INPUT “Y ? *»AK>. 1,6) 

4530 INPUT *Z' ? “»AK( 1,9 ) 
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6540 GOTO 6160 
6550 RETURN 

6560 rem m*mmmm*m**m*m**nm 

6570 rem m**»4*****»»»»****»y*m»* «**»»*****»*».»**** **»»***» 

6580 HOME 

6590 INVERSE t PRINT "CHANGE LOAD CARTS 1 NORMAL 

6600 PRINT : PRINT "THIS ARE THE OPTION AVAILABLE:" 

6610 PRINT : PRINT TAB( 5)5*1. CHANGE A LINE.' 

6620 PRINT TAB( 5)5*2. ADD OR DELETE A LINE." 

6630 PRINT TABC 5 >5’'3. EXIT." 

6640 PRINT 5 INFUT "ENTER OPTICN NUMBER ? "5AX 

6650 IF (At < 1) OR (AX > 3) THEN HOME t PRINT "WRONG OPTION. TRY AGAIN. 

•1 GOTO 6600 

6660 IF' AJ. * 1 THEN INPUT "ENTER NUMBER CF LINE TC BE CHANGED ? "51: GOTO 

6720 

6670 IF At * 3 THEN GDTG 6900 

6680 INPUT "DC YOU WANT TO ADD' OR 'DELETE A LINE ? ' 5 At 

66 90 IF At = ”ADD* THEN Alt = AIX + ill * AU1 GCTC 021 

6700 IF At = "DELETE* THEN HOME IAIt = AIX - i: PRINT "LAST „INE HAS BEE 
N DELETED.": GOTO 6590 

6710 HOME : PRINT "WRONG KEYWORD HAS BEEN ENTERED. TRT AGAIN.' ; GOTO 6590 

6720 HOME 

6730 VTAB (3): PRINT "LOAD OPTIONS. " 

67 40 VTAB (7)! PRINT TAB. 3)5*1. ONLY ROD ELEMENTS ARE USED. 

6750 PRINT TAB< 3)5*2. A MIXTURE GF RODS AND BEAM ElEMENTSS ARE USED. 

6760 VTAB (ID! INPUT "ENTER LOAD OPTICN ? *5117. 

6770 IF (I1X < 1) 0F; (112 > 2) THEN HOME : PRINT WRONG OPTICN NUMBER. T 

RY AGAIN ”1 GOTO 6730 
6780 HOME 

6790 PRINT 'LOAD FGR LOAD CARD * *51 
6800 VTAB (6) 

6810 INPUT "ENTER NODE NUMBER ? *5AJ(I.C) 

6820 INPUT "ENTER F-X ? "5AJ(I.D 
6830 INPUT "ENTER P-Y ? *5AJ;i>2) 

6840 INPUT "ENTER P-Z ? *5AJ(I»3> 

6850 IF (I1X = 1) THEN HOME : GOTO 1720 
6860 INPUT "ENTER M-X T 'fftJl1.4 ) 

6870 INPUT "ENTER M-Y ? "5AJCI.5,' 

6880 INPUT "ENTER M-Z ? "5AJU.6) 

6890 GOTO 6580 
6900 RETURN 

6910 REM **».***'***» **»*»»»»**«»»»»»»***« < * 
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5 HOME : CLEAR 

10 HTAB (10)1 PRINT 'PROGRAM ASTRAl'I PRINT 
20 REM 25 NOVEMBER 1980 

60 HTAB (5): PRINT "ANALISING STRUCTURES WITH APPLE" 

80 GOSUB 1710t REM ALLOCATE ARRAYS 
90 GOSUB 21501 REM READ BATA 
110 GOSUB 57801 REM PRINT INPUT DATA 
155 GOSUB 8610: REM CALCULATE BANDWIDTH 
158 REM. GOSUB 8750; REM STORE BANDWIDTH CF 

160 GOSUB 7860: REM DELETE OLD SSM AND CREATE A NEW NUlL (ALL 0) S5M RE 
GISTER 

200 BP = ADD * 6i REM SIZE OF SSM BP*BP 

210 PRINT ' SIZE OF SSM IS r JBP 

220 FOR 18 = 1 TO AEZS REM * ELEMENT CARD LOOP 

230 GOSUB 3530: REM SECT PROP. 

231 FOR 19 = 1 TO ABZ( 18.3 K REM ELEMT/CARD LOOP 
233 GOSUB 250: REM CALCULATE STRUCTURE SSM 

236 NEXT 19 

235 NEXT 18 

236 PRINT 1 PRINT 'ASSY OF SSM IS COMPLETE": PRINT 

237 GOSUB 731C: REM ASSY OF SS IN D0FIJ:I*2) 

238 PRINT 'SS HAS BEEN STORED IN DOFIJ' 

239 GOSUB 7620: REM STORE LOAD IN DOFIU, I>i) 

260 PRINT "P HAVE BEEN STORED IN DOFIU 

261 GOSUB 3380: REM SAVE DOFIU IN DISK 

262 REM GOSUB 3060:REM PRINT DOFIU 

263 HOME 

266 PRINT "PROGRAM ASTRA! IS COMPLETED. TO CONTINUE LOAD AND RUN PROGRAM 
ASTRA2” 

265 END 

266 REM m**#**.i****)i.**i*>**)r*'**6J»*6 »»»*»* »*>»*****« ***»*< 

268 REM CALCULATION OF SSM 

250 BRZ = BRZ + i: REM ELEMENT COUNTER 

260 GOSUB I860: REM CM FOR VERTICAL ELEMENTS.OBTAIN Nl.KC.NODAL COORD 
S., DIRECTION COSINES. 

265 PRINT : PRINT "ELEMENT * " JBRZS SPEC 2 XT'NODE 1= 'Ifiii 5PC. 2). NODE 

J= r 5N2 

270 GOSUB 2060: REM ACTIVATE D.G.F. IN USE 

280 PRINT ! PRINT "THE DIRECTION COSINES CX.CY.CZ FOR ELEMENT .BRZi ARE 


290 PRINT CXi SPC( 3)9CY? SPC< 3)*CZ 

295 IF (ABZ(IS.l) = 2) THEN GOSUB 60041 REM CALCULATE CG.SG 

300 GOSUB 3260: REM CLEAR MATRICES 

310 REM CALCULATE TRANSFORMATION MATRIX 

320 IF (ABZ(I8.1) = 1) AND (BIZ = 1) THEN GOSUB 5170: REr, VERTICAL ROD 
330 IF (ABZ( 18.1 ) = 1) AND ( BIZ = 0) THEN GOSUB 50:0: REM NON-VERTICAL 


360 IF (ABZ( 18.1) = 2) AND ( BIZ = 1) THEN GOSUB 6210: REr, VERTICAL l 
EAM 

350 IF < ABZ( 18.1 ) = 2) AND (BIZ = 0) THEN GOSUB 4520; REM NOr. VERTICAL 

BEAM 

360 GOSUB 7660: REM STORE ETM 

370 REM GOSUB 6870: REM PRINT TRANS MATRIX 

380 GOSUB 5260: REM FIND INVERSE (TRANSPOSE) OF TRANSFORMATION MATRIX. 
390 REM GOSUB 69601 REM PRINT TRANSPOSE OF ELMT TRANS MATRIX 
600 REM ELEMENT STIFFNESS MATRIX (LOCAL) 
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410 

420 

430 

440 

445 

450 

460 

470 

480 


IF AB2< ISt-i.l = 1 THEN G03UE 5340: REM SPACE TRUSS 

IF ABD<I8*1> = 2 THEN GOSUB 2870: REN SFACE BEAM 

PRINT 

IF AB5X 18,1) = 1 THEN GOSUB 8190: REN STORE ESN ROB 

IF ABKIB.l) = 2 THEN GOSUB 8800: REN STORE ESN BEaM 

REM GOSUB 7050: REM PRINT STIFF MATRIX 

GOSUB 5410: REM ELEMENT STIFFNESS MATRIX < GlOBAl ) 

REM GOSUB 7130: REM PRINT ELHT STIFF MATRIX 

REM 


>4*4# 


490 REM ADD ELEMENTS STIFFNESS TO STRUCTURE STIFFNESS MATRIX. 

500 PR INI. "ASSY OF STRUCTURE STIFFNESS MATRIX 

510 REM 666*666646*6664664444444 644646**# *444# 4# 4*#4.*4***#4#*##4 ***#*# # 

66666666666466616666644464446664464644464i44444 
520 REM R2S IS THE ADDRESS WHERE THE UAlUE IN SSr. UI_L BE STORED 
530 D$ = CHRt <4) 

535 A* = "SSMIX' 

540 PRINT Pt* 1 'OPEN 1 9' * L20’ 

550 FOR I * 1 TC 12 
560 FOR J = 1 TO I 
565 BO = 0 

570 IF BH( I»J ) = 0 GOTO 310 

58C GOSUB 8700: REM CALCULATE ADDRESS CF SSM IN DISN 
680 PRINT [Iti 1 READ ’SAt?', R' *R2» 

690 INPUT BG 

700 BG = BG + BH: I,Ji 


780 PRINT Dt»' WRITE . R' *R25i 

790 PRINT BG 
810 NEXT J 
820 NEXT I 

830 PRINT D$* r CLOSE ’!A$ 

1540 PRINT ''ELEMENT « r * BRD* 1 STIFFNESS MAS BEEN ADDED 70 SSM' 

1545 RETURN 

1700 REM 6444*****44 4 44 *#*** 1444 *** 4 *** 444 * 44444*4 * 4 ***** 4 ************« 

6664 *** 4 ** 4 * 44 *************#** 44**4 ********44********************<**< 


666644644666*44**4* 

REM ALLOCATION OF ARRAY? 
AAIB*< 1 )»ABEGEND>:25.9> 


1710 

1720 DIM 
1730 DIM ACNPCG( 50»3 > 
1740 DIM AGNBOFIX 50*6 * 
1750 DIM AJPNODE:50 *e ) 
1760 DIM AKMSPROP: 10*9) 
1770 DIM D0FIUC300*3) 
1780 DIM BB( 12r 12 )I REM 
1790 DIM BC>, 12*12): REM 
1800 DIM BF< 12.12): REM 
1610 DIM BH( 12*12)1 REM 
1820 DIM BT( 12*12 >: REM 
1830 DIM BH( 3> 3)* BN( 3*3 
A 


TnAmS. MATRIX 

stiff, matrix cerement 

TRANS. MATRIX INVERSE 
STIFF. MATRIX (Ei-EMENT- 

bt=bf4b: used to ca.cu,. 
*BC-. 3*3)* DX\ 3 )* Ex: 3 ■: RE 


1840 RETURN 

1850 REM 6666644444*46644444444444444444444*4. 
1860 REM CK FOR VERTICAL ELEMENTS. OBTAIN Nl.NI 
1870 BID = 0 

1880 N1 = ABDC 16*4) + (ABDCI8.5) 6 >19 - 1). 

1890 N2 = ABXafi.6) + CASK 18.7) 6 <19 - 1); 

1900 IF N1 > N2 THEN N3 = N21N2 = NllNl = N3:N3 
1920 XI = AC< N1 * 1 ) 

1930 Y1 = AC( N1 *2 ) 

1940 Z1 = AC< N1 ? 3 ) 
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1950 X2 = ACt N2,1) 

I960 Y2 = AC(N2»2) 

1970 Z2 = AC(N2,3) 

1980 REM CK FOR VERTICAL ELEMENTS 

1990 biz = 0 

200C IF ( ABS (X2 - XU > .00001 ) GOTO 2010 

2005 IF ( ABS <Z2 - Zi) > .00001 ) GOTO 2010 

2006 BIZ = 1 

2010 REM CALCULATE DIRECTION COSINES 

2012 NL = < (X2 - XI) T 2 + (Y2 - Y1) t 2 + (Z2 - ZU t 2) t .5 

2013 CX = <X2 - XI) / NL 

2014 CY = < Y2 - Yi) / NL 

2015 CZ = (Z2 - Zl i / NL 

2020 RETURN 

2030 REM 

2040 REM ACTIVATING N.D.F. USED 
2050 K = (N1 - 1) * 6 
2060 K3 = (N2 - 1) * 6 
2070 K2 = 3 

2080 IF ABU IS-1 ) = 2 THEN K2 = 6 

2090 FDR I = 1 TO K2 

2100 DOFIU' v K + I )»0) = 1 

2110 DQFILK (K3 + I >«0) = 1 

2120 NEXT I 

2130 RETURN 

2140 REM ttmityiMmiiMOHnniiMtmmnuimiMMitiio 

2150 REM SUBROUTINE READ DATA 

2160 D* = CHRS (4): REM CTRL D 

2161 PRINT : INPUT "ENTER DRAWING CODE NUMBER ? "SAAVO) 

2162 PRINT D$i"OFEN ID r 

2163 PRINT DS?"DELETE ID 

2164 PRINT In> 1 OPEN ID, L12' 

2165 PRINT DS»"WRITE ID,F.l’ 

2166 PRINT AAV 0 ) 

2167 PRINT Dti"CLOSE ID" 

2180 PRINT t PRINT "READING INPUT INFORMATION' 

2190 PRINT D*f 1 OPEN "JAAVCf , i_12‘ 

2200 PRINT DtUREAD "JAAVOU.R 51 
2210 INPUT ADZ 

2220 PRINT Dt J 1 READ " JAAV 0 )?’,R'52 
2230 INPUT AEZ 

2240 PRINT inf'READ ' tAA$i C ) i' *R" 5 3 
2250 INPUT AFZ 

2260 PRINT D*T READ " i AAV 0 m " ,R" i4 
2270 INPUT AHZ 

2280 PRINT D*>"READ ' JAAV 0 )}’ ,R' 55 

2290 INPUT AIZ 

2300 REM READ ELEMENT CARDS 

2310 FOR 1=1 TO AEZ 

2320 J«((I-1)*9)+10 

2330 FOR K = 1 TO 9 

2340 PRINT Dt?"READ "»AA$:0 )»R"r(J + Ki 

2350 INPUT ABZ( I,Ni 

2360 NEXT K 

2370 NEXT I 

2380 REM READ NPCO 

2390 A = (AEZ * 9) + 10 

2400 FOR I * 1 TO ADZ 

2410 J = <<I - 1) * 3) + A 

2420 FOR K = 1 TO 3 
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2430 PRINT DO’READ r SAA$< 0 )?" .R’ ii J + K) 

2440 INPUT ACil.K) 

2450 NEXT K 
2460 NEXT I 

2470 REN READ SUPPRESSIONS 
2480 A = (AEL * 9) + (AD2. » 3 > + 10 
24 90 FOR I = 1 TO AKi 
2500 J = A + <<I-l)*7i 
2510 FOR K = 0 TO 6 

2520 PRINT MJ’READ c f AA$( C )! " ,R‘ ii J r 1 + K5 
2530 INPUT AGXi I j!. ) 

2540 NEXT. K 
2550 NEXT I 

2560 REM READ MAT. i SECT. PROP. 

2570 A = A + (AHL * 7 5 
2580 FOR 1=1 TO AF! 

2590 J=A+((I-i>*10) 

2600 FOF K = 0 TO 9 

2610 PRINT Mr READ ' »AA$( C i J'' ,R‘ 5i J t 1 -r K) 

262C INPUT AKi I»K) 

2630 NEXT K 
2640 NEXT I 
2650 REM READ LOADS 
2660 A = A + iAF2. * 1C) 

2670 FOR I = 1 TO AIL 
2680 J » A + «I - 1) * 75 
2690 FOR K = 0 TO o 

2700 PRINT Mi "READ '• ,AAJ( C ,»’■ »R‘ Pi J + 1 t K 5 
2710 INPUT AJCI.iU 
2720 NEXT 
2730 NEXT I 

2740 PRINT M* r CLOSE’5AA$i 0 5 
2750 RETURN 
2760 REM 

2770 REM READ HEADING 
2780 P$ = CHR* (4; 

2790 PRINT IH»''OPEN HEADING. L20' 

2800 FOR I = 0 TO 5 

2810 PRINT Mr READ HEADING,R,I 

2820 INPUT AA«( I ) 

2830 NEXT 

2840 PRINT M."CLOSE HEADING' 

2850 RETURN 

2860 REM mmmm**m*mmmUtt***»***»»t* *m**»*M*m*»»» 
2870 REM BEAM STIFFNESS MATRIX 
2880 BCi 1.1) = (BE * BA) / NL 
2890 BC( 7.7) = BCi 1.1) 

2900 BCi 7.1) = - BCil.l) 

2910 BCi 2.2) = < 12 * BE * IZ) / i NL t 3) 

2920 BCi 8.8) = BCi 2.2) 

2930 BCi 8.2) = - BCi 2.2/ 

2940 BCi 3.3) = (12 * BE » IY) / iNL T 3 5 
2950 BCi 9.9) = BCi 3,3) 

2960 BCi 9,3) = - BCi 3,3) 

29 70 BCi 4,4) = (BG * BJ) / NL 
2980 BCi 10,10 ) = BCi 4,4) 

2990 BCi 10,4) = - BCi 4,4) 

3000 BCi 5,5) = (4 * BE * 1Y) / NL 
3010 BCi 11,11) = BCi 5,5) 

3020 BCi 6*6) = (4 * BE * IZ) / NL 
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5030 BCv 12*12) = BC<6.6) 

304 0 BCC5.3) = - <6 * BE * IY) / <NL t 2) 

3050 BCill»9> = - BC(5.3> 

3060 BC< 9.5 )= - BC(5»3) 

3070 BC( 11.3) = BC<5.3) 

3080 BCi6.2> = (6 * BE * IZ) / (NL 1 2) 

309 0 BC< 8*6) = - BCC6.2) 

3100 BCC12.8) = - BC( 6?2 ) 

3110 BC( 12.2 ) = BCi6.2> 

3120 BCU1.5) = BC<5.5> / 2 
3130 BC(12r6) = BC<6.6> / 2 
3140 BC<3i5) = BC( 5.3) 

3142 BCC2.6) = BC( 6.2) 

3144 BC( 1.7 ) = BC(7il ■ 

3146 BC(2f8i = BC( 8*2) 

3148 BCi 3*9 ) = BC<9.3 5 
3150 BC^6.6) = BCi 8»6 i 
3152 BCi 5*9 ) = BCi 9,5) 

3154 BCi 4.10) = BCaC.4) 

3156 BCi 3.11 ) = BC> 11.3) 

3156 BC< 2»12 ; = BC< 12.2) 

3160 BCi 5.11 i = BC( 11.5) 

3162 BCi 6.12 ) = BC< 12.6) 

3164 BCi 8.12 ) = BC412.8) 

3166 BCS9.1D = BC< 11.9 ) 

3210 PRINT I PRINT : LOCAL STIFF. OF ELMT *}BR£S' IS COMPLETE!* 1 
3220 RETURN 

3230 rem »mm»*m*»m***mm*»4»mm***m*»*t 
3240 REM SUBROUTINE CLEAR MATRICES 
3250 PRINT "CLEARING BB.BC.BF.BH' 

3260 FOR K1 = 1 TO 12 

3270 FOR K2 = 1 TO II 

3280 BBi Kl.K'2) = 0 

3290 BCi K1.K2 ) = 0 

3300 BFCK1.K2) = 0 
331C BWK1.K2) = 0 
3320 BT>. K1 >K2 ) = C 
3330 NEXT K2 
3340 NEXT K1 
3350 FOR I = 1 TO 3 
3360 FOR J = 1 TO 3 
3370 BNU.J) = 0 
3380 BCU.J: = 0 
3390 NEXT J 

3400 NEXT I 

3410 RETURN 

3420 rem mm*m*»*m»m***mm****m*m»»4*iu* 

3430 REM MATRIX MULTIPLICATION 

3440 REM C=AB ( BOTH MATRICES ARE 12*12) 

3450 FOR K1 = 1 TO 12 
3460 FOR K2 = 1 TO 12 
34 70 FOR K3 = 1 TO 12 

3480 C( K1.K2) = CCK1.K2) + AiKl.KS) * B(K3»K2) 

34 90 NEXT K3 

3500 NEXT K2 

3510 NEXT K1 

3520 rem 

3530 REM SECT. PROF. FOR THE CROUPS 

3540 IF (ABKI8.2) = NG/.) GOTO 3730 
3550 NG2 = A BIX 18.2 ) 
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3560 GOSUB 3750: REM CLEAR SECT. PROF. 

3570 NTS = AKi'NGi.C i 
358C BA = AK> NG/.f 1 • 

3590 BE = AK( NGf.»4 ) 

3600 IF ( ABSi I8t1 - = 1 ) GOTO 3730 
3610 IY = AK<NG?.»2 ) 

3620 IZ = AK( NGSf 3 ) 

3630 BB = AK< NG£ 1 5 ) 

3640 BJ = AK< NGSf 6 ) 

3650 BX = AK( NG£ * 7 ) 

3660 BY * AK<NG£r3) 

3670 BZ =J*KiNC£.9, 

3680 BG » BE / i 2 * (1 + BD)) 

3730 RETURN 

3740 REM ********** ******* **************** ****** ****** 

3750 REM CLEAR SECTION PROPERTIES 

3760 NTS = 0 

3770 BA = 0 

3780 BE = 0 

37 5 0 IY * 0 

3800 IZ = 0 

3810 BIi = 0 

3820 BJ = 0 

3830 BG = 0 

3840 BX = 0 

3850 BY = 0 

3660 BZ = 0 

3870 BG = 0 

3930 RETURN 

4002 rem mmmtmmmm*»«««*M*****mM»»**»*»* 

4004 REM CAL ANGLE OF ROTATION GAMMA ABOUT X AXIS 

4005 FGF; I = 1 TC 3 

4006 FOR J = I TC 3 

4008 BBU.j; = 0 

4009 NEXT J 

4010 NEXT I 

4012 CG = C!SG = 0 
4014 FOR I = 1 TO 3 

4ci6 dx-:d = ctEx-: i > = c 

4018 NEXT I 

4021 IF BIS I > 1 GOTO 4030: REM DC FOR VERTICAL BEAMS 

4022 GOSUB 5170: REM ETM FOR VERTICAL TRUSS 
4028 GOTC 4040 

4030 GOSUB 4810: REM ROT MATRIX R-A4R-B 

4040 REM CALCULATE X-FS 

40 50 HX(1> = BX - XI 

4060 DX(2) = BY - Y1 

4070 DX(3) = BZ - Z1 

4075 IF (BIS = 1) GOTC 4182: REM VERTICAL BEAM 

4080 REM CALCULATE X-PG 

4086 REM CG AND SG FOR NON-VERTICAL BEAMS 

4050 REM MATRIX MULTIPLICATION X-PG=BC*X-FS 

4100 FOR K1 = 1 TO 3 

411C FGF. K2 = 1 TO 3 

4120 EX<K1) = EX-K1) + BBC K1 ifC2) * DXU'2 ) 

4130 NEXT K2 

414C NEXT Kl 

4150 REM CALCULATE COS AND SIN OF GAMMA 
4160 EXC 0) = C EXC 2 > t 2 r EX( 3 ) t 2) t .5 

4170 cg = ex; 2) / Ex;o 
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4180 SG = EX( 3 ) / EX( C ) 

4181 GCTC 4189 

4182 REM GC.5G FOR VERTICAL BEAMS 
4184 BK = (DX<1) t 2 + DXC3) t 2) t .5 
4186 SG = DX( 3) / BK 

4188 CG = - <HX< 1) * CY) / BK 

4189 PRINT I PRINT 'THE FOLLOWING ARE THE DIRECTION COSINES CG.SG ' JCGf SPCc 
3>rSG: PRINT 

419C RETURN 
42CC REN 

4210 REN TRASFORNATIGN MATRIX FOR VERTICAL BEAMS 

4220 BBC 1 »,?) = CY 

4225 BBC 2.1) = - CY * CG 

4230 BB(2.3> = SG 

4235 BBC 3*1) = CY * SG 

4240 BBC3-3) = CG 

4245 GOSUB 4350: REM FILL UP MATRIX 
4250 RETURN 

4345 rem 

4350 REM FILL UP BB MATRIX FOR BEAM 

4360 FOR I = 1 TO 3 

4370 FOR J = 1 TO 3 

438C K = I 4 3»L = J + 3 

4385 M = I + 65N = J + 6 

4390 C = I + 9JP = J 4 9 

4400 BBC K.L ) = BBCI.J) 

4402 BBC M.N ) = BBC I.J) 

4404 BBCO.P) = BBC I.J) 

4410 NEXT J 
44 20 NEXT I 
4500 RETURN 
4510 rem 

4520 REM TRANSFORMATION MATRIX FOR NON-VERTICAL BEAMS 

4530 BK = (CX 1 2 + CZ t 2) t .5 

4540 BBC 1.1) = CX 

4560 BBC 1.2) = CY 

458C BEK 1.3 > = CZ 

4600 BB(2*1) = < ( - CX * CY * CG/ - (CZ » SG;) / BK 
4620 BBC 2.2) = BK * CG 

4640 BBC 2.3) = << - CY * CZ * CG) + (CX * SG >) / BK 

466C BBC 3.1) = < (CX * CY * SG) - C CZ 4CG,) / BK 

4680 BBC 3.2) = - BK * SG 

4700 BBC 3.3) = < C CY * CZ * SG) t (CX * CG)) / BK 
4710 GGSUB 4350: REM FILL UP BB MATRIX 
4790 RETURN 

4800 rem mmmmmm*m»mmm*mm****** 

4810 REM ROTATION MATRIX R-ALPHA4R-BE7A 

4820 BK =■ <CX t 2 + CZ t 2> t .5 

4900 BBC 1.1) = CX 

4910 BBC 1.2) = CY 

4920 BBC 1,3) = CZ 

4930 BBC2.1) = - (CX * CY) / BK 

4940 BBC 2.2) = BK 

4950 BBC 2.3) = - (CY » CZ) / BK 

4960 BBC 3.1) = - CZ / BK 

4970 BBC 3.2) = 0 

4980 BB< 3.3) = CX / BK 

4990 RETURN 

50oc rem mmmmu**mmmmm**m*m***m» 

5010 REM TRANSFORMATION MATRIX FOR NON-VERTICAL SPACE TRUSS 
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5020 GOSUF 4310J REM RCT MATRIX R-AXR-B 

5030 FOR I = 1 TO 3 

5040 FOR J = 1 TC 3 

5100 K = I + 6 

5110 L = J + 6 

5120 BB< K*L ) = BBCI.J) 

5130 NEXT J 

5140 NEXT I 

5150 RETURN 

5160 rem mmmmmmm**mmmm***m**m* 

5170 REM TRANSFORMATION MATRIX FOR VERTICAL SPACE TRUSS 

5180 BB<1&2 > = CY 

5190 BB<7»8) = CY 

5200 BBC 2t1) = - CY 

5210 BB(8*7> - - CY 

5220 BB< 3*3) = 1 

5230 BB<9»9) = 1 

5240 RETURN 

5250 rem )»»»»»*»»»«*»*»* 

5260 REM FINO INVERSE CF TRANSFORMATION MATRIX R 

5270 FOR K1 = 1 TC 12 

5280 FOR K2 « 1 TC 12 

5290 PF< K2»K1) = BBvKl»K2) 

5300 NEXT K2 
5310 NEXT til 
5320 RETURN 

5330 rem #»»•** 

5340 REM STIFFNESS MATRIX FOR RDI> 

5350 BCtlfl) = iBE * BA) / NL 
5360 BCi 7»1 ) = -BCilrl) 

5370 BC< 1»7)= -BCilfl) 

5380 BC(7f7/ = BC( 1*1) 

5390 RETURN 

5400 rem *»* 

5410 REM ELEMENT STIFFNESS MATRIX iGLOBAL) 

54 20 PRINT 'CALCULATION GF ELEMT. STIFF MATRIX i GLOBAL ' 

5430 REM BH=BF»BC4Bf 
5435 FOR A = 1 TC 12 STEP 3 
54 36 B = A + 2 

544 0 FOR K1 = A TC B 

5441 FOR K2 = 1 TC 12 

5442 FOR K3 = A TC B 

5443 BT( K1»K2) = BT(K1*K2> t <BF(K1.K3) * BCi K3.K2 >) 

5444 NEXT K3 

5445 NEXT K2 
544 6 NEXT K1 
5447 NEXT A 

5505 FOR A = 1 TO 12 STEP 3 

5506 B = A + 2 

5510 FOR K1 = 1 TC 12 

5520 FOR K2 = A TC B 

5530 FOR K3 = A TC B 

5540 BH(K1»K2> = BH(K1*K2) 9 <BT<K1*K3) » BBi K3»K2)) 

5550 NEXT K3 
5560 NEXT K2 
5570 NEXT til 
5575 NEXT A 
5580 RETURN 
5590 rem 

5600 REM READ STRUCTURE SUPPRESSIONS INTC S3 
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5610 P* = CHR* (4)! REft CTRL D 
5620 PRINT D*J"OPEN '!AAS<0)>% L6" 

5630 A = (AES * ?) + i ADZ * 3) + 10 
5640 FOR I = 1 TO AHS 
5650 J = A + < < I - 1) * 7 > 

5660 K = 0 

5670 PRINT D$j"REAP ">AA*< 0)»,R" K J + 1 + K> 

5680 INPUT N 

5690 FOR K = 1 TC 6 

5700 PRINT D*i"READ ’ f AA*' 0 ) j 1 . R"»( J + 1 + K> 

5710 INPUT K1 

5720 IF > 0 THEN SS(Nl) = 2 

5730 NEXT K 
5740 NEXT I 

5750 PRINT P*»" CLOSE"! AA*'. 0 ) 

5760 RETURN 
5770 ren 

5780 HONE I PRINT "LIST INPUT INFORMATION" *. PRINT 

5980 REN PRINT ELENENTS 

6000 PRINT "ELEMENT CARLS": PRINT 

6010 PRINT * LINE' . 

6020 HTAB :8)t PRINT "ElMT GRCUF * OF r f 
6040 HTAB (27)1 PRINT "NOPE INC NODE INC 
6050 PRINT " **i 

6060 PRINT SPC< 3 >1" TYPE". SPC 4 J SPC 4);’ELM7 : i SPC 2ti 

6070 PRINT "1"; SPC( 5)i"l'i SPC 4) J ” 2" i SPC 5)i'2‘: PRINT 

6090 FOR I = 1 TQ AES 

6100 HTAB (4)1 PRINT IS 

6110 HTAB (9>{ PRINT ABSU.l). 

6120 HTAB (16)1 PRINT ABSi 1.2). 

6130 HTAB (22)1 PRINT ABSi 1.3 >S 

6140 HTAB (26)5 PRINT ABSi 1.4 ). 

6150 HTAB (33): PRINT ABSi 1.5 K 

6160 HTAB (39): PRINT ABZiI.6)S 
6170 HTAB (1 )! PRINT TAB( 6). ABSi 1.7) 

6200 NEXT 

6210 REM PRINT NPCC 
6220 PRINT ! PRINT 

6230 PRINT "NODE #"S SPC 3)!”X-COCRD' i SPC 4H’Y-CCCRIr ; SPC 4);'Z-CCC 
RD" 

6240 PRINT 

6250 FOR I = 1 TO ADS 
6260 PRINT Ij 

6270 HTAB (11)! PRINT ACI.Ill 
6280 HTAB (21): PRINT ACiI.2>» 

6290 HTAB <31 )J PRINT AC(I»3) 

6300 NEXT 

6310 REM PRINT NODAL SUPPRESSIONS 

6320 PRINT ! PRINT 

6330 PRINT "LINE * NODE # r » SPC 3); 

6340 PRINT *X Y 2 "5 

6350 PRINT "RX RY RZ" 

6360 PRINT 

6370 FOR I = 1 TO AHS 

6380 HTAB (3): PRINT IS 

6390 HTAB ( 10)5 PRINT AGSil.O). 

6400 HTAB ( 17): PRINT AGSd.l )j 
6410 HTAB (21): PRINT AGZ( 1,2)1 
6420 HTAB (26)! PRINT AGS(I»3)» 

64 30 HTAB (30): PRINT AGS( 1.4)1 
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6440 HTAB (35): PRINT AGZ(I,5>5 
6450 PRINT SPC< 4)5AGZ(I,6) 

6460 NEXT 

6470 PRINT I PRINT 

6480 REPS PRINT MATERIAL PROPERTIES 

6490 PRINT 'MATERIAL PROPERTIES’ 

6500 PRINT 

6510 FOR I = 1 TO AF1 STEP 3 
6520 II = I 4 1 
6530 12 = I + 2 

6540 PRINT ’LINE * r J: HTAB (8)5 PRINT 1,11,12 
6550 PRJ(£T ’GROUP"»I HTAB (6): PRINT 1,11,12 
6560 PRINT ’ELMT’i: HTAB (8)1 PRINT AK( 1,0 ),AK( 11,0 >,AK( 12,C ) 

6570 PRINT ’AREA’!: HTAB (8)1 PRINT AK! 1,1 ),AK( 11,1 ,,AK( 12,1 ; 

6580 PRINT ’I-Yr:: HTAB (8)1 PRINT AM. 1,2 ),AK( 11,2 ),AK< 12,2 ) 

6590 PRINT r l-ZZ“;: HTAB (8)1 PRINT AK( 1,3), AK( 11,3 i,AK( 12,3) 

6600 PRINT ’E";: HTAB (8)1 PRINT AK( 1,4 ),AM 11,4 ),Ah< 12,4 ) 

6610 PRINT ’V“i: HTAB (8): PRINT AK( 1,5 ),AK( 11,5 >,AK( 12,5 ; 

66 20 PRINT ’J-,*: HTAB (8): PRINT AK( 1,6 >,AK< 11.6 >,AK( 12,6 ) 

6630 PRINT HTAB <£>: PRINT AK( 1,7 ),AK(I1,7 >, AN-; 12,7 , 

6640 PRINT *Y"'5: HTAB (8): PRINT AK>. 1,8), AM. II, 8), AM. 12,8 > 

6650 PRINT *Z"5J HTAB <8)1 PRINT AK( 1,9 i, AM 11,9), AK< 12,9 ) 

6660 PRINT I PRINT 
6670 NEXT 

6680 REM PRINT ’LOADS" 

6690 PRINT : PRINT ‘LOADS" 

6700 PRINT ’LINE * NCDE *’5 SPC( 3)f 
6710 PRINT ”PX PY PZ "5 

6720 PRINT 'MX MY MZ‘ 

6730 FOR I * 1 TO All 

6740 HTAB (3): PRINT 15 

6750 HTAB (11 SI PRINT AJt 1,0 )? 

6760 HTAB (17)1 PRINT AJ(I,1)5 
6770 HTAB (24 it PRINT AJ< 1,2)1 
6780 HTAB (32)! PRINT AJ(I,3)i 
6790 HTAB (40): PRINT AJ( 1,4)5 
6800 PRINT TAB( 12 )»AJ( 1,5)5 
6810 PRINT TAB; 20)5AJ(I,6) 

6820 NEXT 
6850 RETURN 

6860 rem 

6870 REM PRINT TRANSFORMATION MATRIX 

6880 PRINT "TRANSFORMATION * "5BRI 

6890 FOR I s 1 TO 12 

6900 FOR J = 1 TO 12 

6910 PRINT BBCI.J),' c 5 

6920 NEXT J 

6930 NEXT I 

6940 RETURN 

4950 rem 

6960 REM PRINT TRANSFORMATION MATRIX TRANSPOSE 

6970 PRINT ’TRANSFORMATION TRANSPOSE * “5BRi 

6980 FOR I = 1 TO 12 

6990 FOR J = 1 TO 12 

7000 PRINT BF(I,J)5‘ ’5 

7010 NEXT J 

7020 NEXT I 

7030 RETURN 

7040 rem m*mmm**m*»* m*m*m*»*m**»**» 

7050 PRINT "STIFFNESS MATRIX (LOCAL) * "»BRIL 
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7060 FOR I ^ 1 TO 12 
7070 FOR J = 1 TO 12 
7080 PRINT RCi I rJ)i l ’» 

7090 NEXT J 
7100 NEXT I 
7110 RETURN 

7120 ren mmmummmmm***«M*m**M*****» 

7130 PRINT ’STIFF MATRIX (GLOBAL) * ' iBRX 

7140 FOR I = 1 TO 12 

7150 FOR J = 1 TO 12 

7160 PRINT BH(I,J)f r ’i 

7170 NEXT. J 

7180 NEXT I 

7190 RETURN 

7200 ren 

7210 REN REARANGE DOFIU IN * ORDER OF USAGE 

7220 K = Cl REN RECORB/NUMBER OF DOF IN USE 

7230 FOR I = 1 TO 30C 

7240 IF (DOFIU(ItO) = 0) GOTO 7270 

7250 K = K + 1 

7260 DOFIU(K»3) = I 

7270 NEXT I 

7280 CC = Kt REN RANK OF CONSOLIDATED MATRIX K 
7290 RETURN 

7300 ren mmmmmmmmummmmiimm* 

7310 REN SET AND STORE NODAL SUPPRESSIONS 

7320 FOR I = 1 TO AHX 

7330 RX = (AGX(I.O) - 1) * 6 

7340 FOR J * 1 TO 6 

7350 IF AGX(IrJ) = 0 GOTO 7330 

7360 R1X = RX + J 

7370 DCFIU( R1X»2 > = 1 

7380 NEXT J 

7390 NEXT I 

7400 RETURN 

7410 ren mmm«Mm*mmmmm******m**n»m****** 

7420 REN STORE NODAL LCADS IN STRUCTURE LOAD MATRIX 

7430 FOR 1 = 1 TO AIX 

74 40 FOR J = 1 TC 6 

7450 IF AJ(I*J) = 0 GOTO 749C 

7460 RX = (AJ< I»C) - 1) * 6 

7470 R1X = RX t J 

7480 I»OFIU( RlXf 1) = AJ(I,J) 

74 90 NEXT J 

7500 NEXT I 

7510 RETURN 

7520 ren m*Mmmmnmmwmm*m**m»m****m»***m*»* 

* 

7660 REN SAME ETN (PER ELEMENT) 

7670 REN ONLY TCP-LEFT 3X3 MATRIX IS STORED 

7680 P$ = CHR« ( 4 ) 

7690 B$ = STR< (BRX ) 

7700 A$ = *ETN“ + B* 

7710 R2X = 0 

7718 PRINT Mi”OPEN r iA< 

7719 PRINT Dti”DELETE ’iA$ 

7720 PRINT D<i”OPEN "iA$i“. L20* 

7730 FOR I = 1 TO 3 

77 40 FOR J = 1 TO 3 

77 65 R2X = R2X + 1 




138 


7770 PRINT Dt?"WRITE ■ i At?»R" ?R2i 
7780 PRINT 
7790 NEXT J 

7800 NEXT I 

7810 PRINT Bt?"CLOSE 1 ?At 
7820 RETURN 
7830 ren 

7840 REN CLEAR SSN 

7850 Dt = CHRt (4) 

7880 At = “SSM1X" 

7890 PRINT Dti“OPEN “!At 

7900 PRINJ Dt?'DELETE '?At 

7910 PRINT *SSM HAVE BEEN CLEARED 1 ; PRINT 

7912 REN CREATE A NULL SSN 

7913 PRINT “CREATING A NULi- SSN 1 ’: PRINT 
7916 K3 = 0 

7918 K1 = Cl * k MX t 6! 

7919 PRINT K1»CG,AD7. 

7920 PRINT Dt?"OPEN r ?At?"r L20 11 
7924 FOR I = 1 TO Kl 

7927 PRINT Dt?"WRITE *?At?",R"?I 

7928 PRINT K3 

7930 NEXT I 

7931 PRINT Dt?”CLOSE “?At 

7932 PRINT “LAST ADDRESS OF SSN IS ’SKI 

7933 PRINT i PRINT “NULL SSN HAS BEEN CREATED 1 t PRINT 

7939 RETURN 

7940 ren mmmrnmmmm***m****m*mm* 

7950 REN ERROR MESSAGES 

7960 PRINT "R>. < 0 IN ASSY OF SSN' : GOTO 7930 

7970 PRINT “NUMBER OF NODES EXCEED 50 NODE uIMI'T/# OF DOF EXCEED 300 LIM1 
T“! GOTO 79S0 
7980 STOP 

7990 ren 

805C REN 

8060 REM PRINT DOFIU AND B'J 
8070 PRINT 

8080 PRINT “CHECK PRINTOUT OF DOFIU AND BU 1 

8090 FOR I * 1 TO BP 

8100 PRINT I? SPC< 2)5 

8110 FOR J = 0 TO 3 

8120 PRINT DGFIU(I.J)? SPC< 2)? 

8130 NEXT J 
8140 PRINT 
8150 NEXT I 
8160 RETURN 
8170 rem 

8190 REN STORE ESN ROD (LOCAL) 

8210 Dt = CHRt (4 ) 

8220 Bt = STRt (BRL ) 

8230 At = "ESM 1 ' + Bt 

8245 PRINT DtfOPEN "?At 

8246 PRINT Dt5"DELETE “ ?At 

8250 PRINT Dt?“OPEN "?At?' ! , L20“ 

8295 R2* = 1 

8300 PRINT Dt?"WRITE “?At?“,R "iRZX 
8310 PRINT BC( 1»1) 

8340 PRINT Dt?"CLOSE "?At 
8350 PRINT “ESN SAVED ON DISK' 

8360 RETURN 
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8370 rem ************************************************ 

8380 REM SAVE D0FIUCBP.2) 

8390 D* = CHR$ <4) 

8400 A* = “DOFIU" 

8410 PRINT B*»"OPEN ”?A* 

8420 PRINT D*f"DELETE ‘iM 
8430 PRINT D$i"OPEN ’ »A*i = » L20" 

8450 PRINT D$i"WRITE *?AS?“,Rl ,: 

8460 PRINT BP! REM RANK OF SSM 
8480 PRINT DSTWRITE “ JA4J 1 »R2 r 
8490 PRINT BR2I REM NUMBER OF ELEMENTS 

8496 PRJhjJ Dt?"WRITE "fA*i ! ,R3* 

8497 PRINT CFt REM BANDWIDTH IN NODES 

8500 FOR I = 1 TO BP 

8510 FOR J = 0 TO 2 

8520 K = (I * 3) + J + 1 
8530 PRINT D* 5 11 WRITE " ? At? “ ,R’?K 
8540 PRINT DQFIUC I»J) 

8550 NEXT J 
B560 NEXT I 

8570 PRINT IK?"CLOSE ”?A$ 

8580 PRINT *DGFIU SAVED IN DISK'I PRINT I PRINT 
8590 RETURN 

8600 rem ******************************** 

8610 REM CALCULATE BANDWIDTH 

8620 FOR I = 1 TO AE2 

8625 FOR J = 1 TO AB2CI.3C 

8630 N1 = ABIC 1*41 + (A3SC 1*5 > * CJ - 1)) 

8635 N2 = ABIC 1*6) + (ABICl.7) * (J - 1); 

8640 B = C ABS C N2 - N1) i 
8645 IF B > CF THEN CF = B 
8650 NEXT J 

8655 NEXT I 

8656 Cl = (CF + 1) * 6 

8660 PRINT : PRINT "1/2 BANDWIDTH HAS BEEN CALCULATED AS <. CF+1 i*e* '.CIS PRINT 
8665 K1 = CCI * CADI * 6)) 4 20 

8670 IF K1 > 60000 THEN PRINT "THE CAPACITY OF THE DISK STORAGE HAS BEEN 
EXCEEDED AS C BANDWIBTHtBGF)#2C>6C00C. CHECK TO REDUCE BANDWIDTH OR 
NUMBER OF NODES”: PRINT 
8673 IF K1 > 60000 THEN STOP 

8675 RETURN 

8676 rem *************************************** 

8700 REM CALCULATE ADDRESS OF THE ELEMENT SSM IN DISK 
87 02 IF I < =6 THEN Til = N1 

8703 IF I > 6 THEN T12 = N2 

8704 IF J < -6 THEN T22 = N1 

87 05 IF J > 6 THEN T22 = N2 

87 06 IF I < =6 THEN T32 = I 

87 07 IF I > 6 THEN T32 = CI - 6) 

8708 IF J < =6 THEN T42 = J 

8709 IF J > 6 THEN T42 = (J - 6) 

8710 R32 = (CT12 - i> * 6) + T3%: REM ADDRESS OF THE DIAGONAL 

8715 R42 = (CT22 - 1) * 6) + T421 REM ADDRESS OF TERM C ROW i 

8720 R12 = R32 * Cl 

8725 R22 = R12 - (R32 - R42) 

8730 RETURN 

8735 REM ********************* 

8750 REM SAVE BANDWIDTH CF 
8753 A* = "ID” 
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8755 PRINT IHf"OPEN "iAir'r LI2" 

8760 PRINT Mr"URITE “rAtr"rR2' 

8765 PRINT CF: REN BANDWIDTH 
8770 PRINT Mr "CLOSE "rAi 
8775 RETURN 

8780 ren ******************* 

8800 REN STORE ESN BEAN <LOCAL) 
8810 Di * CHRi ( 4 ) 

8820 Bi = STRi < BRI;) 

8830 At = "ESN" + R$ 

8845 PRINT Mr"OPEK 'rAi 

8846 PRINJ Mr*“DELETE 'rAi 

8850 PRINT Dir"OPEN 'rAiTr L20’ 

8860 FOR I = 1 TO 6 
8865 R2* = I 

88 7 0 PRINT MT'WRITE " iAir 1 »R‘ rR21 
8875 PRINT BC( I r I ) 

8880 NEXT I 

8885 PRINT MT'WRITE ’rASr'rR7" 

8890 PRINT BCC 5r3 ) 

8895 PRINT Dir"WRITE 'rAiT .RS" 

8897 PRINT BC(6r2) 

8900 PRINT Dir"CLOSE "rAi 
8905 PRINT "ESN SAVED ON DISK' 

8910 RETURN 

8915 ren ************************* 
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5 HOME : CLEAR 

10 HTAB i10 >i PRINT "PROGRAM ASTRA2" : PRINT 
20 REH 27 OCT 1980 

90 DIM AAIDt<l),CD< 150 i»CE(150>»BUi 150 ..DOFIU! 150,3 > 

IOC GOSUB 12201 REM READ DOFIU 
110 GOSUB 620: REM READ DATA 

120 GOSUB 1430: REM REARANGE DOFIU IN' ORDER OF USAGE 
130 PRINT "DOFIU HALE BEEN REARANGED' 

135 GOSUB. 61001 REM DELETE DuFIU 

140 GOSUB 4105 REM CLEAR KA.KB.KC.LA.LB 

150 GOSUB 2280: REM CALCULATE U-AA AND u-AB 

155 GOSUB 60001 REM STORE REARANGED DOFIU 

160 GOSUB 2160: REM PRINT DOFIU AND BU 

165 GOSUB 6150: REM SALE BU 

17C GOSUB 2710: REM ASST OF K-AA AND k-Ar 

175 A* = "SSM1X": REM GOSUB 6310: REMCLEAR SSM 

177 IF BO < = 43 THEN GOSUB 70001A$ * "KA1' : GOSUB o3iO;Ai = "LAI': GOSUB 
6310! GOTO 250: REM CAl LA.LB « RAM-CLEAR AA.LA 
180 GOSUB 3750: REM DECOMPOSE K TO FORM l 
185 At = "KAl r J GOSUB 6310: REM ClEAR KA 
190 GOSUB 4770: REM CAlC L-BNLERSE 
195 At = "LAI": GOSUB 6310! REM CLEAR _A 
250 PRINT "RUN ASTRA3 TO COMPLETE THE JOB 
260 END 
270 rem 

410 PRINT 'CLEARING MATRICES KA.KB.rC.i-A.wE 
44C At = "KA1" 

450 PRINT DtrOFEN "iAt 
460 PRINT Dti"DELETE "iAS 
470 At = "KB1" 

480 PRINT Dt!"OPEN ".At 
490 PRINT Dt;"DELETE "iAt 
50C At = "KC1" 

510 PRINT Dti'OFEN "iAt 
520 PRINT Dt*"DELETE "iAt 
530 At = "LAI" 

540 PRINT Dti' OFEN ’ iAt 
550 PRINT Dti"DELETE "iAt 
560 At = "LBl" 

570 PRINT Dti' OPEN ' iAt 
580 PRINT Dti"DELETE "iAt 
600 RETURN 
61C REM 

620 REM SUBROUTINE READ DATA 
630 Dt = CHRt s 4 >: REM CTKi. D 
632 PRINT Dti "OPEN ID. L12' 

634 PRINT Dti'READ ID.Rl 
636 INPUT AAt( 0 > 

638 REM PRINT Dti"READ ID.Si" 

640 REM INPUT CF 
643 PRINT Dti "CLOSE ID" 

650 PRINT : PRINT "READING INPUT INFORMATION" 

660 PRINT Dti'OFEN ' iAAt’.C '/»'» w12' 

670 PRINT Dt."READ "»AAi<,0)i"»R’i1 
680 INPUT AD* 

690 PR I NT Dt i" READ " .AA$( 0 )V >R i2 
700 INPUT All 
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710 PRINT Bt i "READ "JAA*<0>i' ,Ri 3 
720 INPUT API 

730 PRINT Dt5 r READ " iAAt(0)»’ .R 1 »4 
740 INPUT AHI 

75C PRINT DtrREAIi ' iMHOi' ,F.' >5 
760 INPUT All 

765 PRINT Mi'CLOSE r !AA$vO) 

766 PRINT r INPUT DATA HAS BEEN READ' 

860 RETURN 

870 REP m*****«i»*»****»****»»»;»»**JM*f*» 

1220 PRINT 'READ DQFIIK BP-2 1 REP, * ELEMENTS. RANI, CP SSP 
1230 Dt =_ CHR* < 4 ; 

1240 A* = "DOFIU' 

1250 PRINT III!■ OPEN ' iA$!' , L2C' 

1270 PRINT Dti‘READ r »Atr »R1 : 

1280 INPUT BP: REP RANK OP SSP 
1300 PRINT Dt T READ r iA*r.R2' 

1310 INPUT BRi: REP NUPBER CP ELEMENTS 

1314 PRINT Dt•'' READ ' ,R3’ 

1315 INPUT CPI REP BANDWIDTH IN NODES 

1316 CG = <CF t I i »N REn BANDWIDTH IN DOF 

1320 FOR I = 1 TO BF 

1330 FOR -J = C TO 2 

1340 K = <1 * 3> + J t 1 
1350 PRINT D*•"READ ' »A*»',R ; JM 
1360 INPUT DOFIU 1 .' I • J) 

1370 NEXT J 

1360 NEXT I 

1390 PRINT I«;’’CLOSE '!A* 

1400 PRINT 'DOFIU READ PROP DISK' 

1410 RETURN 

1420 REP ****************** **».******* ************************** 
1430 REP REARANOE DOFIU IN * ORDER OF USAGE 
1440 K = 0: REP RECGR.D/NUPBER GF DOF IN USE 
1450 FOR I = 1 TO 150 

1460 IF (DOFIU( I«0; = C) GOTO 1450: REP DGF NGT USED 
1470 K = K + 1 
1480 DOFIU< K»3 ) = I 
1490 NEXT I 

1500 CC = K: REP RANK OF CONSOLIDATED MATRIX K 

1510 RETURN 

1520 rep mm****m*»»i***mm*****»*m»iM'******** 

1530 REP SET AND STORE NODAL SUPPRESSIONS 

1540 FOR I = 1 TO AH2 

1550 RZ = (AGlXIfO) - 1) * o 

1560 FOR J = 1 TO 6 

1570 IF AGZ(I*J) = 0 GOTO 1600 

1580 RlZ = RX + J 

1590 DOFIU! R1Z»2) = 1 

1600 NEXT J 

1610 NEXT I 

1620 RETURN 

1630 rep **************** 

1640 REP STGRE NODAl LOADS IN STRUCTURE LOAD MATRIX 

1650 FOR I = 1 TO All 

1660 RZ = (AJC.C) - 1) * 6 

1670 FOR J = 1 TO 6 

1680 IF AG21I.J) = C GOTO 1710 

1690 RlZ = RZ + J 

1700 DOFIU(R1%»1 ) = AJ(I»J) 
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1710 NEXT J 
1720 NEXT I 
1730 RETURN 

215C rem mmimm *»**********»»**********«****:***»*»»********** 

2160 REM PRINT DOFIl AND BU 
2170 PRINT 

2180 PRINT "CHECK PRINTOUT OF DCFIU AND FU' 

2190 FOR I = 1 TO BP 
2200 PRINT 15 SPC'. 2) 

2210 FOR J = 0 TO 3 
2220 PRINT DOF I ID I»-J > i SF'Ct 2,5 
2230 NEXT. J 
2240 PRINT BUCI ) 

2250 NEXT I 
2260 RETURN 
2270 REM 

2280 REM PARTITIONING MATRICES U-A AND U-B 
2290 BO = 0! REM RANK OF UA 
2300 FOR I = 1 TO CCt REM CONSOLIDATED SSM 
2310 K = D0FIU' If3) 

2320 IF DOFIUtK»2> = 1 GOTO 2350: REM DOF IS SUFPRESSED 
2330 BV = BV + 1 
2340 BIKBV) = K 
2350 NEXT I 

2360 PRINT 'THIS IS BV '' 5BV 
2370 K1 = BV 

2380 FOR I * 1 TO CCI REM CONSOLIDATED SSM 
2390 K = DOFIU( I»3 > 

2400 IF DOFIU< K»2) = 0 GOTO 2430: REM DOF IS FREE 
2410 K1 = K1 + 1 
2420 BLK K1) = K 
2430 NEXT I 

2440 BW = i CC - BV ): REM RANK: OF U-B 

2450 PRINT : PRINT ’U-A AND U-B HAVE BEEN CALCULATED 

24 60 PRINT : PRINT r RANK OF CONSOLIDATED MATRIX ;V IS 5CC 

2470 PRINT I PRINT "RANK OF PARTITIONED MATRIX K-AA IS ' SBL 

2480 PRINT : PRINT ’RANK OF PARTITIONED MATRIX r.'-BD IS 5B- 

2490 PRINT 'MATRICES U-A AND u-B HAVE BEEN CALCULATED’ 

2500 RETURN 

2510 REM ***************************************** 

2550 NEXT J 

271C PRINT 1 PRINT n ASSY OF PARTITIONED SSM K-AA 1 1 PRINT 
2712 REM K-AA IS STORED AS lOL'ER SYMMETRIC MATRIX 
2720 Dt = CHRS (4) 

2725 KlL = 01 REM ADDRESS OF K-AA 
2728 CG = < CF + 1) * 6 

2730 FOP. I = 1 TO BV 

2740 FOR J = 1 TO I 

2770 REM R22 IS THE ADDRESS WHERE THE VAlUE IN SSM IS STORED 
2780 N17. = BU( I ) 

2790 N2Z = BU< J ) 

2795 BG = 0 

2800 IIS: = NIL * CG 

2804 12/. = NIL - N2L 

2805 I3Z =122+1 

2806 REM PRINT r BU>: I >-BU( J ) IS "5121 

2807 IF 132 > CG GGTO 29901 REM OUTSIDE BANDWIDTH 
2810 R2L = I1L - 122. 

2870 A* = " SSMIX’ 

2880 PRINT D*5’OPEN “5AS5’ r L2C" 




2890 PRINT Dt TREAD r !AM', R";R2X 

2900 INPUT BQ 

2910 PRINT It*!' CLOSE ' fAt 

2990 REN WRITE K-AA 

2995 Kl% = K12 + 1 

3000 C* = ’NAT 

3090 PRINT DM 1 OPEN "rCM'i L20‘ 

3100 PRINT DM'URITE "JCM", R' !K12 

3110 PRINT BG 

3120 PRINT D*5”CLOSE "iC* 

3130 NEXT J 

3140 NEXT. I 

3170 PRINT "MATRIX K-AA STORED IN DISK 
3180 REM 

3190 PRINT ‘ASSY OP PARTITIONED SSM K-ABT'I PRINT 
3200 REM K-ABT IS STORED A5 A COMPLETE MATRIX BW*BV 
3210 D* = CHRS (4) 

3220 REM BW IS RANK OF U-B 

3230 REM CC IS THE RANK OF SSM CONSOLIDATED 

3235 Ko = BV + 1 

3240 FOR I = KB TO CC 

3250 FOP. J = 1 TO BV 

3270 REM CALCULATE ADDRESS OF SSM CORRESPONDING TO K-AA 

3280 REM R2X IS THE ADDRESS WHERE THE VAlUE IN SSM WIuL BE 

3290 Nil. = BU< I ) 

3300 N22 = BU< J ) 

3310 IF BU( J ) > BU< I > THEN Nil = BU-:J>1*22 = EL; I , 

3320 REM CALCULATE THE ADDRESS OF K-ABT 
3330 112 = NIX * CC 

3335 122 = N12 - N22 

3336 132 =122+1 

3337 BG = C 

3339 IF 132 > CG GCTG 35301 REM OUTSIDE BANDWIDTH 

3340 R22 = 112 - 122 
34 00 A* = "SSMIX" 

3410 PRINT D*;”CPEN ' J A* ? ’ » L2C" 

3420 PRINT DM'• READ ‘JAti' i R-.R2X 

3430 INPUT BG 

3440 PRINT DM’ CLOSE "»A* 

3445 REM WRITE K-ABT 
3530 C* = "KBl r 

3535 K12 = (U - KS i * BU) + Jl REM ADDRESS OF K-ABT 

3620 PRINT DM "OPEN " »C*? r » L2C 

3630 PRINT D*;"WRITE '5C*S’ » R'rKlX 

3640 PRINT BG 

3650 PRINT DM"CLOSE "fC* 

3670 NEXT J 
368C NEXT I 

3710 PRINT ’MATRIX K-ABT STORED IN DISK' 

3720 RETURN 

3740 REM *********** ***»*.*; ****** ***************».»******** 
3750 PRINT *. PRINT "CALCULATION OF L-AA" 

3790 D* = CHR* (4) 

3800 FOR I = 1 TO BV 
3810 FOR J = 1 TO I 
3820 REM PRINT "CHECK I= r ;iK 
3830 IF I = J GOTO 3990 
3840 XI = Cl REM AXL 
3850 M = J - 1 
3855 IF M = C GOTO 3912 


J=‘ iJ 
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3860 FOR K = 1 TC H 

3862 REtt AXL=AXL+«LA< J»K)*LA( I«K) 

3866 T3 = J!T4 = K 

3868 GOSUB 6400: REtt ADDRESS (J»K) 

3871 AS = -LAI” 

3872 GOSUB 64501 REtt lA(I»K) 

3874 X2 = T8 

3876 T3 = IIT4 = KS GOSUB 64CC: RE.-, ADDRESS (l.Ki 
3878 T22 = T52 

3880 GOSUB 6450: REM LA< K»J ) 

3882 X3 = T8 

3891 XI = JC1 + <X2 * X3) 

3910 NEXT K 

3912 T3 = I5T4 = J 

3914 GOSUB 6400: REtt CAL (I.J) 

3915 AS = ”KA1" 

3916 GOSUB 6450: REtt CAL KhiltJ) 

3917 REtt X2*riA\ I»J) 

3918 X2 = T6 

3920 T3 = JIT4 = J 

3922 GOSUB 64DC: REtt CAL (J»J) 

3923 AS = "LAI” 

3924 GOSUB 6450: REtt CAL LA(J»J) 

3926 X3 = T8I REtt 

39 30 REtt BG=XL(IfJ5 

3940 BO = (X2 - XI i / X3 

3950 REtt CALCULATE ADDRESS a.J) 

3952 T3 = i:T4 = JJ GOSUB 6400 
3955 R22 = T52 

3960 GOSUB 46901 REtt WRITE LA 

3970 GOTO 4140 

3980 REtt I=J 

3990 XI = 0: REtt AXL 

4000 IF I = 1 GOTO 4060 

4010 « = I - 1 

4020 FOR K = 1 TO tt 

4030 T3 = I5T4 = KJ GOSUB 640C: REtt ADDRESS 5I»K) 
4032 AS = "LAI 

4034 GOSUB 6450: REtt READ LA(I.K) 

4036 X2 = TS 

4051 XI = XI + <X2 12) 

4060 NEXT K 

4080 T3 = I*.T4 = 15 GOSUB 6400: REtt ADDRESS iM'. 
4084 AS = “KA1‘ 

4086 GOSUB 6450: REtt READ KAi I • l ) 

4088 X2 = T8 

4090 BQ = (X2 - XI) t .5 
4094 R22 = T52 

4096 GOSUB 46901 REtt WRITE lA<1,I> 

4140 NEXT J 
4150 NEXT I 

4180 PRINT "DECOMPOSITION OP K INTO L IS COMPLETE - 

4190 RETURN 

4680 REtt ********* 

4690 REtt WRITE MATRIX LA 
4695 AS = r LAl” 

4700 PRINT DSi 1 OPEN ■ ;AS » L2C 
4710 PRINT DSi”WRITE -»AS»*t R £ FR2X 
4720 PRINT BQ 
4730 PRINT DSi"CLOSE 1 »AS 
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4740 RETURN 

47so rem ******************************************* 

4770 PRINT t PRINT "CALCULATION OF L-AA INVERSE’ 

4790 Dt = CHRt (4) 

4800 FOR 1=1 TD BV 

4810 REM PRINT "CHECK 1=’ i 15" J=‘ il 

4820 T3 = IJT4 = It GCSUR 64CC: REM ADDRESS il.I) 

4825 A* = "LAI": GOSUB 6450: REM READ LA(I»I) 

4827 XI = T6 

4830 BO = 1 / XI 

4850 GOSUB 5520*. REM WRITE LB 

4860 IF.L = 1 GCTC 5020 

4870 M = I - 1 

4880 FOR J = 1 TC M 

4890 REM PRINT "CHECK I=‘!Ii' J* r iJ 

490C XI = c: REM AXL 

4910 FOR K = J TC M 

4920 T3 = I5T4 = KI GOSUB 6400: REM ADDRESS U.K/ 

4925 A* = "LAI": GOSUB 6450: REM READ LAi I »K) 

4927 X2 = TS 

4930 T3 = K;T4 = j: GCSUB O40C: REM ADDRESS (K.J) 

4935 A* = "LB1” J GCSUB 6450: REM READ LBCK.J) 

4940 X3 = T6 

4941 XI = XI + 1X2 * X3) 

4960 NEXT K 

4970 T3 = i:T4 = I! GOSUB 6400: REM ADDRESS tI»I, 

4975 At = "LAI": GCSUB 6430: REM READ LA( I«I> 

4980 X2 = T6 

4985 BG = - XI / X2 

4990 T3 = I1T4 = j: GCSUB 6400: REM ADDRESS II.J) 

5000 GOSUB 5520: REM WRITE lB 
5010 NEXT J 
5020 NEXT I 

5050 PRINT "CALCULATION OF L INVERTED IS COMPLETE" 

5060 RETURN 

5510 REM ********* 

5520 REM WRITE MATRIX LB 
5525 At = "LB1" 

5530 PRINT Dti"OPEN "iAii'» L2C 
5540 PRINT Dti”WRITE "!A*i L i R"»T52. 

5550 PRINT BG 

5560 PRINT Dti"CLOSE "iAt 

5570 RETURN 

5580 REM ******************************************* 
6000 REM WRITE D0FIWI.3) 

4005 At = "DOFIU" 

6007 PRINT Dti' OPEN "SAti 1 , L1C 
4008 PRINT Dt!"WRITE r iAti' ,R1" 

6009 PRINT BP: REM RANK CF SSM 
4010 PRINT Dti"WRITE " iAS:’ rS2' 

6011 PRINT BR%: REM NUMBER OF ELEMENTS 

6012 PRINT Dti"WRITE • iASi' f R3' 

4013 PRINT CC: REM RANK CF CONSOLIDATED MATRIX 

4015 PRINT Dti"WRITE "iAtiSRC 

6016 PRINT BV! REM RANK OF Ufi 

6016 FOR I = 1 TC BP 

6C20 FOR J = 0 TC 3 

6030 K = (I * 4) + J 

4040 PRINT Dti"WRITE "iAti' ,RiK 

6050 PRINT DCFIUU.J) 
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6060 NEXT J 
6070 NEXT I 

6080 PRINT DtP CLOSE "5 At 
6090 PRINT "DGFIU SAVED IN DISK’ 

6095 RETURN 

6097 re« mtmmmmm 

6100 REM DELETE DOFIU 
6105 PRINT DtPGFEN DOFIU’ 

6107 PRINT DtP DELETE DOFIU’ 

6108 RETURN 

6150 REM WRITE BU 
6160 A* r.-BL” 

6170 PRINT DtPGFEN ’ JA* ?", L4" 

6180 FOR I = 1 TO BF 

6190 PRINT DtP WRITE "f A*i' ,R"51 

6200 PRINT BU< I) 

6210 NEXT I 

6220 PRINT DtP CLOSE r iAt 
6230 PRINT "BU SAVED IN DISK" 

62 4 0 RETURN 

6250 REM m*********»»**»»i 
6310 REM CLEAR SSK.LA 
6340 PRINT It*!' OPEN ' fAt 
6350 PRINT D*i"DE lETE '»At 
6370 RETURN 
639 9 REM 

6400 REM CALCULATE ADDRESS CF KArKBrLAflB 

6410 T51 = 0 

6415 T6 = T3 - 1 

6417 IF To = 0 GCTC 6435 

6420 FOR K6 = 1 TG To 

6425 T5Z = T5J + K3 

6430 NEXT K3 

6435 T57. = T51 + T4 

64 40 RETURN 

64 4 5 rem mm**wm»*i** 

6450 REM READ MATRICES (KAiLA»LB 
64 52 PRINT Dt*"OPEN "fAtP, L2C" 

6455 PRINT Dti"READ "JAtP, R : »T5L 
6460 INPUT TS 
6465 PRINT Dt5"CLOSE ".At 
64 70 RETURN 

6475 REM m*m*»**mm** 

7000 PRINT "CALCULATION QF LA.LB WHEN BV\*aS USING RAM 1 ; PRINT 

7002 IF BV > 48 THEN STOF 

7003 Dt = CHR't ( 4 ! 

7005 REM CALCULATE SIZE OF LA>LB 

7010 Cl = 0 

7015 FOR I * 1 TC BV 

7020 Cl = Cl + I 

7025 NEXT I 

7026 REM ALLOCATE ARRAYS 

7027 DIM KA1( Cl )»LAli Cl )*LB1( Cl) 

7029 rem mmtmmmt* 

7030 REM READ KA 

7031 Ki = 0 
7035 At = "KA1" 

7037 Dt = CHRt (4) 

7040 PRINT DtPOPEN SAtP, L20" 

7045 FOR I * 1 TC BV 
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7050 FOR J = 1 TG I 
7055 K1 = K1 + 1 

7060 PRINT BtfREAD 'iA*i n . R” SKI 
7065 INPUT KA1( Ki ) 

7070 NEXT J 

7075 NEXT I 

7080 PRINT Mi" CLOSE '»A* 

7082 ren 

7085 REN CALCULATE LA 

7090 FOR I * 1 TC BV 

7095 FOR J = 1 TC I 

7100 REN _ PRINT "CHECK I= : !Ir’ J="»J 

7105 IF I = J GOTO 719C 

7110 XI = Cl REN AXl 

7115 N = J - 1 

7120 IF N = 0 GCTC 7150 

7125 FOR K = 1 TC N 

7126 REN AXL=AXL+< LAI' J»K )*LAh KiJ )'/ 

7127 T3 = J!T4 = KI GC3UB 64CC; REN ADDRESS Mitt I 
7130 X2 = LAHT52) 

7133 T3 = I1T4 * K*. GCSUB 64COS REN ADDRESS U»K) 

7135 X3 = LAKT5Z) 

7140 XI = XI 4 as * X3> 

7145 NEXT K 

7150 T3 = JST4 * JS GOSUB 64001 REN CAL <J,J> 

7153 X3 = LA1( T5S ) 

7155 T3 = IJT4 * JS GCSUB 6400! REN ADDRESS U»K) 

7160 X2 = KA1CT5Z) 

7170 LAHT5S) = (X2 - XI ) / X3 

7175 GOTO 7255 

7180 REN Iso- 

7190 XI = OS REN AXL 

7200 IF I = 1 GOTO 7232 

7205 N = I - 1 

7210 FOR K = 1 TG N 

7215 T3 = IST4 = K: GCSUB 64CCS REN ADDRESS (I.K; 

7220 X2 = LAUT57.) 

7225 XI = XI + (X2 t 2) 

7230 NEXT K 

7232 T3 = IST4 * IS GCSUB 64001 REN ADDRESS > I»I> 

7235 X2 = KAKT5Z) 

7250 LA3< T5Z) = (X2 - XI > t .5 
7255 NEXT J 
7260 NEXT I 

7263 GOSUB 7520S REN STGRE LAI 

7265 PRINT "DECOMPOSITION CF K INTO L IS CONPLETE' 

7299 REN mmmmsmmimmmm********** 

7300 PRINT 'CALCULATE LB USING RAN"I PRINT 
7310 FOR I = 1 TO BV 

7315 REN PRINT 'CHECK I=”»14J=' »I 

7320 T3 * IS T4 = IS GOSUB 64CCS REN ADDRESS (I»I) 

7335 LBKT5Z) = 1 / LAKT5Z) 

7340 IF I = 1 GCTC 7415 

7345 M = I - 1 

735C FOR J = 1 TC N 

7355 REN PRINT "CHECK 3=" SIS' J='iJ 

7360 XI = OS REM AXL 

7365 FOR K = J TO M 

7370 T3 * I! T4 = KS GOSUB 6400! REN ADDRESS (I»KJ 
7375 X2 = LAKT57.) 
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7380 T3 = K'.T4 = j: GOSUB 64CC1 REM ADDRESS ; K.J ; 
7385 X3 = LBKTSZ) 

7390 XI = XI + C X2 * X3 i 
7395 NEXT K 

7400 T3 - IIT4 = 11 GOSGB 6400: REM ADDRESS CI»I> 

7401 X2 = LAKT52) 

7403 T3 = i:T4 = J! GOSUB 6400: REM ADDRESS ;1»J> 
7405 LBKT5Z) = - XI / X2 
7410 NEXT J 

7415 NEXT I 

7420 REM STORE LB IN DISK 

7422 At rJTBl" 

74 24 PRINT Dti" OPEN - »A*J' » L2C 
7426 FOR I = 1 TO Cl 
7428 PRINT Dt*"URITE "iAti 1 , R 1 91 
7430 PRINT LBKI) 

74 32 NEXT I 

7434 PRINT Dti"CLOSE "iAt 

7440 PRINT "CALCULATION OF L INVERTED IS COMPLETE' 
7445 RETURN 

7500 rem m*****»m»mm 

7520 REM STCRE LA IN DISK 
7522 At = "LAI" 

7524 PRINT Dti "OPEN * : iAii'» L2C 
7526 FOR I = 1 TO Cl 
7528 PRINT Dti‘’WRITE ' iAti", R* ?I 
7530 PRINT LAKI) 

7532 NEXT I 

7534 PRINT Bti"CLCSE "iAt 
7540 PRINT "LAI STORED IN DISK- 
7545 RETURN 
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10 Rtf’! PRGGF:AM ASTRA3 

20 REM 25 NOVEMBER 1980 

30 HOME 1 CLEAR 
35 DIM AAIDS*.' 5) 

♦0 GOSUB 6201 REM READ DATA 

50 DIM CD(30C),CE<;3C0i,Btf 300),DCFItf 300.3, 

40 GOSUB 122C: REM READ DOFIU 

70 GOSUB 74001 REM READ BU 

115 REH_ .PRINT BP .BRil.CC.BV 

120 IF BV < * 43 THEN GOSUB 30001 GOTO 2031 REM CAlC NC USING RAM 

200 GOSUB 55901 REM CALCULATE LA INVERSE 

203 REM GOSUB 31001 REM CLEAR LB 

210 GOSUB 61601 REM CALCULATE NODAL DEFLECTIONS 

220 GOSUB 67001 REM CALCULATE REACTIONS 

230 GOSUB 72201 REM PRINT REACTIONS 

240 GOSUB 90001 REM STORE NODAl DEFLECTIONS 

250 PRINT " RUN ASTRA4 TO COMPLETE JOB' 

260 END 
270 REM 

620 REM SUBROUTINE READ DATA 
630 Di = CHRS (4)1 REM CTRL D 

640 PRINT D**OPEN ID. L12 r 

641 PRINT P*!’READ ID.ftl' 

642 INPUT AAS< C ) 

643 PRINT D*5"CLOSE II" 

450 PRINT 1 PRINT 'READING INPUT INFORMATION' 

660 PRINT D*»"OPEN '' »AAS< 0 )i» L12' 

670 PRINT PS •'' READ "»A A S v 0 )i" ,R' i i 
680 INPUT AD* 

69C PRINT H*f ’ READ ' 5AASC 0 >? ! ,R- .2 
700 INPUT AEL 

710 PRINT Dtf' READ ' .AA$( 0 : .R »3 

72C INPUT AF* 

730 PRINT Dt. ’ READ ' »AA$\ C )» 1 »R‘ >4 
740 INPUT AM* 

750 PRINT D*»"READ r »AA$( 0 »R‘ 55 

760 INPUT AH. 

800 PRINT P*!"CLOSE JAA«< C ) 

860 RETURN 

870 REM **»*****».*»**»*» »>*»*»» »*»*****>»**!> 

1120 REM READ DGFIU 
1130 D$ = CHRS \4 , 

1140 A* = "DOFIU' 

1150 PRINT DtrCPEN r ;A*r, L1C 
1170 PRINT PSi"READ "iASi' ,Ri : 

1180 INPUT BP1 REM RANK OF S5M 
1190 PRINT D*." READ ' . AS5' .F.2‘ 

1192 INPUT BR*1 REM NUMBER OF ELEMENTS 

1193 PRINT D*S'READ 1 iAST,P.3' 

1194 INPUT CC1 REM RANK OF CONSOLIDATED MATRIX 

1195 PRINT D*«"READ ".AS.'.RO" 

1196 INPUT BV1 REM RANK CF UA 

1197 PRINT B*»"CLOSE ' IAS 

1198 RETURN 

1199 rem *mm*mmm»*mm**o»»*m***iM»*-M 

1220 REM READ DOFItfBP,3 )1REM * ELEMENTS, RANK OF SSM 
1230 GOSUB 1120 





151 


124C PRINT DSi"OPEN r iASi', L1C 

1320 FOR I = 1 TO BP 

1330 FOR J = 0 TO 3 

1340 K = (I * 4) + J 

1350 PRINT Mi r READ ’ iASi',R : iii 

1360 INPUT DOFILK I»J > 

1370 NEXT J 
1380 NEXT I 

1390 PRINT PSi' CLOSE '.A* 

1400 PRINT “DGFIU REAP IN DISK’ 

1410 RETURN 

1420 res. *♦#»***«*#**»*♦ *«* 

2060 PRINT “ERROR-FIND ADDRESS OF KA’; GOTO 2CSC 
2070 PRINT r ERROR-FIND ADDRESS CF KB 1 I GOTO 20SC 
2080 STOP 
2090 REN 

2160 REN PRINT DGFIU AND BU 
2170 PRINT 

2180 PRINT -CHECK PRINTOUT OF DuFIU AND BU' 

2190 FOR I = 1 TO BP 
2200 PRINT Ii SPU 2> 

2210 FGR J = 0 TO 3 
2220 PRINT DOFIUdiJ)! 5FC\ 2 If 
2230 NEXT J 
2240 PRINT BU I ) 

2250 NEXT I 
2260 RETURN 

2270 REN ****** **»****»*]>***>»** **.*>4 ***** *»»>.#;»*»*»>«***»* *** 

3000 REN CLEAR L-AA 
3010 FOR I = 1 TO 3 
3015 BS = SIRS •; I ) 

302 0 AS = ''LA' + B» 

3025 PRINT DSi’OPEN ’ ! AS 
3030 PRINT Mi"DELETE ’ iAS 
3035 NEXT I 
3040 RETURN 
3090 REN 

3100 ren clear l-aa inv. 

3120 AS = ’LBl" 

3125 PRINT DS5 1 OPEN 'iAS 
3130 PRINT Mi "DELETE ' iAS 
3140 RETURN 

3145 ren m*******m»»*m* 

4000 REN FIND ADDRESS CF LB•KC 

4010 T5>. = C 

4015 T6 = T3L - 1 

4017 IF T6 = 0 GC T C 4035 

4020 FOR K8 = 1 TO T6 

4025 T5Z = T5I + KS 

4030 NEXT K3 

4035 T5I = T5S + T45. 

4040 RETURN 

4050 ren *************** 

4100 REN READ MATRIX LB 
4105 DS - CHRS (4) 

4110 PRINT PS?"OPEN "iASS'f L2C r 
4120 PRINT Mi “READ ' iASi' , F.‘ iT5L 
4130 INPUT TS 
4140 PRINT Mi "CLOSE "iAS 
4150 RETURN 
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5580 REM 

5590 PRINT 'CALCULATION OP K-AA INVERSE' 

5605 REM KC< K1»K2 >«KCC Ki .K2 )+< LBT; K1*1C3 )»LBi K3.K2, - 

5607 R2X = C 

5610 FOR K1 = 1 TC PV 

5620 FOR K2 = 1 TC K1 

5630 REM PRINT 'CALCULATE K-A INVERSE v 5f.l» 1 > 1 *K2 i' '/ 

5660 BG = 0 

5650 FOR K3 = 1 TO BV 

5652 IF K1 > K3 GOTC 5760: REM BECAUSE CF 0 TERMS IN ^B 

5654 IF IQ > K3 GOTC 576C1 REM BECAUSE CF G TERMS IN _B 

5660 REIj „ FIND ADDRESS CF i_P TRANSFUSE 

5665 T3Z = K3JT4X = Kit GOSUB 4000: REM ADDRESS LBT=uBi N3>Nl / 

5670 Al = r LBl“: GOSUB 4ICC’. REM READ LBT 
5675 X2 = T6 

5710 REM FIND ABRESS CF LB 

5715 IF K2 = K1 THEN BO = BG + (X2 : 111 GCTC 57oC 
5740 T3X = K3JT41 = K2: GOSUB 4000: REM ADDRESS -B.KSf.-Q; 

574 5 A$ = P LB1''! GOSUB 41CC: REM READ LB 
574 7 X3 = T6 

5750 BQ = BG + < X2 * X3) 

5760 NEXT K3 

5800 R27. = R2‘. + 1! REM ADDRESS KC 
5810 GOSUB 6055: REM WRITE rX 
5820 NEXT K2 
584 0 NEXT K1 

5845 PRINT r K-A INVERSE HAS BEEN CALCULATED' ; PRINT 

5860 RETURN 

5670 REM *********** 

6055 REM URITE KC 

6056 At = "NCI' 

6060 PRINT D»rCFEN ' iAS J' » L20" 

6070 PRINT B*rWRITE SAIP'r R' »R2» 

6080 PRINT BO 

6090 PRINT Dir CLOSE r * AS 

6100 RETURN 

6170 REM **mm***m****m**miM ******* ************** ****** 
6180 PRINT "CALCULATE NODAL DEFLECTIONS' 

6190 REM MARK NCDES SUPPRESSED 
6200 FOR I * 1 TC BF 

6210 IF DDFIUd.2) = 1 THEN DCFIUdrD) = 9999 
6220 HCFIUd«3) = C 
623C NEXT I 

6240 REM STORE LOADS IN MATRIX CE 
6250 FOR I = 1 TO BV 
6260 K = BU-: I) 

6270 CEd) = BGFIUO K»I > 

6280 NEXT I 

6290 REM CALCULATE DEFUECTICN5 

6300 FOR I = 1 TC BV 

6310 BO = 0 

6340 FOB J = 1 TC- BV 

6350 GOSUB 6470: REM READ KC 

6380 NEXT J 

6390 FOR L = 1 TO BV 

64 00 BO = BG + i CIK u) * CE: L ) i 

6410 NEXT L 

6415 K = BU( I ) 

64 20 DOFIU(K.D) = BG 
6430 NEXT I 
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645C RETURN 

646C ren *********** 

6470 REM READ KC 

6472 IF I > * J THEN T3i = IJT41 = J 

6473 IF I < J THEN T3Z = J5T4I = I 

6474 GOSUB 40CCt REM FIND ADDRESS 
6540 A* = 'KC1" 

6550 PRINT DO "OPEN "{fits'. L20" 

6560 PRINT D$?"READ ’{fit?", R"{T5i: 

6570 INPUT CD(J) 

6580 PRINT Dt»"CLGSE "{fit 
6610 REIljRN 

669C REM ************ **************** ***»*.****************♦ 

6700 PRINT "CALCULATE REACTIONS' 

6710 REM STCRE DEFlECTICKS IN MATRIX CE 
6720 FOR I = 1 TO BV 
6730 K = BU< I) 

67 40 CE'; I ) = DQFI0\K.2) 

6750 IF CE\ I> = 9999 THEN CE;I) = C 
6760 NEXT I 

6770 PRINT "DEFLECTIONS HALE BEEN STORED 

6780 REM CALCULATE REACTIONS 

6785 KS = BV + 1 

6790 FOR I = KS TC CC 

680C BG = C 

6810 K = I - BV 

6830 FOR J = 1 TC BV 

6840 GOSUB 6970; REM READ KB 

6670 NEXT J 

6880 for L = 1 TO BV 

6890 BG = BG + (CD; L ) » CE'. L)) 

6900 NEXT L 
6910 K = BU< I S 
6920 D0FIUCK.3 5 = BG 
6930 NEXT I 
6950 RETURN 

6960 REM *********** 

6970 REM READ KB 

6980 R2% = UK » 1) * BV i + J 

7040 fit = "KB1" 

7060 PRINT DtUCFEN "{fit!", L2C 
7070 PRINT Dti’ REfiD "{fit{ ! , R" {R2« 

7080 INPUT CDCJ) 

7090 FRINT DtS"CLOSE ’ {fit 
7110 RETURN 

7210 REM **■*».*** *********** ********** ************** ******************** 

***** 

7220 PRINT "PRINT NODAL tGAD. DEFLECTIONS AND REACTIONS 
7230 K = 0 

7240 PRINT : PRINT "'NODE # ! { SPCv lCJJ'DOF' { SPC; 1C>{"NCDE ..CAD' . SFC-. 1 
0>»"DEFLECTIONS { 3PC\ 10)?' REfiCTICNS": PRINT 
7250 FOR I = 1 TO BF STEF o 
7260 K = K + It REM NODE NUMBER 
7270 FOR J = 1 TC 6 
728C J1 = (K - 1) * 6 + J 
7290 PRINT Ki 
730C HTAB (18 )t PRINT J{ 

7310 HTAB <3C)t PRINT DOFIU(Jl.l). 

7320 IF D0FIU<J1.2> = 9999 THEN HTAB ; 40 it PRINT 3PU 5 )»"**** i 
7330 IF PGFIUv J1.2 ) < > 9999 THEN HTAB <40 )t FRINT SFU 5>«POFIU Jl.2 , 
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7340 HTAB (30! PRINT SPC( 20 >!DCFIU< J1.3 > 

735C NEXT J 
7360 PRINT 
7370 NEXT I 
7380 RETURN 
7399 rem 

740C REM READ BU 
7405 A* = 'BU" 

7410 Dt = **: REM CTRL D 
7415 PRINT B*5"0PEN 'fAtT. L4' 

7420 FOB J = 1 TC BP 

7425 PRINT D*J”REAI> '*A$»'»R"»I 

74 30 INPUT BU( I ) 

74 35 NEXT I 

7440 PRINT Mr CLOSE ' »A* 

7445 PRINT "BU READ FROM DISK" 

7450 RETURN 

7455 REM ***»»***m»*«m*mM *»*>»*** 

8000 PRINT : PRINT ' CALCULATE KC USING RAM' 

8005 Cl = 0: REM SIZE OF ARRAY UB 
8010 FOR I » 1 TC BV 
8015 Cl = Cl + I 
80 20 NEXT I 

8029 DlMKCHCIi 

8030 DIM LBKCI) 

8031 Dt = CHRS 14 i 
8035 REM READ LB 
8C37 At = "LBl" 

8040 PRINT PtrCPEN ‘.At?'. i_2C 
8045 FCR I = 1 TC Cl 
60 50 PRINT DtrifiAD " (At 5" jR" il 
8055 INPUT LB1( I ) 

8060 NEXT I 

8065 PRINT Dt?"CLOSE r iAt 

8070 REM CALCULATION OF KC USING RAM 

807 5 REM KC'. K,1.K2 >=KC>, K1 *K2 M LBT< K1»K3 >*LB; Of K2 / i 

8080 R21I = 0 

8081 PRINT "TEST BV= "JBV 

8085 FCR K1 = 1 TC BV 

8090 FOR K2 = 1 TC Ki 

8095 REM PRINT "CALCULATE K-A INVERSE i"iKi r *"»K2R’ )‘ 

8100 BG = 0 

8105 FGR K3 = 1 TC BV 

8110 IF Kl > K3 GCTC SloCt REM BECAUSE OF 0 TERMS IN LB 

8115 IF K2 > K3 GCTC 81601 REM BECAUSE CF 0 TERMS IN LB 

8120 REM FIND ADDRESS CF LB TRANSPOSE 

8125 T32 = K31T4Z = K15 GCSUB 4CCCI REM ADDRESS LBT=LB\K3» Kl 
8130 X2 = LBKT52) 

8135 REM FIND ADRESS OF LB 

8140 IF K2 = Kl THEN BG = BG i <X2 t 2)1 GCTC S16C 

8145 IF K2 = Kl GCTC 5760 

8150 T32 = K31T4S * K21 GCSUB 40001 REM ADDRESS LBiK3»K2 i 

8155 X3 = LBKT5X) 

8157 BG = BG + < X2 * X3) 

8160 NEXT K3 

8165 R22 = K22 +11 REM ADDRESS KC 
8170 KCKR22) = BG 
8175 NEXT K2 
8180 NEXT Kl 
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8200 PRINT "STORE KOI IN DISK": PRINT 
8205 Dt = CHR* ■: 4 ) 

8210 PRINT D*»" OPEN KC1' 

8215 PRINT H*»”DELETE KCl" 

8220 PRINT D$i ‘ OPEN KC1. L20 1 
8225 FOR I = 1 TO R2i 
8230 PRINT D$*'WRITE KC1. RSI 
8235 PRINT KC1< I) 

8240 NEXT I 

8250 PRINT D*f"CLOSE KOI’ 

8255 PRINT J PRINT "K-A INVERSE HAS BEEN CALCULATED' : PRINT 

8260 REIURN 

8990 ren mmm*4»****m*mmumm***m** 

9000 PRINT "STORE NODAL DEFLECTIONS" 

9001 D* = CHR* v 4 ) 

9003 BP - ADS * 61 REN * OF DCF 

9004 PRINT D*5"OPEN ND" 

9005 PRINT D*DELETE ND‘ 

9006 PRINT D*5"OPEN NIn L20 l 
9010 FOR 1 = 1 TO BP 

9015 PRINT IisrURITE ND. RSI 
902C PRINT DOFIUd.2) 

9030 NEXT I 

9040 PRINT D*i"CLOSE ND" 

9050 RETURN 
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10 HOME I CLEAR l PRINT ’'PROGRAM A57RA4"! PRINT 
20 REM 25 NOVEMBER 1930 
70 PRINT 

90 GOSUB 2801 REM ALlGCATE ARRAYS 

110 GOSUB 620". REM READ DATA 

230 GOSUB 90001 REM READ NODAL DEFLECTIONS 

240 GOSUB 74001 REM CALCULATION OF MEMBER FORCES 

250 HOME : PRINT "END OF JOB" 

260 END 

270 REM. - ****#».M***M **********> *1*1 ******** *»»»«**»».*4M>**M *****»»*♦ 
tttttttttti »******»» *»m **#*#»***»****»»** ♦**»»*»»*********»**#**#»»*» 

280 REM ALLOCATION OF ARRAYS I 
290 DIM AAIDV: 5)»ABEGENiC25.9 ) 

300 DIM NIK 300 ): REM NODAL DEFLECTIONS 

330 DIM BB(12»12>: REM TRANS. MATRIX 

340 DIM BC( 12.12): REM STIFF. MATRIX , ELEMENT-LCCAl , 

350 DIM BF<12): REM ELEMENT FORCES AT THE NODES 

360 DIM BH(12): REM HOLD NODAl DEFLECTIONS 

370 DIM BT-:i2): REM BT=BB*BH USED TO OA.OU.ATE ELEMENT FORCES 

380 RETURN 

390 REM 

620 REM SUBROUTINE READ DATA 
630 D* = CHR* ■: 4 >: REM CTRL D 
640 GOSUB 880 

650 PRINT : PRINT "READING INPUT INFGRMA'IGX' 

660 PRINT Dt."OPEN ” .AA** 0 )*' . 

670 PRINT Dt ."READ " »AA*(0 )»" .R". 1 
680 INPUT ADS 

690 PRINT Dt?"READ "iAAti 0 )i ".R" !2 
70C INPUT AEL 

710 PRINT Dti"READ ' .AAtC 0 " »R‘ .3 

720 INPUT AFL 

730 PRINT Dt; ‘ READ ' .AAtiG)?" »R' .4 
740 INPUT AHL 

750 PRINT Dt.' READ ' iAAti 0 5!" »R‘ 55 
760 INPUT All; 

770 REM READ ELEMENT CARDS 
780 FOR 1=1 TC AE1 
790 J * << I * 1; * 9) + 15 
800 FOR K = 1 TO 9 

810 PRINT Dt."READ " .AAV 0 )i" »R" K J t K ) 

820 INPUT ABIil.N) 

830 NEXT K 
840 NEXT I 

850 PRINT Dt i"CLOSE’fAA*>: C / 

860 RETURN 

870 REM *».***»****>*»»**»* *»»*#*»»** **»*»*»* 

880 REM READ HEADING 
B9C Dt = CHR* (4) 

900 PRINT Dt.' OPEN ID. L12" 

920 PRINT Dti’READ ID. Rl" 

930 INPUT AAti C ) 

950 PRINT Dt."CLOSE ID' 

960 RETURN 

980 REM »**»**«*****»******»****»**»#*»*»»»*»*»»*»**** 

100C REM READ ESM BEAM 
1010 Bt = STRt (BRLi 
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1015 A$ = "ESW r B« 

1020 PRINT D* 5 "OPEN ’iA*!", L2C' 

1030 FOR I = 1 TO 6 
1035 R2Z = I 

1040 PRINT D$5"READ R’ iR2S 

1045 INPUT BCi 1.1) 

1050 NEXT I 

1060 PRINT Dtf’REAB ' » R7' 

1065 INPUT BCi 5*3) 

1070 PRINT B*> 1 READ ' ifii;', RS’ 

1075 INPUT BCi or 2) 

1077 PRINJ D*;'ClCSE iA* 

1080 FOR I = 1 TO o 
1085 J = I + 6 
1090 BC(JrJ) = BCi I * I j 
1095 BCi Jrl) = - BC(M< 

1100 NEXT I 

1105 BCi 11,5) = BCi 5,5) / 2 
1110 BCi 12,6) = BCi o,o) / 2 
1115 BC< 9,5)= - BCi 5,3) 

1120 BO; 8,6) = - BCiofC; 

1125 BCi 11,3) = BCi 5o) 

1130 BCi 12,2) = BCi 6,2.' 

1135 BCi 11,9) = - 30,5,3) 

1140 BCi 12,6.i = - BCi 6-2) 

1145 FOR I =■ 2 TO 12 
1150 K = I - 1 
1155 FOR J = 1 TO K 
1160 BC> J,I) = BC<I,J) 

1165 NEXT J 

1170 NEXT I 

1175 RETURN 

1180 REri MAx*4rm«t *«»*** 

2000 REM CLEAR MATRICES BB,BC 

2010 FOR I = 1 TO 12 

2012 BF ( I= OIBHil . = CtBT(I) = C 

2015 FOR J = 1 TC 12 

2020 BB< I,J) = C:BC I,J) = 0 

2025 NEXT J 

2030 NEXT I 

2035 RETURN 

7390 REN I**************.* > ft **»>*»»*»i»*)r *».»*» 

7400 PRINT ’CALCULATION OF ELEMENT FORCES’: PRINT 
7410 BR2 = Cl REN ELEMENT COUNTER 

7415 PRINT "THE MEMBER FORCES WlLw BE PRINTED CUT IN FREE FORMAT In THE 0 
RDER SHOWN BELOW. THIS IS BECAUSE OF FORMAT LIMITATION WITHIN THE Bn 
SIC PROGRAMMING LANGUAGE’ I PRINT I PRINT 
74 20 PRINT ’ELEMT *’» SPCi 2)1"NODE #■ 1 SPCi 5,i’FX i SFC, i£ >J’ rT i SFCi 

lorpz", spc( ioi 

7430 PRINT *KX"i SPCi 10)?"MY l ; SPCi lC)i"MZ 1 PRINT 

7435 rem mmm*m****m*******»» 

7440 FOR 18 = 1 TO AE5I: REM * OF ELEMT CARDS 

7445 NTZ = ABZi 16,1 ) 

7450 FOR 19 = 1 TO ABZi 18,3)1 REM ELEMT/CARD 
7460 BR2 = BR2 + 1 

7470 N1 = ABXi 18,4 ) + (ABXiIB,5) * i 19 - 1) i 
7480 N2 = ABZi 18,6) + (ABSiI6,7) * <19 - D) 

7485 IF N1 > N2 THEN N3 = N11N1 = K2JN2 = N3 

7487 GOSUB 200C: REM CLEAR MATRICES BB,BC 

7490 IF ABSX 18,1) = 1 THEN GOSUB 35001 REM READ ESM RCI 
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74 95 IF ABZUSfl) = 2 THEN GOSUB 1000: Rtf'. READ ESr, BEAT, 

7500 GOSUB 7600: REN READ ETf! 

7510 GOSUB 794C: REN READ NODAL DEFLECTIONS 

7520 GOSUB 8070: REN CALCULATE FORCES 
7530 GOSUB 8270: REN PRINT ELEMENT LOADS 
7540 NEXT 19 

7560 NEXT 18 

7580 RETURN 

7590 REN m#*********m>'4»*i*»»**:M 
7600 REN LOAD ETfi 

7605 REN ONLY THE TCF 3X3 ETN IS STORED 

7610 B$ =_ STR« ( BRS i 

7615 R22 = 0 

7620 A$ = 'ETN' + BJ 

7622 PRINT DSMCPEN ' iA*«’ r L2C 

7625 FOR I = 1 TO 3 

7630 FOR J = 1 TC 3 

7655 R23i = R2S + 1 

7700 PRINT Dt i 1 READ "iAii’ , R' }R22 

7710 INPUT B&a.J) 

7715 NEXT J 
7720 NEXT I 

7725 PRINT D*>' CLOSE 5 A* 

7730 FOR I = 1 TC 3 

7735 FOR J = 1 TO 3 

7739 11 = I + 3112 = I + 6JI3 =1+9 

7740 J1 = J 4 35J2 = J + 6M3 = J + 9 

7745 BB< I2.J2) = BB(I.Ji 

7750 IF AMC8.1) = 2 THEN BBII1»J1) = BEKI»J)5BB< I3»J3; = BBMM, 
7755 NEXT J 
7760 NEXT I 

7839 RETURN 

7840 REN ***»iM*'M?M..'».»«*»**»*** 

7940 REN READ NODAL DEFLECTIONS 

7950 FOR I = 1 TO 6 

7960 K = <<N1 - 1) * f) + I 

7970 IF ND«K> = 9999 THEN BH< I) = 0: GCTG 7990 

7980 BH< I) = NIK K) 

7990 NEXT I 

8000 FOR I = 7 TC 12 

8010 K = <<N2 - 1 ) * 6) + I - 6 

8020 IF NIK N ) = 9999 THEN BHOI i = 01 GCTC 8040 

8030 BH< I > = NIK K) 

8040 NEXT I 
8050 RETURN 

8060 REN, *****************;'**** i **********X******* 

8070 REN, CALCULATION CF LOADS 
8080 REN BT=BB*BK 

8140 FOR K1 = 1 TC 12 

8150 FOR K3 = 1 TC 12 

8160 BT(K1) = BT< K1 ) + ( BBUClfKS) * BK( K3 ;; 

8170 NEXT K3 

8180 NEXT K1 

8190 REN, BF=BC*?T 

8200 FOR K1 = 1 TC 12 

8210 FOR K3 = 1 TG 12 

8220 BF(K1) = Br< K1) + (BCIK1.K3) * BTOKSi’j 

8230 NEXT K3 

8240 NEXT K1 

8250 RETURN 
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8260 rem mmmmmmmmm***************** 

8270 REM PRINT ELEMENT FORCES 

8280 PRINT BRSJ SPC( 2)5 

8290 FOR I = 1 TO 12 STEP 6 

8300 IF I < =6 THEN PRINT Nl>: GOTO 8320 

8310 PRINT SPC< 4JJN25 

8320 PRINT SPC< 2 >»PF; I)» 

8330 PRINT SPC< 2)5BF(I + 1)5 
8340 PRINT SPC< 2)JBF(I + 2)i 

8350 PRINT SPC< 2)5BF(I + 3)5 

8360 PRINT SPC< 2)SBF(I + 4)5 

8370 PRINJ SPCt 2 )5BF( I + 5; 

839 0 NEXT I 
8395 PRINT 
8400 RETURN 

8410 rem mmmmmmmm**m*mm*m**m»*»**»«*»«* 

8500 REM LOAD ESN FOR ROD 
8510 B$ = STRS ( BR5,) 

8515 At = r ESM' + Bt 

8520 PRINT Uti "OPEN 'PAtr. L2C’ 

8530 PRINT DS."READ ‘.Ati’, RI" 

8540 INPUT BCil.l) 

8550 PRINT lit;’ CLOSE SAt 
8560 BC< 7.7 ) = BC< 1.1) 

8570 BO- 7,1) = - BC* 1.1) 

8580 BC( 1.7)= - BC( 1 * 1 / 

859 0 RETURN 

8595 rem m#**.*m»******>******»*** m*4*»»*********** 

9000 PRINT ’■READ NODAl DEFLECTIONS' 

9005 BP = ADS. * 61 REM * OF DCF 
9C10 Dt = CHRt ' 4 ) 

9015 PRINT IitrCFEN ND» L2C 
9020 FOR I = 1 TO BP 
9030 PRINT Dt»"READ KB. R‘ f I 
9040 INPUT NIK 11 
90 50 NEXT I 

9060 PRINT DS*, 1 CLOSE ND" 

9070 RETURN 

9060 REM ***»*»**m*»m*mf>»mf **»**>»*»* *»*jm 
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EXAMPLE PROBLEMS 

The six example problems solved using ASTRA and 
chosen to verify its accuracy are illustrated in 
Figures 13 through -18. Problems 1 and 2 and 3 represent 
a plane truss, shown in Figure 13, with problem 1 
illustrating the "hand" solution and 2 the corresponding 
computer solution. Problem 3 is the IASTRA run to 
prepare the input data for ASTRA. The structures 
analyzed in example problems 4 through 7 include a plane 
truss, a space truss, a plane frame and a space frame 
and are illustrated in Figures 14 through 18 
respectively. Solutions to these problems were 
calculated using ASTRA and verified with solutions 
obtained previously using other finite element computer 
programs. The figures show the structures to be 
analyzed, including structure dimensions, node and 
element numbering sequence, applied loads and support 
restraints. The computer outputs are listed in the 
following pages and includes the input data, calculated 
results, as well as program comments to let the user 
known how far along the computer solution has progressed. 
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Example Problem 1; Hand Solution 

In this example problem, the truss shown in 
Figure 13a will be analysed using hand calculations and 
the matrix displacement method. In order to construct 
the idealized structure, the truss is divided into the 
three elements shown in Figure 13b joined at three nodes. 

To obtain the structure stiffness matrix [K s ] we 
start with the element stiffness and assume two degrees 
of freedom (DOF) at each node: two translations u and v 
in the X and Y directions. The stiffness matrix for the 
elements can be found in Table III. 



The force matrix [Q] and displacement matrix [S] 
are constructed from Equation (3.24) 


[Q] = [k] [S] 
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The element transformation matrices are 
assembled from Equations (4.42) and (4.43) by deleting 
the degrees of freedom that are not needed. 

For element 1, the direction cosines Cx and Cy 
and element transformation and stiffness matrices are 
calculated as follows: 
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The stiffness matrices [K] for elements 2 and 3 are 
calculated in a similar fashion. 



where B = 1 / (2 \T2) 


The structure stiffness matrix is then calculated from 
Equation (3.27) 
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From Equations (3.27) thru (3.30) the structure 
stiffness matrix is partitioned to form the submatrices 
[K qq ] and [k a/3 ] T . 
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and the unkown displacements and reactions are calculated 
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The element forces are then calculated 
[Q] = [K] [R] [S] 


For element 1 
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1 0-1 Oil 


[q] (1) = 


0 

-1 


0 0 0 
0 1 0 


0 0 0 0 


0 1 
0 0 
0 0 


0 0 
1 0 
0 1 


0 | 
2.0 | 
4.83j 


0.0 | 
2.0 ! 
0.0 I 


element forces for elements 2 and 3 are calculated using 
tne same procedure. 

Thus we obtain the deflections for the structure 
at node 2 


X2=2.0 inches 
Y2=4.8 inches 

and the reactions at nodes 1 and 3 are 
Rxl=-2.0 lb 
Ryl= 0. lb 
Rx3- 1.0 lb 
Ry3=-1. lb 

The element forces are calculated to be: 
2.0 lb in tension for element 1, 

0.0 lb for element 2, and 
1.4 lb for element 3. 



167 


3RUN 

PROGRAM IASTRA 

ANALISING STRUCTURES UITr APPLE 
TKI5 PROGRAM IS USES TC PREPARE THE INPUT DATA 

CHOOSE ONE OF THE AVAILABLE OPTIONS: 

1. CREATE A NE- INPUT BATA FI-E. 

2. LIST INPUT DA'A. 

3. MODIFY AN EXISTING DATA FILE. 

4. COPY OPTION. 

5. EXIT FROM THE PROGRAM. 


EKTER OPTION NUMBER ?1 
SUBROUTINE HEADING 
DUG ID CODE ?PRi 



A MIXTURE OP rEA.*: AND ROD ELEhEf.TS ARE USED. 
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ENTER OF TIG;. NUMBER 7 i 

ENTER NODE NUMBER ANT V GR N To THE QUESTIONS. 

NCGAw SUPPRESSIONS CARL ♦ ... i 

KOBE NUMBER ? 1 

SUFFRE33 X 7 Y 

SUF FREES Y ? Y 

SUr'r.ESS Z ? Y 

ENTER NODE NUMBER AND Y GR N Tu THE GUZSTIGNS. 


NODAL SUPPRESSIONS CARD * ... 2 

NODE NUMBER ? 2 

SUPPRESS X ? N 

SUPPRESS Y 7 K 

SUPPRESS Z ? Y 

ENTER N'GLE NUMBER AND Y GR N 70 


NODAL SUF t RES3LC 
NODE NUMBER T 3 
SUPPRESS X *■ Y 
SUPPRESS r ? Y 
SUPPRESS Z ? Y 
DC YOU WANT TO A 
INPUT MAT. i SEE 
NUMBER OP Gfiuw-3 
ENTER OPTION NUM 
1. ONLY RGB EL 






















1* CREATE A NEw INPUT DATA FILE* 
2. LI5T INPUT DATA* 

3* NOTIFY AN EXISTING DATA FIL.E* 
4* COPY OPTION. 

5* EXIT FRCr, THE PROGRAM* 

ENTER OPTION NUMBER TO 
END OF JOB 
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3RUK ASTRA1 

PRCGRAK ASTRA1 

ANALISING STRUCTURES WITH APPLE 

ENTER DRAWING CODE KUKBER ? RCi 

READING INPUT INFORMATION 
LIST INPUT INFORMATION 

ELEMENT CARDS 

LINE ELKT CROUP # OF NODE INC NODE INC 

1 TYPE ♦ ELKT 112 2 

11 12 10 2 1 

2 1 2 1 2 C 3 0 

NODE # X-CCCRB Y-CCCRI Z-CCCRD 

1 0 C C 

2 1C C C 



LOADS 

LINE # NODE * PX PY PZ HX KY KZ 

1 2 1 1 0 0 0 C 

1/2 BANDWIDTH HAS BEEN CALCULATED AS ;CFtl>*6= IS 

SSK HAVE BEEN CLEARED 


CREATING A NULL SSK 








324 C 

LAST ADDRESS OF SSM IS 324 


3 
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HULL SSM HAS BEEN CREATED 

SIZE OF SSM IS IS 

ELEMENT 4 1 NODE 1= 1 NODE J= 2 

THE DIRECTION COSINES CXfCYrCZ FOR ELEMENT 1 ARE 
1 0 C 

CLEARING BB»BCfBF»BH 
ESM SAVED ON DISK 

CALCULATION OF ELEMT* STIFF MATRIX (GLOBAL) 

ASSY OF STRUCTURE STIFFNESS MATF.IX 
ELEMENT 4 1 STIFFNESS HAS BEEN ADDED TO SSM 

ELEMENT 4 2 NODE 1= i NODE J= 3 

THE DIRlu i ION CCoaKeld uX»utf uz. r»ipi Ei-EMlN* ARE 

C 1 c 

CLEARING BBfBCfBFfBn 
ESM SAVED ON DISK 

CALuLLh t aUN Or c.i_c.M >* 3 i xrr Mp* i p»i.A ^ uLuBAl ) 

ASSY OF STRUCTURE STIFFNESS MATRIX 
ELEMENT 4 2 STIFFNESS HAS BEEN ADDED TO SSM 

ELEMENT 4 3 NODE 1= 2 NODE J* 3 

THE DIRECTION COSINES CX?CYrCZ FOR ELEMENT 3 ARE 
-.7071u6/dl */0/106/oi 0 

CLEARING BBiBCfBFfBH 

ESM SAVED ON DISK 

CALCULATION OF ELEMT* STIFF MATRIX (GLOBAL) 

ASSY OF STRUCTURE STIFFNESS MATRIa 
ELEMENT 4 3 STIFFNESS HAS BEEN ADDED TO SSn 

ASSY OF SSK IS COMPLETE 

SS HAS BEEN STORED IN DGFIU 
F HAVE BEEN STORED IN DGFIU 
DGFIU SAVED IN DISK 


PROGRAM ndlfthl i3 UbMr l-L. j wjJ ♦ 1 U wUfx > XtiOC t-Oh*.' HAD Kvjh r P»ooi\rtM Ad i Rh«. 

3RUNLOAD ASTRAL 
3RUK »- 

PROGRAM Ao t 

READ DOFIIK BPr2) 

DOFXU READ FROM DISK 

RslA.^'.Nw INFu i xNr uRKmT aO«i 
INPUT DATA HAS BEEN READ 
ILF +. u HAVt. Bc.£.N Rc.nRhNuC.i- 
CLEARING MATRICES KA#KBfKCfLAiLB 
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THIS IS BV 2 

U-A AND U-B HAVE BEEN CALCULATED 
RANK OF CCKouLiOATt.^ MATi\aX K IS 9 
RANK OF PARTITIONED MATRIX K-AA IS 2 

RANK OF PARTITIONED MATRIX K-BB IS 7 
MATRICES U-A AND U-B HAVE BEEN CALCULATED 
BOFIU SAVED IN DISK 

CHECK PRINTOUT OF DOFIU AND BU 
110 117 

2 10 12 8 

3 1 C 1 3 1 

4 0 C C 7 2 

5 0 C 0 8 3 

6 C 0 C 5 9 

7 1 1 C 13 13 

8 1 1 0 14 14 

9 1 C 1 15 15 

1C C 0 C C C 

11 C C C C C 

12 C 0 0 0 C 

13 1 0 1 C C 

14 1 C 1 0 0 

15 1 C 1 0 C 

16 C C C C C 

17 C C C C C 

IB 0 C 0 C 0 

BU SAVED IN DISK 

ASSY OF PARTITIONED SSK K-AA 

MATRIX K-AA STORED IN DISK 
ASST OF PARTITIONED SSK K-ABT 

MATRIX K-ABT STORED IN DISK 

CALCULATION OF LA.LB UHEN BV >45 USING RAM 

LAI STORED IN DISK 

DECOMPOSITION OF K INTO L IS COMPLETE 
CALCULATE LB USING RAM 

CALCULATION CP L INVERTED IS COMPLETE 
RUN ASTRA3 TO COMPLETE THE JOB 

3LCAB ASTRA3 
3 RUN 

READING INPUT INFORMATION 
DCFIU READ IN DISK 
BU READ FROM DISK 

CALCULATE KC USING RAM 

TEST EV= 2 

STORE KC1 IN DISK 
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K-A INVERSE HAS BEEN CALCULATED 

CALCULATE NODAL. DEFLECTIONS 

CALCULATE REACTIONS 

DEFLECTIONS HAVE BEEN STORED 

PRINT NODAL. LOAD. DEFLECTIONS AND REACTIONS 


NODE * 


DOF 


3LCAD ASTRA4 
3RUN 

PROGRAM ASTRA4 


NODE LOAD 


Uii 

»u» 

uu 


tut 

uu 

uu 


READING INPUT INFORMATION 
READ NODAL DEFLECTIONS 
CALCULATION OF ELEMENT FORCE 


L'ILL BE PRINTED OUT IN FREE FORMAT IN THE ORDER SHOWN BELOW. 
WITHIN' THE BASIC PROGRAMMING LANGUAGE 


ELEMT i NODE # PX 


MX M'l MI 


1 1 -2 0 0 0 C 0 

2 2 ,0 0 C 0 0 

2 1 0 0 0 0 0 0 

3 0 0 0 C 0 0 


3 2 1.41421356 0 0 0 0 0 
3 -1.<1421356 00000 


END OF JOB 









3 RUN ASTRA1 

PROGRAM ASTRA1 

ANALISING STRUCTURES UITH APPLE 
ENTER DRAWING CODE NUMBER ? RC2 


READING INPUT INFORMATIC 
LIST INPUT INFORMATION 

ELEMENT CARDS 

LINE ELMT CROJF * OF 
* TYPE # ELMT 

11 12 
2 1 2 2 

3 1 3 2 


NODE # X-CDORD Y-CC 

1 C SC 

2 60 8C 

3 0 0 

4 60 0 


LINE * NODE # X Y 

1 1 CO 

2 2 CO 

3 3 11 

4 4 11 


NODE INC NODE INC 
112 2 

12 2 2 
3 111 

12 4-2 


Z-COORD 

0 

0 

0 

0 


Z RX RY RZ 

10 0 0 

10 0 0 

10 0 0 

10 0 0 


MATERIAL PROPERTIES 

LINE * 1 2 

GROUP 1 2 

ELMT 1 1 

AREA 6 8 

I-YY 0 0 

I-ZZ 0 0 

E 1COOO 10000 

V 0 0 

J 0 0 

X' 0 0 

r o c 

Z' o . o 


3 

3 

1 

10 

0 

0 

10000 

0 

0 

0 

0 

0 


LOADS 

LINE * NODE # PX PY PZ MX 

1 2 20 10 0 0 


1/2 BANDWIDTH HAS BEEN CALCULATED AS < CF41>*6= 24 




175 

SSH HAVE BEEN CLEARED 

CREATING A NULL SSH 

57c 0 4 

LAST ADDRESS OF SSH IS 576 

NULL SSH HAS BEEN CREATED 

SIZE OF SSH IS 24 

ELEHENT • 1 NODE 1= 1 NODE J= 2 

THE DIRECTION COSINES CX»Cr*CZ FOR ELEHENT 1 ARE 
,999999999 0 0 

CLEARING BBiBC»BF,Bri 

ESH SAVED ON DI3N 

CALCULATION OF ELEHT. STIFF HATRIX ^GlGBAl) 

ASSY OF STRUCTURE STIFFNESS HATRIX 
ELEHENT 4 1 STIFFNESS HAS BEEN ADDED TO SSH 

ELEHENT * 2 NODE 1= 3 NODE J= 4 

THE DIRECTION COSINES CX»C!»CZ FOR ELEHENT 2 ARE 
.999999999 C C 

CLEARING BBfBCrBFjBri 

ESH SAVED ON DISf, 

CALCULATION OF ELEHT. STIFF HATRIa (GLOBA.) 

ASSY OF STRUCTURE STIFFNESS HATRIX 
ElEHEnT 4 2 STIFFNESS HAS BEEN ADDED TO Safi 

ELEHENT I 3 NODE 1= 1 NODE J= 3 

THE DIRECTION COSINES CX>CY«CZ FOR ELEHENT 3 ARE 

ClEARINO Be f BC >BFiBn 

ESH SAVED ON DISr, 

CALCULATION CF ELEHT. STIFF HATRIX (GlGBA.; 

ASSY OF STRUCTURE STIFFNESS HATRIX 
ELEHENT 4 3 STIFFNESS HAS BEEN ADDED TO SSn 

ELEHENT 4 4 NODE 1= 2 NODE J= 4 

THE DIRECTION COSINES CXiCYrCZ FOR ELEHENT 4 ARE 
0 -1 0 

CLEARING BBfBCiBFfBh 
ESH SAVED.ON DISK 

CALCULATION OF ELEHT. STIFF HATRIX (GLOBA.) 

ASSY OF STRUCTURE STIFFNESS MATRIa 
ELEHENT * 4 STIFFNESS HAS BEEN ADDED TO SSH 

ELEHENT 4 5 NODE 1= 1 NODE J= 4 

THE DIRECTION COSINES CX»CY»CZ FOR ELEMENT 5 ARE 
.6 — • 8 C 

CLEARING BB»BC»BF»BH 
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ESrl SAVED OK DISK 

CALCULATION OF ELEMT. STIFF MATRIX (GLOBAL) 

ASSY OF STRUCTURE STIFFNESS MATRIX 
ELEMENT 4 5 STIFFNESS HAS BEEN ADDED TO SSri 

ELEMENT 4 6 NODE 1= 2 NODE J= 3 

THE DIRECTION COSINES CX»CY.CZ FOR ELEMENT 6 ARE 
-.6 -.8 ~ 0 
CLEARING BBiBCrBF>BH 

ESM SAVED OK DISK 

CALCULATION OF ELEMT. STIFF MATRIX (GLOBAL) 

ASSY OF STRUCTURE STIFFNESS MATRIX 
ELEMENT ♦ o STIFFNESS HAS BEEN ADDED TO SSri 


ASSY OF SSri IS COMPETE 

SS HAS BEEN STORED IN BGFIU 
P HAVE BEEN STORED IN DOFIU 
DOFIu SAVED IN DlSr, 


PROGRAM ASTRA! IS COfiri_E ic.D• TO CONTINUE lQAD AND RUN FROGRAri ASTRA2 

3LQAD ASTRAL 
IRON 

PROGRAM ASTRA2 


READ DGFIO: 


READING INPUT INFORMATION 
INPUT DATA HAS BEEN READ 
DOFIU HAVE BEEN REARANGID 
CLEARING MATRICES KA»KB>KC»LAfLB 
THIS IS BV 4 


li-A AND U-B HAVE BEEN OAi-OULATED 
RANK OF CONSOLIDATED MATRIX K IS 12 

ran;; of partitioned matrix k-aa is 4 

RANK OF PARTITIONED MATRIX K-BB IS S 
MATF.IOES U~A AND U~B HAVE BEEN CALCULATED 
DOrlu SAVED IN D15r 

CHECK PRINTOUT OF DOFIU AND BU 
110 0 11 
2 1 0 0 2 2 

3 10 13 7 

4 0 0 0 7 8 

5 0 0 0 8 3 

o 0 0 0 9 9 

7 1 2C 0 13 13 

S 1 10 0 14 14 

9 1 C 1 15 15 

10 0 0 0 19 19 
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11 0 0 0 20 20 

12 0 0 0 21 21 

13 1 0 1 0 0 

14 10 10 0 

15 1 0 1 0 0 

16 0 0 0 0 0 

17 0 0 0 0 0 

18 0 0 0 0 0 

19 1 0 1 0 0 

20 1 0 1 0 0 

21 1 0 1 0 0 

22 0 0 0 0 0 

23 0 0 0 0 0 

24 0 0 0 0 0 

BU SAVED IN DISK 

ASSY OF PARTITIONED SSn K-AA 

MATRIX K-AA STORED IN DISK 
ASSY OF PARTITIONED SSn K-AB7 

MATRIX K-AB7 STORED IN DISK 

CALCULATION OF LA, LB WHEN Bv >48 USING RAN 

LAI STORED IN DISK 

CALCULATE LB USING RAN 

CALCULATION Or L INVERTED IS COMPLETE 
RUN ASTRA3 TO COMPLETE THE JOE 

3L0AD ASTRA3 
OR UN 

READING INPUT INFORMATION 
DOFIU READ IN BISK 
BU READ FROM DISK 

CALCu-ATE KC UoINu RA* 

TEST Bv= 4 
STORE KC1 IN DISK 


K-A INVERSE HAS BEEN CALCULATED 

CALCULATE NODAL DEFLECTIONS 

DEFLECTIONS HAVE BEEN STORED 

PRINT NODAL LOAD, DEFLECTIONS AND REACTIONS 

REACTIONS 

0 
0 

0 

0 
0 
0 


,0350427345 

.0102564101 

**** 
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2 2 

2 3 

2 4 

2 5 

2 6 

3 1 

3 2 

3 3 

3 4 

3 5 

3 6 

4 1 

4 2 

4 3 

4 4 

4 5 

4 6 

STORE NODAL DEFLECTIONS 
RUN ASTRA4 TO CONFUTE JOB 

3LOAD ASTRAL 
TRUK 

PROGRAM ASTi\n4 


READING INPUT INFORMATION 
READ NODAL DEFLECTIONS 
CALCULATION OF ELEMENT FORCES 

THE MEMBER FORCES WILL BE PRINTED OUT IN FREE FORMAT IN THE ORDER SHOWN BELOW. 
FORMAT LIMITATION WITHIN THE BASIC PROGRAMMING LANGUAGE 


ELEMT # NODE * PX FY F'Z MX MY M2 

1 1 -7.69230771 OOOOO 

2 7.69230771 OOOOO 

23000000 
4 0 0 0 0 0 0 

3 1 -10.2564101 OOOOO 

3 10.2564101 OOOOO 

4 2 6.41025625 OOOOO 

4 -6.41025625 OOOOO 

5 1 12.820512c OOOOO 

4 -12.8205126 OOOOO 

6 2 -20.5125203 OOOOO 

3 20.5125203 OOOOO 


-6.41025625E-03 

mi 


tut 

mi 

mt 


-12.3076922 

-26.6666663 


mt 

mt 

mt 


-7.69230756 

16.6666663 


END OF JOB 
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3LGAD ASTRAi 
3RUN 

PROGRAM ASTRAI 

ANALISING STRUCTURES WITH APPLE 

ENTER DRAWING CODE NUMBER ? RC3 

READING INPUT INFORMATION 
LIST INPUT INFORMATION 

ELEMENT CARDS 

LINE ELKT CROUP * OP NODE INC NODE INC 
* TYPE * ELHT 112 2 

11 13 10 2 1 

2 1 1 2 2 0 3 1 

3 1 113 0 4 0 


NODE 4 X-CQuRD Y-COuftD Z-CuuRD 

10 0 0 

2 75 0 0 

3 0 100 0 

4 0 0 100 


LINE * NODE A X Y Z RX RY RZ 

1 1 1110 0 0 
2 2 0 1 1 0 0 0 

3 3 111000 

4 4 000000 


MATERIAL PROPERTIES 


LINE * 1 2 

GROUf 1 2 

ELNT 1 • 0 

AREA 10 0 

I-YY 0 0 

I-ZI C 0 

E 10000 0 

V 0 0 

v u 0 

X' 0 0 

Y' 0 0 

r • o . o 


5 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


LOADS 

LINE * NODE * PX PY PZ MX MY KZ 

1 4 30 40 0 0 0 0 

1/2 BANDWIDTH HAS BEEN CAuCUk-ATED AS iCFTl)46= 24 







SSh HAVE BEEN CLEARED 
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CREATING A NULL SSH 

576 0 4 

LAST ADDRESS OF SSh IS 576 

NULL SSH HAS BEEN CREATED 

SIZE OF SSH IS 24 

ELEMENT * 1 NODE 1= 1 NODE J= 2 

THE DIRECTION COSINES Cx.CYfCZ FOR ELEMENT 1 ARE 
.999999999 0 C 

CLEARING BBiBCiBF.BH 

ESh SAVED ON DISK 

CALCULATION OF ELEHT. STIFF MATRIX ( GlOBAl. i 
ASSY OF STRUCTURE STIFFNESS MATRIa 
ELEMENT 4 1 STIFFNESS HAS BEEN ADDED TO SSh 

ELEMENT 4 2 NODE 1= 1 NODE J= 3 

THE DIRECTION COSINES CX»Ci»CZ FOR ELEMENT 2 ARE 
0 .999999779 0 

CLEARING BBiBCfBFiBn 

ESM SAVED ON DISK 

CALCULATION Or ELEHT. STIFF MATRIa (GUSBAu; 

ASSY OF STRUCTURE STIFFNESS MATRIX 
ELEMENT 4 2 STIFFNESS HAS BEEN ADDED TO SSn 

ELEMENT * 3 NODE 1= 1 NODE J= 4 

THE DIRECTION COSINES CX.CuCZ FOR ELEMENT 3 ARE 
0 C .99997799? 

CLEARING BBiBCiBF.BH 

ESH SAVED ON DISK 

CALCULATION OF ELEHT. STIFF MATRIX (GLOBAL) 

ASSY OF STRUCTURE STIFFNESS MATRIX 
ELEHENT 4 3 STIFFNESS HAS BEEN ADDED TO SSH 

ELEMENT 4 4 NODE 1= 2 NODE J* 3 

THE DIRECTION COSINES CXiCiiCI FOR ELEHENT 4 ARE 
-.6 .799999997 C 

clearing bb.bc7DF.bh 

ESH SAVED, ON DISK 

CALCULATION OF ELEHT. STIFF MATRIX (GLOBAL) 

ASSY OF STRUCTURE STIFFNESS MATRIX 
ELEHENT 4 4 STIFFNESS HAS BEEN ADDED TO SSH 

ELEHENT 4 5 NCDE 1= 2 NODE J= 4 

THE DIRECTION COSINES CX.CY.CZ FOR ELEHENT 5 ARE 
-.6 0 .799999999 

CLEARING BB>BCiBF >Bn 
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ESM SAVED Ok DISK 

CALCULATION Or EuEMT. STIFF MATRIX (ClOBAl) 
ASSY OF STRUCTURE STIFFNESS MATRIX 
ELEMENT ♦ 5 STIFFNESS HAS BEEN ADDED TO SSn 

ELEMENT # 6 NODE 1= 3 NODE J= 4 


THE DIRECTION COSINES CXfCnCZ FOR ELEMENT 6 ARE 
0 -.707106781 .707106781 

CLEARING BB.BC.BF.BH 

ESM SAVED ON DISK 

CALCULATION OF ELEMT. STIFF MATRIX < GlOBAl. > 

ASSY OF STRUCTURE STIFFNESS MATRIX. 

ELEMENT ♦ 6 STIFFNESS HAS BEEn ADDED 70 SSri 


ASSY OF SSn IS CGnFuETE 


SS HAS BEEN STORED IN DOFIu 
P HAVE BEEN STORED IN DOFIu 
DOFIU SAVED IN DISK 


PROGRAM ASTRA1 IS COMPLETED. TO CONTINUE LOAD AND RUN PROGRAM AS7RA2 

3L3AD ASTRA 22 
DRUM 

PROGRAM ASTRA2 


BGFIU READ FF.Cr, 


READING INFuT INFORMATIGN 
INPUT BATA HAS BEEN READ 
DuFIU KA.u RuAr.ANvuu 




1NED MATRIX. K-Ar 


RANK CF PARTITIONED MA7F 
MATRICES U”A AND U-c riAV 
DOFIu SAVED IN DISK 

CHECK PRINTOUT OF DOFIU 
110 117 

2 1 0 i 2 It 

3 1 0 1 3 20 


7 1 0 C 13 3 

6 1 0 1 14 8 
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11 0 0 0 20 14 

12 0 0 0 21 15 

13 1 0 1 0 0 

14 1 0 1 0 0 

15 1 0 1 0 0 

16 0 0 0 0 0 

17 0 0 0 0 0 

18 0 0 0 0 0 

19 1 30 0 0 0 

20 1 60 0 0 0 

21 1 0 0 0 0 

22 0 0 0 0 0 

23 0 0 0 0 0 

24 0 0 0 0 0 

BU SAVED IN DISK 

ASSY OF PARTITIONED SSri K-AA 

MATRIX K-AA STORED IN DISK 
ASSY OF PARTITIONED SSrt K-ABT 


MATRIX K-AET STORED IN DISK 

CALCULATION OF LA,LB WHEN BV:=46 USING RAN 

LAI STORED IN DISK 

DECOMPOSITION OF K INTO L IS COMPLETE 
CALCULATE LB USING RAN 


CALCULATION OF L INVERTED IS COMPLETE 
RUN ASTRA3 TO COMPLETE THE JOE 

3L0AD ASTRAS 
3RUN 

READING INPUT INFORMATION 
DCFIU READ IN DISK 
BU READ FROM BISK 

CALCULATE KC USING RAM 

TEST BV= 4 

STORE KC1 IN DISK 


K-A INVERSE HAS BEEN CALCULATED 


CALCULATE NODAL DEFLECTIONS 

CALCULATE REACTIONS 

DEFLECTIONS HAVE BEEN STORED 

PRINT NODAL LOAD, DEFLECTIONS AND REACTIONS 


NODE t 


NODE LOAD DEFLECTIONS 


REACTIONS 


1 o tm 

2 o *m 

3 C *U» 

4 0 0 

5 0 0 

6 0 0 


0 

0 

0 


-24.6710526 

0 

-79.9999998 









READING INPUT INFORMATION 
READ NODAL. DEFLECTIONS 
CALCULATION OF ELEMENT FORCES 

THE MEMBER FORCES WILL EE PRINTED OUT IN FREE FORMAT IN THE ORDER SHOWN BELOW. 
FORMAT LIMITATION WITHIN THE BASIC PROGRAMMING LANGUAGE 


ELEMT 4 NODE 4 PX F'Y PZ MX MY MZ 

1 1 -24.6710526 C C 0 0 0 

2 24.6710526 OOOOO 

21000000 

3 0 0 0 0 0 0 

3 1 -79.9999999 00000 

4 79.9999999 00000 

4 2 -8.88157897 00000 

3 8.88157897 00000 

5 2 49.9999993 00000 

4 -49.9999995 00000 

6 3 56.5685425 0 0 0 0 0 

4 -56.5635425 00000 


END OF JOB 
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3L0AD A5TRA1 
□RUN 

PROGRAM A5TRA1 


ANALISING STRUCTURES WITH APPLE 
ENTER DRAWING CODE NUMBER ? RC4 


READING INPUT INFORMATION 
LIST INPUT INFORMATION 

ELEMENT CARDS 

LI« ELKT CROUP * OF NODE INC NODE INC 
* TYPE * ELMT 112 2 



3 

3 

3 


NODE * X-CCGED 


Y-CQORI’ Z-COOAD 


5 


0 C 0 

240 0 0 

0 144 0 

240 144 0 

4SC 144 0 

C 253 0 

240 253 0 

450 255 0 

0 ^ 523 0 


LINE 4 NODE 4 


9 

10 
11 


8 

9 

AO 

11 


1111 
0 111 

1111 
OC11 
0 0 11 

0 0 11 

0 0 11 

0 0 11 

0 0 11 

0 0 11 

0 0 11 


1 0 
1 0 
1 0 
1 0 
1 0 
1 0 
1 0 
1 0 
1 0 
1 0 


MATERIAL PROPERTIES 

LINE 4 12 3 
GROUP 1 2 3 
ELMT 22 0 
AREA 27.7 14.7 0 
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I-YY 124 56.4 0 
I-ZZ 3270 395 0 
£ 29GGC000 S900000C 


V .25 .25 

i 50 5C 

X' 0 -900 

V 900 0 

T 0 0 


0 

0 


LOADS 

LINE « NODE 4 
1 4 

3 6 

4 7 

5 8 

6 9 

7 10 

8 11 


PX PY PZ 

50000 -50GGG 0 0 

0 -100000 0 0 

0 -50000 0 0 

75000 -50000 0 0 

0 -150000 0 0 

0 -10C000 0 c 

100000 -1CCOOO 0 0 

0 -100000 0 0 


NX 


N't NZ 

0 -2000000 

0 0 

0 2000000 

0 -2000000 

0 -2000000 

0 4000000 

0 -4000000 

0 4000000 


1/2 BANBUISTH HAS BEEN CALCULATED 




SSr. HAVE BEEN CLEAREI 
CREATING A NULL 55rl 


LAST ADDRESS Or SSri IS 155* 


Nili-L 3SM HAS BEEN CREATES 


SIZE GT SSr, IS 66 


ELEMENT 4 1 KOBE 


NODE J= 5 




ELEMENT 1 ARE 


THE FOlLCLING ARE THE SI REST 10;. COSINES CGrSG .995555555 0 

CLEARING BBiBCrBRiBH 

LOCAL STIFF. Or ELMT 1 IS CCMFtETES 


ESn SAVES ON SIS;. 

CALCULATION OF ELEMT. STIFF MATRIX ‘.GluBAwj 
ASSY OF STRUCTURE STIFFNESS MATRIX 
ELEMENT 4-1 STIFFNESS HAS BEEN ASSES TO SSr, 

ELEMENT 4 2 NODE 1= 5 NODE J= 6 

THE DIRECTION COSINES CZ»CY*CZ FOR ELEMENT 2 ARE 


THE FGlLGUING ARE THE DIRECTION COSINES CG.SG .999995990 0 

Clearing bb»bc»bfi3h 















186 


LOCAL STIFF. OF EI.HT 2 IS COMPLETED 


ESP. SAVED ON DISK 

CALCULATION OF ELEMT. STIFF MATRIX <GLOBAL > 
ASSY OF STRUCTURE STIFFNESS MATRIX 
ELEMENT * 2 STIFFNESS HAS BEEN ADDED TO S5M 

ELEMENT * 3 NODE 1= 7 NODE J= S 


THE DIRECTION COSINES CX.CY.CZ FOR ELEMENT 3 ARE 
1 0 0 

THE FOLLOWING ARE THE DIRECTION COSINES CG.SG .995999995 0 


CLEARING BBiBC.BF.3H 

LOCAL STIFF. OF ELMT 3 IS COMPLETED 

ESM SAVED ON DISK 

CALCULATION OF ELEMT. STIFF MATRIX (GLOBAL) 

ASSY OF STRUCTURE STIFFNESS MATRIX 
ELEMENT 4 3 STIFFNESS HAS BEEN ADDED TO 3SM 

ELEMENT * 4 NODE 1= S NODE J= 9 

THE DIRECTION COSINES CX.CY.CZ FOR ELEMEfnT 4 ARE 
1 C C 

THE FOLLOWING Ar.E THE DIRECTION COSINES CG.SG .999999995 0 

CLEARING BB.BC.BF.BH 

i-GCnu STIFF. OF ELMT 4 IS C3MFLETED 

ESM SAVED ON DISK 

CALCUkA i iGf. ur Ek.tM.. S. ST MA.F.IX iGbuBAi./ 

ASSY GF STRuGTuRE STIFFNESS MATRIX 
ELEMENT 4 4 STIFFNESS nAS BEEN ADDED TG 5Sr, 

ELEr.c.v. 4 5 Nui’c I s lv NuDic J = 11 

THE DIRECTION COSINES CX.CY.CZ FOR ELEMENT 5 ARE 

1 c c 

THE FOLLOWING ARE THE DIRECTION COSINES CG.SG .999999995 0 


CLEARING BB.BC.BF.Bh 


LOCAL STIFF. OF ELMT 5 IS COMPLETED 
ESM SAVED*ON DISK 

CALCULATION OF ELEMT. STIFF MATRIX i GLOBAL i 
ASSY OF STRUCTURE STIFFNESS MATRIX 
ELEMENT 4 5 STIFFNESS HAS BEEN ADDED TO SSri 

ELEMENT 4 <s NODE 1= 1 NODE J= 4 


THE DIRECTION COSINES CX.CY.CZ FOR ELEMENT 6 ARE 
G 1 C 
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THE FOLLOWING ARE THE DIRECTION COSINES CGfSG .999999995 0 

CLEARING BB.BCfBFfBH 

LOCAL STIFF. OF ELMT 6 IS COMPLETED 

ESM SAVED ON DISK 

CALCULATION OF ELEMT. STIFF MATRIX <GLOBA.., 

ASSY OF STRUCTURE STIFFNESS MATRIX 
ELEMENT • 6 STIFFNESS HAS BEEN ADDED TO SSK 

ELEMENT # 7 NODE 1= 4 NODE J= 7 

THE DIRECTION COSINES GX.CY.CZ FOR ELEMENT 7 ARE 
0 1 C 

THE FOLLOWING ARE THE DIRECTION COSINES CG.SG .999999995 0 

CLEARING BB»BC»BF»Bn 

ESM SAVED ON DISr. 

CALCULATION OF ELEMT, STIFF MATRIX (GlGBAl > 

ASSY OF STRUCTURE STIFFNESS MATRIX 
ELEMENT # 7 STIFFNESS HAS BEEN ADDED TC SSM 

ELEMENT * 5 NCDE 1= 7 NODE J= 10 

THE DIRECTION COSINES CX.CYiCZ FOR ELEMENT S ARE 


THE FOlLOwING ARE THE DIRECTION COSINES CG»5G .999999995 0 

CLEARING 8SfBC»£F>BH 

ESM SAVED CN DCS,. 

CALCULATION OF ELEMT. STIFF MATRIX iGuOSA,./ 

ASSY OF STRUCTURE STIFFNESS MATRIX 
ELEMENT ♦ S STIFFNESS HAS BEEN ADDED TO SSm 

ELEMENT ♦ 9 NODE 1= 2 NODE J= 5 

THE DIRECTION COSINES CXtCY»CZ FOR ELEMENT 9 ARE 


THE FOLLOWING ARE THE DIRECTION COSINES CG*5G .999999995 0 

Clearing bbiBC»efiBh 

LOCAL STIFF. OF ELMT 9 IS COMPLETED 
ESM SAVED On DISK 

CALCULATION CF ELEMT. STIFF MATRIX < GlOBA;.) 

ASSY OF STRUCTURE STIFFNESS MATRIX 
ELEMENT * 9 STIFFNESS HAS BEEN ADDED TO SSM 
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ELEMENT 4 10 NODE 1= 5 NODE J= 8 

THE DIRECTION COSINES CX»CY»CZ FOR ELEMENT 10 ARE 
0 1 0 

THE FOLLOWING ARE THE DIRECTION COSINES CG.SG .999999993 0 

CLEARING BB.BC.BF,Bn 

LOCAL STIFF. OF ELMT 10 IS COMPLETED 

ESM SAVED ON DISK 

CALCULATION OF ELEKT. STIFF MATRIX i GLOBAL) 

ASST OF STRUCTURE STIFFNESS MATRIX 
ELEMENT * 10 STIFFNESS HAS BEEN ADDED TO Sari 

ELEMENT * 11 NODE 1= S NODE J= 11 

THE DIRECTION COSINES CX.CY.CZ FOR ELEMENT 11 ARE 
0 1 0 

THE FOLLOWING ARE THE DIRECTION COSINES CG»SG .999999996 0 

CLEARING BB»BC>BF>Bn 

LOCAL STIFF. OF ELMT 11 IS COMPETED 

ESM SAVED ON DISK 

CALCULATION OF ELEMT. STIFF MATRIX (GLOBAL; 

ASSY OF STRUCTURE STIFFNESS MATRIX 
ELEMENT 4 11 STIFFNESS HAS BEEN ADDED TO 55M 

ELEMENT 4 12 NODE 1= 3 NODE J= 6 

THE DIRECTION COSINES CX»CY»CZ FOR ELEMENT 12 ARE 
C 1 0 

THE FOLLOWING ARE THE DIRECTION COSINES CG»5G .597999999 0 

CLEARING BBiBC»BFfBH 

LGCAt. STIFF. OF ELMT 12 IS COMPLETED 

ESM SAVED ON DISK 

CALCULATION OF ELEMT. STIFF MATRIX ({GLOBAL) 

ASSY CF STRUCTURE STIFFNESS MATP.IX 
ELEMENT 4 12 STIFFNESS HAS BEEN ADDED TO San 

ELEMENT 4 13 NODE 1= 6 NODE J= 5 

THE DIRECTION COSINES CX»CY»CZ FOR ELEMENT 13 ARE 

C 1 0 

THE FOLLOWING ARE THE DIRECTION COSINES CG»SG .999955599 0 

CLEARING BBfBC»BF«Bh 

LGCA_ STIFF. CF ELMT 13 IS COMPLETED 

ESM SAVED ON DISK 
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CALCULATION OF ELEMT. STIFF MATRIX ( GuOBAl ) 
ASSY OF STRUCTURE STIFFNESS MATRIX 
ELEMENT * 13 STIFFNESS HAS BEEN ADDED TO SSfi 

ASSY OF SSK IS CCMPLETE 

SS HAS BEEN STORED IN DOFIU 
f HAVE BEEN STORED IN DOFIU 
DOFIU SAVED IN DISK 


PROGRAM ASTRA1 IS COMPLETED. TO CONTINUE LOAD AND Ruk PROGRAM A5TRA2 

3LOAD ASTRA2 
3RUN 

PROGRAM ASTRAC 

read dcf:u:b?.2 ; 

DOFIU READ FROM DISK 

INF UT DATA HAS BEEN READ 
DCr -HA.z BEiiN KcArvAN»i*‘ 

CtEARI'.j MATF.IOEi KA»;\r»i\C»4-A»LB 
TnlS IS 3. IT 


RANK OF FAPT2T2GKED MATRIX K-AA IS 27 

RANK OF PARTITIONED MATRIX K-BE IS 39 
MATRICES U-A AND U -E HAVE BEEN CALCULATED 
DOFIU SAVED Ik DI5.\ 



3 1 C 1 3 IS 


4 i c : x :9 

5 1 0 1 5 2C 

6 1 0 1 6 24 

7 i 0 0 7 25 

8 1 0 1 8 26 

9 1 C 1 9 30 

1C 1 0 1 1C 31 

11 1 0 1 11 32 

12 1 0 0 12 36 

13 1 0 1 13 37 

14 1 0 rX 14 33 

15 1 0 1 15 42 

16 1 0 1 16 43 

17 1 0 1 17 44 

IS 1 C 0 IS 46 



23 1 C 1 23 56 
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24 1 -2000000 0 24 60 

25 1 C 0 25 61 

26 1 -100000 0 26 62 

27 1 0 1 27 66 

28 1 0 1 28 1 

29 1 0 1 25 2 

3C 1 0 0 3C 3 

31 1 0 0 31 4 

32 1 -500 00 0 32 5 

33 1 0 1 33 6 

34 1 0 1 34 8 

35 1 0 1 35 9 

36 1 2000000 0 36 10 

37 1 75000 C 37 11 

38 1 -50000 0 38 13 

39 1 0 1 39 14 

40 1 C 1 40 15 

41 1 0 1 41 16 

42 1 -2000C0C 0 42 17 

43 1 C 0 43 21 

44 1 -150000 0 44 22 

45 1 0 1 45 23 

46 1 0 1 46 27 

47 1 0 1 47 26 

4S 1 -2CC00CC 0 4E 23 

4? 1 0 0 45 33 










3L0AD AS7RA3 
3RUN 
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READING INPUT INFORMATION 
BDFIU REAI IN DISN 
BU READ FROM DISK 

CALCULATE KC USING RAM 
TEST BV= 27 
STORE KC1 IN DISK 


K-A INVERSE HAS BEEN CALCULATED 

CALCULATE NODAL REFLECTIONS 

DEFLECTIONS HAVE BEEN STORED 

PRINT NODAL uGAD» DEFlxCTICnS AND REAwTIONS 



DEFLECTIONS 


**»» 

*»*i 

*»»» 

*»»» 

»d 


5.13; 

»m 

tut 

ttii 

***» 


2.332201S6E-03 


**»» 

*»*» 

tut 

tut 


-.046995464c 


4.77075564 

1.9611161SE-03 

tut 

ttii 


4.8014944 Zi 
-.161939334 

m* 

utt 


2.33220209E-03 


4.8084506S 

-.0704249636 

**** 


REACTIONS 

-179910.96c 

-5605.72102 

0 

I390ec5~.9 

0 

479555.92D 
0 
0 
0 


-45053.9665 



0 


0 

0 


0 


0 













192 


6 

6 


5 

6 


o tut o 

2900000 -6.1B51267E-C3 0 


7 

7 

7 

7 

7 

7 


2 

3 

4 

5 

6 


75000 4.654BB978 

-50000 -6.82469405E-03 

o utt o 

o *m o 

0 tut 0 

-2000000 -4»03368303E-03 


0 


0 


2 

3 

4 

5 

6 


0 6.42744454 

-150000 -.277798596 

0 ttt* 

0 ttt* 

0 tut 

-2000000 -3.67690258E-03 


0 

0 

0 


0 


0 


0 


2 

3 

4 

5 

6 


0 6.6070372 

-100000 -.120725053 

0 ttt* 

0 ttt* 

0 ttt* 

4000000 1.832754IE-03 


0 

0 

0 


0 


0 


0 


10 

10 

10 

10 

10 

10 


3 

4 

5 

6 


100000 

-100000 

0 

0 

0 

-4000000 


12.7303939 

-.0347633772 


tttt 

ttt* 

-7.87082248E-03 


0 

0 

0 


0 


0 


0 


11 

11 

11 

11 

11 

11 


2 

3 

4 

5 

6 


0 12.7136094 

-100000 -.362456679 

0 tttt 

0 tttt 

0 tttt 

4000000 4.05028969E-03 


0 

0 

0 


0 


0 


0 


STORE NODAL DEFLECTIONS 
RUN ASTRA4 TO COMPLETE JOB 


JLQAIi ASTRA4 
3RUN 

PROGRAM ASTRA4 


READING INPUT INFORMATION 
READ NODAL DEFLECTIONS 
CALCULATION OF ELEMENT FORCES 

THE MEMBER FORCES WILL BE PRINTED OUT IN FREE FORMAT IN THE ORDER SHOWN BELOW. 
FORMAT LIMITATION WITHIN THE BASIC PROGRAMMING LANGUAGE 


ELEMT t NODE * PX PY PZ MX MY HZ 

1 4 -102334.687 -81815.35 000 -15473324.5 
5 102834.687 81815.36 000 -4162366.71 



2 5 -23283.084 -45103.132 000 -20469*6.32 

6 23283.084 45103.132 000 -8777785.38 

3 7 91794.5587 -73616.4274 000 -9765194.18 

B -91794.5587 73616.4274 000 -7902748.44 

4 B 68372.0075 -31146.6633 000 -5914602.7 

9 -88372.0075 31146.6633 000 -1560596.51 

5 1ft 56179.0959 -50373.9138 000 -9170832.11 

U -56179.0959 50373.9138 000 -291B907.24 

6 1 -5805.72101 179910.987 000 13906657.9 

* 5805.72101 -179910.987 000 12000524.2 

7 4 26009.659 27026.3331 000 1472800.25 

7 -26009.659 -27026.3331 000 2416991.76 

8 7 49626.086 43820.9776 000 5346202.41 

10 -49626.086 -43820.9776 000 5170832.14 

9 2 479555.923 -1.2409687E-04 000 -.0361375809 

5 -479555.923 1.2409687E-04 000 9.67025757E-04 

10 5 342S43.674 796C1.6273 000 .6209333.11 

8 -342843.674 -79601.6273 000 5253301.25 

11 8 150373.917 56178.988 000 6564049.86 

11 -150373.917 -56178.983 000 6918907.19 

12 3 226249.795 45038.967 000 -.0222396851 

6 -226249.795 -45088.967 000 6492811.36 

13 6 131146.664 68372.0168 000 4284973.94 

9 -131146.664 -68372.0168 000 5560596.5 


END CF JOB 
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3L0ADRUN ASTRA1 

PROGRAM ASTRA1 

ANALISINC STRUCTURES UITH APPLE 

ENTER DRAWING CODE NUMBER T RC5 

READING INPUT INFORMATION 
LIST INPUT INFORMATION 

ELEMENT CARDS 

LINE ELHT GROUP * OF NODE INC NODE INC 

4 TYPE 4 ELHT 1 12 2 

12 1110 4 0 

2 2 1 1 2 0 4 0 

3 2 1 1 3 0 4 0 


NODE 4 X-COORB Y-COORD Z-COORD 

10 0 0 
2 60 0 120 

3 120 0 0 

4 60 60 60 


LINE 4 NODE 4 X Y Z RX RY RZ 

1 1 111111 

2 2 111111 

3 3 111111 


MATERIAL PROPERTIES 

LINE 412 3 

GROUP 1 2 3 

ELMT 20 0 

AREA 30 0 0 

I-YY 20 0 0 

I-ZZ 30 0 0 

E 29000000 0 0 

V .25 0 0 

J 50 0 0 

X' 0 0 0 

Y' 0 0 0 

V 1000 0 0 


LOADS 

LINE 4 NODE 4 PX PY PZ MX MY MZ 

1 4 10000 0 0 0 0 0 

1/2 BANDUIDTH HAS BEEN CALCULATED AS < CF+1>*6= 24 


SSM HAVE BEEN CLEARED 
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CREATING A NULL SSH 

576 0 4 

LAST ADDRESS OF SSH IS 576 

NULL SSH HAS BEEN CREATED 

SIZE OF SSH IS 24 

ELEttNT * 1 NODE 1= 1 NODE J* 4 

THE DIRECTION COSINES CX.CY.CZ FOR ELEMENT 1 ARE 
>577350269 .577350269 .577350269 

THE FOLLOWING ARE THE DIRECTION COSINES CC.SC -.499999999 .866025403 

CLEARING BB.BC.BF.BH 

LOCAL STIFF. OF ELHT 1 IS COMPLETED 

ESH SAVED ON DISK 

CALCULATION OF ELEHT. STIFF MATRIX (GLOBAL) 

ASSY OF STRUCTURE STIFFNESS MATRIX 
J.EMENT * 1 STIFFNESS HAS BEEN ADDED TO SSH 

ELEMENT I 2 NODE 1= 2 NODE J= 4 

THE DIRECTION COSINES CX.CY.CZ FOR ELEMENT 2 ARE 
0 .707106781 -.707106731 

THE FOLLOWING ARE THE DIRECTION COSINES CC.SG .995383406 -.0959785031 

Clearing bb.bc.bf.bh 

LOCAL STIFF. OF ELHT 2 IS COKPLETED 
ESH SAVED ON DISK 

CALCULATION OF ELEHT. STIFF MATRIX (GLOBAL) 

ASSY OF STRUCTURE STIFFNESS MATRIX 
ELEMENT * 2 STIFFNESS HAS BEEN ADDED TO SSH 

ELEMENT * 3 NODE 1= 3 NODE J- 4 

THE DIRECTION COSINES CX.CY.CZ FOR ELEMENT 3 ARE 
-.577350269 .577350269 .577350269 

THE FOLLOWING ARE THE DIRECTION COSINES CC.SC -.592136907 -.805837378 

CLEARINC BB.BC.BF.BH 

LOCAL STIFF. OF ELHT 3 IS COMPLETED 

ESH SAVED ON DISK 

CALCULATION OF ELEHT. STIFF MATRIX (GLOBAL) 

ASSY OF STRUCTURE STIFFNESS MATRIX 
ELEHENT 4 3 STIFFNESS HAS BEEN ADDED TO SSH 

ASSY OF SSH IS COMPLETE 


SS HAS BEEN STORED IN DOFIU 
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P HAVE BEEN STORED IN DOFIU 
DOFIU SAVED IN DISK 


PROGRAM ASTRA1 IS COMPLETED* TO CONTINUE LOAD AND RUN PROGRAM ASTRA2 

3LQAD ASXRA2 
3RIUN 

PROGRAM ASTRA2 

READ SVIUBP.2) 

DOFIU READ FROM DISK 

READING IW>UT INFORMATION 
INPUT DATA HAS BEEN READ 
DOFIU HAVE BEEN REARANGED 
CLEARING MATRICES KA.KB.KC.LA.LB 
THIS IS BV 6 

U~A AND U-B HAVE BEEN CALCULATED 

RANK OF CONSOLIDATED HATRIX K IS 24 

RANK OF PARTITIONED HATRIX K-AA IS 6 

RANK OF PARTITIONED MATRIX K-BB IS 18 
MATRICES U-A AND U-B HAVE BEEN CALCULATED 
DOFIU SAVED IN DISK 

CHECK PRINTOUT OF DOFIU AND BU 

1 1 0 1 1 19 

2 1 0 1 2 20 

3 1 0 1 3 21 

4 1 0 1 4 22 

5 1 0 1 5 23 

6 1 0 1 6 24 

7 10 17 1 

8 10 18 2 

9 10 19 3 

10 1 0 1 10 4 

11 1 0 1 11 5 

12 1 0 1 12 6 

13 1 0 1 13 7 

14 1 0 1 14 8 

15 1 0 1 15 9 

16 1 0 1 16 10 

17 1 0 1 17 11 

18 1 0 1 18 12 

19 i 10000 0 19 13 

20 1 0 0 20 14 

21 1 0 jO 21 15 

22 1 0 0 22 16 

23 1 0 0 23 17 

24 1 0 0 24 18 

BU SAVED IN DISK 

ASSY OF PARTITIONED SSK K-AA 

MATRIX K-AA STORED IN DISK 
ASSY OF PARTITIONED S3M K-ABT 
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matrix k-abt stored IN SISK 

CALCULATION OF LArLB WHEN 1V<=48 USINC RAM 
LAI STORES IN BISK 

IECOMPOSITION OF K INTO L IS COMPLETE 
CALCULATED! USINC RAM 

CALCULATION OF L INVERTED IS COMPLETE 
RUN A6TRA3 TB COMPLETE THE JO! 

3LOABS ASTRAS 
3RUN 


READING INPUT INFORMATION 
BOFIU READ IN DISK 
BU READ FROM DISK 

CALCULATE KC USINC RAM 
TEST BV= 6 
STORE KC1 IN DISK 


K-A INVERSE HAS BEEN CALCULATED 


CALCULATE NODAL DEFLECTIONS 

CALCULATE REACTIONS 

DEFLECTIONS HAVE BEEN STORED 

PRINT NODAL LOAD, DEFLECTIONS AND REACTIONS 


NODE t 


DOF NODE LOAD DEFLECTIONS REACTIONS 


1 0 m* 

2 o m * 

3 0 tut 

4 0 tut 

5 0 Utt 

6 0 UU 


-4993.81773 

-4992.59791 

-4986.27833 

195.154713 

-248.380568 

176.333181 


2 1 0 

2 2 0 

2 3 0 

2 4 0 

2 5 0 

2 6 0 


-12.9486432 

.231835623 

.465738176 

27.7020922 

367.322795 

559.986885 


3 

3 

3 

3 

3 

3 


1 0 UU 

2 0 UU 

3 0 tttl 

4 0 tttt 

5 0 Ult 

6 0 UU 


-4993.23362 

4992.36608 

4985.81259 

-195.03738 

-239.649604 

165.840146 


2 0 

3 0 

4 0 

5 0 

6 0 


1.7884384E-03 
4.96813565E-08 
2.68683797E-08 
-1 •68373674E-07 
3.14832586E-06 
-2•50380362E-05 


0 

0 

0 

0 

0 

0 




STORE NODAL DEFLECTIONS 
RUN ASTRA4 TO COMPLETE JOB 
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3LOAD ASTRA4 
3RUN 

PROGRAM ASTRA4 


READING INPUT INFORMATION 
READ NODAL DEFLECTIONS 
CALCULATION OF ELEMENT FORCES 

THE ICNBER FORCES UILL BE PRINTED OUT IN FREE FORMAT IN THE ORDER SHOWN BELOW. 
FORMAT LIMITATION WITHIN THE BASIC PROGRAMMING LANGUAGE 


ELEMT * NODE * PX 


PY 


PZ 


MX 


MY 


MZ 


1 1 -8644.48891 

4 8644.48891 


5.65790312 -.862542649 
-5.65790312 .862542649 


71.0760484 165.704804 313.626806 

-71.0760484 -76.0667442 274.359735 


2 


-.165394081 

.165394081 


1.73377339 -12.8415225 -136.234083 650.021027 90.50797e 
-1.73377339 12.8415225 136.234083 439.61829 56.6075752 


3 


8643.74893 -5.60311923 -1.23664234 69,9909736 186.466418 -288.627799 
-8643.74393 5.60311923 1.23664234 -69.9909736 -57.9507782 -293.665432 


END OF JOB 




